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THEORETICAL STUDIES OF SOLAR PUMPED LASERS 


By 

W. L. Harries 
1 . INTRODUCTION 

One concept for collecting solar energy is to use large solar 
collectors and then use lasers as energy converters whose output beams act 
as transmission lines to deliver the energy to a destination. The 
efficiency of the process would be improved if the conversion could be done 
directly using solar pumped lasers, and the purpose of this report is to 
study the possibility of making such lasers. The conversion from laser 
beam energy back to useful work is not discussed. 

There are many applications for such lasers, and these are dealt with 
in Section 2. By including the applications first, the requirements for 
the lasers will be more evident. They are especially applicable to the 
Space program, and include cases where no other methods of delivering power 
seem possible. Using the lasers for conveying information and surveillance 
is also discussed. 

Many difficulties confront the designer of an efficient system for 
power conversion. These involve the nature of the solar spectrum, the 
method of absorbing the energy, the transfer of power into laser beams, and 
finally, the far field patterns of the beams. The requirements of the 
lasers are discussed in Section 3. 

Specific laser configurations are discussed in Section 4. The main 
thrust of this study is into gas laser systems, because for space appli- 
cations, the laser could be large, and also the medium would be uniform and 
not subject to thermal stresses. Dye and solid lasers are treated 
briefly. 
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For gas lasers, a chart of the various possibilities will be shown in 
Section 4, and the various families of gas lasers divided according to the 
mechanisms of absorbing solar radiation and of lasing. Several specific 
models will be analyzed and evaluated according to the criteria of Section 
3. 

Overall conclusions for the program are summarized in Section 5, and 
the performances of the lasers related to the requirements of various 
applications . 


2. APPLICATIONS OF SOLAR PUMPED LASERS 
The main thrust of this study is to examine solar pumped lasers as 
energy converters to supply power. The possibility of providing communi- 
cation and/or surveillance at much lower power levels is also discussed. 

2.1 Solar Pumped Lasers as Power Converters 
The use of very large solar collectors on orbiting space stations with 
lasers sending beams to earth to augment the electrical utilities does not 
seem feasible as the required power levels (many billions of MW) are far 
too high. There are also the problems of directing the beams into popu- 
lated areas and of transmission through the atmosphere. 

The use of solar pumped lasers for space applications however intro- 
duces exciting possibilities as more reasonable power levels are required 
(< 100 KW) , there are unobstructed paths for the beams and atmospheric 
absorption is eliminated in many instances. With planetary atmospheres, 
the laser frequency can be tailored to penetrate an atmosphere such as CO 2 
on Mars, for instance. 

For space applications, we may consider vehicle- to-vehicle power 
transfer whilst in flight, and also orbiting vehicle- to -ground power 
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transfer . In the first case, one vehicle in a mission of several vehicles 
would carry a large solar collector and use a laser to send power to the 
others which would not need to carry such large collectors, or could also 
be in the shadow of a planet. 

The second possibility of an orbiting vehicle sending power to ground 
has the important application of providing power to rovers exploring a 
planetary surface. The use of beam power in this instance is one of the 
few methods possible for prolonged excursions. The NASA Office of 
Aeronautics and Space Technology, and the Office of Exploration sponsored a 
workshop at Langley to consider Beamed Space - Power and our studies on a 
Martian Rover are summarized below. 

2.2 Laser Powered Martian Rover 

The power requirements considered for a Martian Rover were 4.5 KW for 
an unpressurized skeleton vehicle capable of a range of 10 Km and about 
75KW for a vehicle visualized as a mobile habitat with full life support 
systems and unlimited range. The power was to be supplied from orbiting 
satellites with solar collectors a height h above the vehicle. Some 
calculations showing numbers are useful. 

We assume the rover has a receiving dish of diameter D r mounted on the 
vehicle. The transmitter must then be able to point at the receiver within 
an accuracy Atf 

M < D r /h (2-D 

The maximum attainable accuracy in practice is A0 — 0.2" arc * 10 
radian. Thus, if D r - 2m (for the smaller vehicle) and h - 2 x 10^m, 
(geosynchronous orbit for Mars) then A 6 — 10 ^ rad, an impossible value, 
as the present limit of pointing accuracy is 10"^ rad. However, if h were 
reduced to 2 x 10^m, accurate pointing would be possible. 
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Accurate pointing is a necessary but not sufficient condition if the 
rover is in motion as then it must be tracked by the orbiting vehicle. As 
the rover is controlled from the ground, its motion can be regarded as 
random. The orbiting station S at height h above the rover, (Fig. 1) 
senses the position AB, and the sensing signal from AB takes a time t - h/c 
to reach S where c is the velocity of light. The laser beam from S takes a 
similar time, so the total delay between being sensed and receiving the 
laser energy is 2h/c. If the rover moves with velocity v on the ground 
surface, it moves a distance 2hv/c in this time, from AB to A'B'. We 
require BB' « AB - D r or BB' - aD r where a is a precision factor (e.g. 
0.1). Hence we require 

2hv/c < aD r (2-2) 

If v - lOKm/hr - 2.8 ms“^, c-3x lO^rns*^, a - 0.1, D r - 2m then the 
condition is not satisfied at geosynchronous orbit, but is satisfied if h - 
2000 Km. 

Calculations showed that three orbiting satellites 2000 Km symmetri- 
cally above the Martian surface with a period of 3 hours 19 minutes and a 
velocity of 2821 ms"^ could cover 55% of the Martian surface with a view 
time of 56 mins 40s and dead time of 9 mins and 54s, A number greater than 
3 would provide more than adequate coverage. 

To beam the energy at the vehicle so that the majority is collected by 
the receiver dish, the sizes of both transmitting and receiving dishes must 
be greater than certain minimum values. If the diameter of the dishes are 
D t (transmitting), D r , respectively, and A the wavelength of the signal 
then from the condition of diffraction limitation 

D t D r > . ) h (2-3) 

ir 
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where the factor 4/w assumes the best possible case with a Gaussian inten- 
sity profile. With a laser wavelength of 1/im, h - 2 x 10^ m, then D t D r ~ 
2m^ and D t « 1 m which are reasonable values. An important conclusion from 
Eq. 2-3 is that if instead of a laser beam, microwaves with A - 1 cm were 
used, then D t D r — 2 x 10^; for D r — 2m then D^- - lOKm, impossibly large. 

Thus, the orbiting lasers could supply essentially continuous power, 
and no storage requirements are necessary during the Martian nights. 

Powers of up to 75 KW are however not attainable at present, but the con- 
cept seems feasible and is probably the only way to power a rover. 

The concept of laser beams powering a static ground station to which a 
rover would return to recharge its batteries is not viable as the solar 
collector dishes could be placed directly on the ground station. As the 
ground station could only collect solar energy at night, energy storage at 
both station and rover would be required. 

2.3 Communication and Surveillance 

Solar powered lasers could be used for line - of - sight communication and 
surveillance either in a space environment or left unattended in remote 
terrestrial locations. The requirements of high power are unnecessary 
here, hence lasers of low power such as have been investigated at Langley 
could be used. For operation on a continuous basis, even if the solar col- 
lector is in shadow part of the time, it is possible to use "black body 
pumped" lasers where a cavity is heated, and the heat energy stored to con- 
tinue pumping the laser while the satellite is in the planet's shadow. For 
example, a geosynchronous satellite at 3*6 x 10 7 m above the earth's sur- 
face is in the earth's shadow for about 1 hour and 10 minutes every 24 
hours. Terrestrial solar pumped lasers could only operate during sunlight, 
unless they were black body lasers capable of storing heat over a period of 
over 10 hours . 
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3. GENERAL REQUIREMENTS FOR SOLAR PUMPED LASERS 
AS SOLAR ENERGY CONVERTERS 

We consider the general requirements for solar pumped gas lasers 
acting as energy converters before we deal with specific systems. 

The laser must be efficient in converting the solar energy into an 
output beam. The efficiency can be subdivided into a product of effi- 
ciencies according to the sequence of physical processes occurring. 

The criteria for a solar pumped laser are as follows: 

1. There must be broadband absorption; 

2. Peak absorption should occur near the peak of the solar 
spectrum; 

3. The absorbing medium (gas) must be at sufficient pressure to 
absorb most of the photons , otherwise a multipass arrangement 
must be used. 

4. High quantum yield into a long-lived (metastable) state which 
serves as the upper laser level is necessary. 

5. In general, quenching of the excited state should be small, but 
this condition is alleviated if 

6. The lower level is rapidly depopulated to maintain inversion; 

7 . The upper and lower levels must be sufficiently separated to 
yield a reasonable quantum efficiency; 

8. The process must be reversible; if not the components must be 
reconstituted by flow methods. 

9. Operation at high temperature would be advantageous to reduce the 
weight of the radiator 

10. The beam must diverge as little as possible. 

Items 1 and 2 are included in a "solar utilization efficiency", item 3 is 
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included in an "absorption efficiency", items 4, 5 and 6 are included in a 
"kinetic efficiency", while item 7 is included in the "quantum efficiency". 
All of these terms are defined more fully in Appendix A, a paper in the 
AIAA J of Propulsion and Power. 

In addition to efficiency another criterion for space applications is 
the output power per unit weight, as the lasers have to be lifted off a 
launch pad. Items 8 and 9 are of concern here. This also is considered in 
Appendix A where the effects of high temperature operation are outlined. 

When the beam eventually leaves the laser, the far field pattern is of 
importance, and the divergence of the beam must be small. This aspect is 
item 10 above. The fraction of the output power that strikes a given 
target is termed the "collection efficiency". A study of far field pat- 
terns is given in Appendix B, Technical Report PTR-89-6. This study is 
different from any previous studies in that for large lasers it examines 
the effect of non-homogeneity of the medium, as would occur in a gas laser 
pumped from the outside. This study was not completed when the grant was 
terminated, but the preliminary results imply that lasers of large diameter 
and non-uniform media may suffer considerable divergence of their beams. 

The general criteria explained in Section 3 will next be applied in Section 
4. 

4. SPECIFIC LASER SYSTEMS 

This study has concentrated almost entirely on gas lasers because for 
space applications, there are no stringent requirements on size (in 
contrast to weight) and also the medium would be uniform. In addition to 
gas lasers however, a limited study was made of dye lasers in liquid or 
vapor form which is included in section 4.3. 
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For gas lasers, investigations concentrated both on different lasing 
materials and also on different laser systems. The families of lasers 
investigated are shown in Figure 2. They are subdivided according to how 
the radiation is absorbed and utilized and then further subdivided 
according to specific mechanisms. The first subdivision is whether part of 
the solar spectrum is absorbed, (Section 4.1), or whether the whole solar 
spectrum is utilized (Section 4.2). Utilizing the whole spectrum may not 
be advantageous as we shall see. 

4.1 Solar Lasers Which Utilize Only Part of the Solar Spectrum 

Lasers of this type can be further subdivided according to whether the 
lasing medium is dissociated or not. If dissociation occurs then one of 
the species must be in an excited state, which forms the upper lasing 
level, and this is dealt with in Section 4.1.1. The non-dissociative types 
can be metal vapors such as Na2 , Cs2 etc, where absorption of a photon 
results in the molecule being transformed into an electronically excited 
state-Section 4.1.2. 

4.1.1 Dissociative Laser Media 

We further subdivide dissociative lasers according to whether they 
possess a separate absorbing medium which forms an excited state, and that 
medium in turn then collides with and hands over energy to a second medium 
which acts as the lasant. The processes of absorption and lasing are thus 
separated, and each tailored to give the highest laser efficiency. Studies 
of this type are summarized in section 4.1. 1.1 below. Following that 
section is a study of solar lasers that use the same material to both 
absorb the solar radiation and lase. 
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GAS-FXLLED 
SOLAR PUMPED 
LASERS 



X- 3 halogen, Y 3 different halogen atom 
R 3 complex radical 


Figure 2. Various families of solar pumped lasers. 
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4.1. 1.1 Dissociative Lasers Which Use a Separate Absorber to Hand Over 
Energy to the Lasant 

For this type of laser an investigation was made both of the physical 
processes involved and especially of what would be the best materials for 
lasing. The results are summarized in Appendix C, the paper "Solar-Pumped 
Electronic- to -Vibrational Energy Transfer Lasers" by W.L. Harries and J.W. 
Wilson, Space Solar Power Review 2 . 367 (1981). 

4. 1.1. 2 Solar Lasers Which Use One Material To Absorb and Lase 

The results of the investigation in Section 4. 1.1.1 implied that for 
all the materials considered, where a different material absorbed the 
radiation and another lased, that the efficiencies were below 5 x 10 . It 

was thought that the efficiency might be improved if the hand-over process 
was eliminated and a medium chosen which performed both the process of 
absorption and lasing. Reasons for the choice of material and the laser 
mechanisms were investigated and the results are summarized in Appendix D, 
the paper "Kinetic Modeling of an IBr Solar Pumped Laser' by W.L. Harries 
and W.E. Meador, Space Solar Power Review 4, 189 (1983). 

4. 1.1. 3 Summary of Dissociative Lasers 

The above studies indicate that as far as the materials investigated 
that for lasers with a separate absorbing medium the efficiencies were 5 x 
10 * 3 or less. Where the same material (IBr) performed both absorption and 
lasing, the efficiencies could approach about 0 . 012 . 

Referring to Figure 2, the types X 2 , and XY where X is a halogen and Y 
a different halogen atom have been covered above. Not covered is the type 
RX where R is a complex radial e.g. R - C 3 F 7 or C 4 Fq etc. and X is an 
Iodine atom. Considerable theoretical and experimental work was performed 
at Langley on this type of solar pumped laser, and the lasers operated 
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successfully driven by a solar simulator. Their absorption bands are 
included in Fig. 6 in Appendix C, but both theory and experiment at Langley 
were done by others and they are not included in this report. 

Lasers of this type suffer the disadvantage that during lasing, the 
compounds break up into constituents and a flow system has to be introduced 
where the constituents have to be recombined to form the lasant. Also 
their absorption band occurs in a narrow band in the ultra violet (see Fig. 
6 in Appendix C) and their solar efficiencies are low. 

4.1.2 Non Dissociative Laser Media 

Hitherto, the types of molecular lasers considered have been those 
where photodissociation occurred, but absorption of a photon can also yield 
an excited molecule in an upper electronic state (Fig. 3). The transition 
occurs according to the Franck-Condon principle, and the molecule could 
return to the lower electronic level by emission of a stimulated photon. 

A comprehensive article by Wellengehausen (IEEE J of Quantum 
Electronics, QE-15, 1108 (1979)) reported on the development of optically 
pumped cw dimer lasers^) and listed a number of experiments where lasing 
had occurred in Na 2 , Te 2 , Br 2 , Li-2> K 2 , and S 2 , as well as in the halogens 
12 and Br2» The pumping in the experiments was done by lasers, either 
frequency doubled Neodymium (532.5 nm) or an argon ion laser (472.7 nm) . 
These frequencies lie near the peak of the solar spectrum. The effi- 
ciencies quoted were up to 15%, but this has to be further reduced by a 
"solar efficiency" for solar pumping. 

If the dimer were solar pumped (Fig. 4), the Franck-Condon transitions 
would result in vibrational levels much higher than the lasing level, shown 
here near the upper level minimum. If it were possible for the molecule to 
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Figure 3 . Laser cycles between bound electronic molecular states, (a) 
three-level cycle and (b) four-level cycle with radiative or 
nonradiative transitions between the levels 4 and 3. 
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Figure 4. Energy level diagram for a solar pumped dimer molecule M2. The 
vibrational levels only are shown. A molecule raised to the nth 
level of the upper electronic A state can either descend to the 
(n - l), (n*2) ...» levels by collisions, (continuous arrows) and 
radiation (dotted arrows) or else radiate a photon (wave arrow) 
and then drop to the lower electronic X state. 



cascade down in energy without falling to ground level, then there would be 
a " funneling" to the upper laser level. The fraction that arrives at the 
upper laser level depends on the probability at each level of spontaneous 
emission with relaxation to the lower electronic level, and on the 
probability of descending to a lower vibrational-rotational level in the 
upper electronic state by collisions or radiation. For homonuclear 
diatomic molecules, radiative transitions between rotational vibrations of 
the same electronic state are not allowed. 

If the "funneling" is efficient, then this type of laser would be very 
attractive. However, the wider the absorption band, the higher will be the 
initial vibrational level in the upper electronic state, and the less the 
probability of arriving at the upper lasing level without radiating. Also, 
the higher the initial vibrational level, the less will be the effective 
quantum efficiency. Hence, compromises will have to be made between solar 
kinetic and quantum efficiencies. 

The possibility of lasing by solar pumping in metallic vapors such as 
those mentioned seemed hopeful from the point of view of efficiency. Also, 
the metal vapors would be at high temperatures, and high temperature lasers 
would need a smaller radiator, which would decrease the weight and hence 
improve the power to weight ratio. Lastly and very importantly was the 
possibility that these lasers might emit in the visible. 

Considerable study was made of metal vapor lasers, but unfortunately 
as was pointed out by W.E. Meador, although experiments where lasers were 
used to pump Na2 seemed to imply reasonable thresholds , the threshold 
levels for solar pumping were much too high for practical use. Normally 
for solar pumped lasers, the solar collector concentrates the radiation. 

The upper limit of concentration is believed to be around 20,000-- a 
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limitation of the optics based on the fact that the sun is not a point but 
distributed source. Accordingly, the investigation of lasing by metal 
vapors was deemphasized around 1987, but as several aspects of metal vapors 
as lasing media had been completed and as they might be useful for future 
studies of optical pumping, they are included below. 

The problem of self-absorption of stimulated emission is a general one 
and the case of sodium is treated in Section 4. 1.2.1. Another important 
question is the ratio of dimers to monomers in metallic vapors which 
depends on temperature. Although curves existed in the literature, we were 
unable to find a reasonable mechanism or physical picture and a simple 
treatment and formulae are given in Section 4. 1.2. 2. The problem of 
metastability of a dimer molecule is also a general one in lasers and is 
treated in Section 4.1.2. 3 

At the same time as these studies were performed, an experimental 
program under Dr. N. Jalufka of NASA/Langley was undertaken to study 
absorption and emission in Na 2 vapor. A new electronic transfer process 
was suggested by Jalufka, and the experiments and theory were written up as 
a paper entitled "Electronic Energy Transfer in Molecular Sodium" by N.W. 
Jalufka and W.L. Harries and submitted to the Journal of Molecular 
Spectroscopy. This is included in Appendix E. Calculations on a possible 
sodium laser are shown in Section 4. 1.2. 4, A and on Lithium in Section 
4. 1.2.4, B. Although it is felt that these metallic vapors may not be 
useful for solar pumped lasing, these calculations might be useful in other 
laser applications. 
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4. 1.2.1 Stimulated Emission and Self Absorption in Sodium Vapor 

A. Introduction 

In any lasing medium, the emission wavelength should be chosen where 
there is little self absorption. As emission and absorption spectra for 
metallic vapors did not seem available, therefore, estimates were made of 
these cross sections for sodium vapor as functions of wavelength. Although 
absolute values were not obtained, information on where the emission wave- 
length should occur became evident. 

In Section 4. 1.2. IB, the method of obtaining quantities proportional 
to the cross sections versus wavelength is outlined. A further comparison 
based on alternative expressions for the absorption and emission cross 
sections over a limited wavelength range is made in Section 4. 1.2. 1C which 
supports the evidence of Section 4. 1.2. IB. 

B. Absorption and Emission Probabilities for Na2 Vapor 

The probabilities can be estimated if the Franck-Condon factors F(v , 
v * ) are known for transitions from the various vibrational levels v in the 
ground state X 1 2* of Na 2 , to vibrational levels v' in an electronically 
excited state. An excited state which would correspond to absorption in 
the visible is the B 1 n u and fortunately, a very complete table of Franck- 
Condon factors was given by Kusch and Hessel.( 3 > The amount of absorption 
from any level v" would depend on the relative population of the level to 
the other levels in the same electronic state. The levels would be filled 
according to Boltzmann statistics. Similarly, the emission would depend on 
the relative population of the upper vibrational levels . 

For the ground electronic state the energy of the level v" in cm-1 was 

E(v" ) - W (v n + - ) - W x ( v " + -) 2 (4 ' 1) 

e' „ e e 2 
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where - 159.1 cm"^, and - 0.7254 cm"^ and the lowest energy - 

E(v" - 0) . ( 4 > 

The probability of filling level v M was 

E(v" ) "E q 

P(v") - exp (4-2) 

kT 

when k is Boltzmann's constant and T the temperature in degrees Kelvin. 
The fraction of molecules in this state was 

f(v-) - n P(vn) (4-3) 

E P(v") 
v" - 0 

where n was taken as 45 . 

The minimum of the potential curve for the B^II U state was at 20,319 
cm’l, and the energy of the v' state was 

E(v') - w ' (v' + 1 ) - w x ' (v* + - ) 2 + 20,319 (4-4) 

e 2 e e 2 

where w e ' - 124.9 cm’^, and w e x e ' - 0.6999^), and the lowest vibrational 
energy level E q* - E(v' - 0). The fraction of molecules filling level v' , 
namely f(v') was calculated as in equations (2) and (3), except n was now 
taken as 28. The number n for both ground and excited electronic states 
was determined by the availability of Franck-Condon factors. Both numbers 
are sufficiently large so that inclusion of any higher levels would not 
change the results appreciably. 

A quantity proportional to the absorption cross section, which is a 
function of wavelength A is^) 

Q(A) - F(v" , v') f(v")i/ (4-5) 

where F(v" ,v') is the Franck-Condon factor for the transition v"-*v' , u is 
the frequency of the absorbing photon determined by E(v') - E(v"). Plots 
of Q(A) vs. A for transitions from level v" to v' corresponding to 
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0 < v 


" < 45 , or 0 < v' < 28 are shown in Figure 5. The overall outline 

resembles a Gaussian absorption curve. 

The values of Q(A) were then averaged over an interval of 5 ran to give 
smoother curves and are shown in Figure 6 for temperatures of 500, 650 and 
1000° C. The curves now resemble Gaussian absorption curves and increasing 
the temperature decreases the peak value and widens the wings , as it 
should. 

An estimation of a quantity proportional to a stimulated emission 
cross section can be made in a similar manner. The probability of filling 
a vibrational level in the upper electronic state was again assumed to 
occur according to Boltzmann statistics. This assumption, however, may 
lead to error because the levels in the upper state are filled from the 
ground state, and their population densities are dependent on the rate of 
arrival from the ground state. If they remain in the upper state for many 
collisions, a Boltzmann distribution should result, but if they are quickly 
quenched such a distribution is unlikely. 

Assuming the fraction in level v' is f(v') as previously calculated, 
then a quantity R( A) proportional to the stimulated emission cross section 
is given by^^ 

R( A) - F(v" , v' ) f (v' ) v 2 (4 ' 6) 

Plots of Q(A) , R(A) vs. A are shown in Figure 7. The values have 
been smoothed from those of Figure 5 by averaging over a 5 nm bandwidth. 

It should be noted that for both quantities, the ordinate is arbitrary, and 
that they are only proportional to the absorption and stimulated emission 
cross sections. Hence, if Figure 7 is regarded as a plot of the cross 
sections then they are not normalized with respect to each other. The 
temperature values are 500° C and 1000° C. 
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Plots of Q(x) vs \ for different temperatures where the values have been averaqed over 
5 nm intervals. 







It is evident that the R(A) curve is shifted to the right of the Q(A) 
curve, suggesting a region between 500 and 540 run where stimulated emission 
would dominate over self absorption. 


C. Alternative Method of Comparing Emission and Absorption Cross 
Sections 

The absorption cross section for a molecule is given by(^) 


where A ^3 is the Einstein 
absorption wavelength, and 
section is given by 


13 .2 

a - A 

a 3 a 

8*c 

coefficient for the transaction 
c the velocity of light. The 


1 


(4-7) 

-►3, A a is the 
emission cross 


a 

e 


A 4 

e 32 


4w 2 c AA 


(4-8) 


e 

where A e is the emitted wavelength, A 32 the Einstein coefficient for the 

transition from level 3 , the upper lasing level to level 2 the lower lasing 

level AA e here is the emission bandwidth. If we consider the possibility 

o 

of absorbing the wavelength A , then A - A and the ratio — is then 

e a e 


a 


o 

e 


2c 2 A 32 


A 2 
e 


a n A - 0 AA 
a 13 e 


(4-9) 


and increases with A 2 . The expression may hold in the region between 5 and 

e 

5.4 nm in Figure 7, but clearly does not hold elsewhere. 

Assuming A e - 525 nm, A 13 - A 32 - 1.9 x 10 7 s-1 and that AA e is deter- 
mined by Doppler broadening, (Ai/ - 3 Ghz, AA e — 9.2 x 10 13m)( 7 ) then 

a 

a - 6.9 x 10 ' 14 cm 2 , a - 1.3 x 10 ' 12 cm 2 and the ratio — - 19 . 
a e 0 

a 
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D. Conclusions 


The above independent check of the values of a e and a a in the 

a 

wavelength range 5 to 5.4 run showed — - 19. The plots of Q(A) and R(A) 

a 

a 

R( A ) 

are not normalized, but at 5.25 nm show the ratio L to be about 3, 

Q(A) 

suggesting the R(A) curve should be about six times higher. The form of 
the curves show that the only region where lasing should be attempted in 
Na2 is between 5.2 and 5.4 nm. This conclusion is consistent with the fact 

that lasing at 5.25 nm has been reported by Wellegehausen. 

4. 1.2. 2 Densities of Dimers and Monomers in Metallic Vapors 

A. Summary 

The dimer/monomer ratio is derived for metallic vapors as a function 
of temperature. The basic physical mechanisms are dissociation caused by 
collisions with energetic atoms, and three body recombination in the gas 
phase, and the creation of a saturated vapor at the gas -liquid interface. 
Increasing the temperature increases the collision and recombination rates 
so the dimer/monomer ratio rises. The effect of Na 2 acting as a third body 

in the recombination process is included for the first time. 

B. Introduction 

The possibility of using metallic vapors in lasers for solar energy 
conversion is appealing as they could operate at high temperatures. The 
lasing medium would be the dimers, which in many instances consist of only 
a small fraction of the gas particles. Plots of this fraction as a func- 
tion of temperature were given by Lapp and Harris, from results of Stull 
and Sinke(^) but without explanation. The purpose of this paper is to 
understand the processes which are occurring and to duplicate these 
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results. The effect of the dimers acting as third bodies are included for 
the first time. 

C. Physical Processes Occurring in Metallic Vapors 

The vapor consists of a mixture of monomers and dimers boiled from the 
liquid phase, with a metallic wick connecting the regions of high and low 
temperatures. In the vapor phase, it is assumed that the reactions of 
dissociation and recombination can occur as shown in Table I. 

Table I. Reactions Occurring in the Vapor Phase 


Type of reaction 


Dissociation 

a 



Na 2 + Na' -*• Na + Na + Na 

Three body recombination 

b 

k 

Na 4 * Na + Na 5 Na^ + Na 

— i 

c 

k 

Na + Na + Na 2 * Na 2 + Na 2 


It is assumed that dissociation of Na2 occurs according to reaction a 
in Table I where kf is the rate coefficient, and Na' is an atom with 
kinetic energy exceeding D 0 , the dissociation energy of Na 2 , and hence in 
the tail of a Maxwellian distribution: 

Na'/ Na “ ex P (-D /kT) (4.10) 

Recombination in the gas occurs by three body collisions, and it will 
be assumed first that the predominating third body is Na, an assumption 
that works well at low temperatures as will be seen later. Hence, reaction 
b of Table I dominates. The rate equation for Na2 is then 

d(Na 0 ) 3 

- -kf Na£ Na' + k^ (Na) (^.H) 

dt r 


25 



In the steady state, equations (4-10) and (4-11) yield 


- Jl exp (D /kT) - K (4.12) 

2 k ° 

(Na) K f 

where K is a constant for a given temperature. 

Fortunately, reference (8) gave values of the partial pressures of 
dimers and monomers at a fixed temperature for several metals and if D 0 is 
known, then the quantity k r /kf can be calculated. Equation (4-12) can then 
be used to calculate Na 2 /(Na)^ for any temperature. 

The saturated vapor pressures of metals are well known and usually 
given in the form(10) 


L °g 10 P 


0.05223 u 
a + b 

T 


(4.13) 


where p is in torr, T in degrees K and a and b are characteristic 
of the metal. The expression is valid from 180°C to 883°C. 

The vapor pressure is also 


p - (Na 2 + Na) kT (4.14) 

and on rearranging equations (4-13) and (4-14) 

Na + Na - — C exp(-Q/kT) (4.15) 

2 kT 

where C and Q are obtained from equation (4-13). Equations (4-12) and 
(4-15) then give solutions to Na, Na2 , and Na 2 /Na. 

Data for four metals are given in Table II, and plots of Na 2 /Na are 
shown in Figure 8. The plots showed excellent agreement with the curves of 
Stull and Sinke.(^) For all the vapors, the dimer-monomer ratio increases 
with temperature. The reason is that at higher temperatures the densities 
of the particles increase and hence, the recombination rate which is depen- 
dent on three body collisions increases faster than the dissociative rate. 
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Table II. Data for Metals 


Vapor pressure 

constants 

Partial Pressures at 700K (torr)(l>2) 

0 o (eV) 


a 

b 

Monomer 

Dimer 


Na 2 

103300 

7.553 

0.69 

2.2 x lO- 2 

0.730 

Rb 

76000 

6.976 

19.40 

0.454 

0.490 

Cs 

73400 

6.949 

27.60 

0.625 

0.450 

K 

84900 

7.183 

7.24 

8.2 x lO’ 2 

0.510 


Three body collision where Na2 is the third body have been neglected above. 
Also dissociative collisions from Na2 molecules in the tail of the 
distribution. The error introduced is small within the temperature range 
of Figure 8, as at most they constitute less than 8% of the total 
particles . 


D. Effect of Na 2 as a Third Body 

The Na2 molecules acting as third bodies contribute an additional 
recombination reaction (c) to reaction (b) in Table I. A rate coefficient 
k s is introduced, different from k r , but unfortunately, the values of k r , 
k s are not known. The rate equation for Na 2 now becomes 
d(Na ) 

_ - - k~Na 0 Na' + k Na + k Na Na„ (4-16) 

dt r 2 r s 2 


and in the steady state, using equation (4-10), we obtain 


Na„ 


2 k s 

Na + — NaNa„ 


— exp (D q /KT) 
k f 


(4-17) 


Equations (4-17) and (4-15) then given solutions of Na, Na2 and Na 2 /Na. 
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Plots of Na 2 /Na vs. temperature are shown in Figure 9 for the case of 
sodium where the effect of the dimers as third bodies is shown in terms of 
a parameter a defined as a - k s /k r . Thus, when a — 0 we revert to the 
case of Na2 in Figure 8. The temperature scale is shown to be 1000 °C 
although the range of validity of Equation (4-13) which leads to Equation 
(4-15) is only up to 883°C. As the ratio k s /k r is not known, the behavior 
of k s can only be illustrated by the parameters a - k s /k r . 

E. Conclusions 

The dimer/monomer ratio calculated as a function of temperature for 
four different metal vapors showed excellent agreement with the curves of 
Stull and Sinke,( 9 ) indicating that dissociation, three body recombination, 
and the creation of a saturated vapor were the basic processes occurring. 
Inclusion of the dimers at third bodies, in recombination collisions, 
increased the dimer/monomer ratio slightly at higher temperatures. 

4. 1.2. 3 Metastability of N 2 , Na2 , and CO 

N 2 and CO when vibrationally excited are long lived and metastable, 
whereas the vibrational levels in Na2 are short lived. An investigation 
showed that polarizability was not the cause of the short lifetimes- -the 
Na2 had the smallest polarizability. A comparison of the vibrational 
frequencies, w e , showed that the metastable N 2 and CO had much greater 
values by a factor of over 10, to the Na 2 ; the reduced masses, /i, were of 
the same order for all three. The spring constants of the molecules, 
proportioned to were about 100 times greater for N 2 and CO than for 

Na2 . 
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Fig. '9. Effect of Na? as a third body on the dlmer/monomer ratio vs temperature. 
The quantity a Is equal to k s /k r where k. Is the rate coefficient for 
three body recombination where Na 2 Is the third body, k r where Na Is 
the third body. The case a a 0 corresponds to the curve for sodium in 
figure 2.1. The curves were valid only up to 883°C. 
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The spring constant determines the shape of the potential energy 
curves for the molecules; these are very flat with closely spaced vibra- 
tional energy levels for Na 2 . In contrast they are very steep sided with 
widely spaced vibrational energy levels for N 2 and CO. In the case of Na 2 , 
the vibrational levels have a spacing of the same order as the thermal 
energy of the molecules, -0.02 eV and thermal collisions can deactivate the 
Na 2 . With N 2 and CO, the spacing is of order 0.3 eV--much greater and 

consistent with the long lived behavior. 

The short time of transfer from one vibrational level to the other in 
Na 2 may be advantageous as the energy from levels which are several quantum 
numbers above the upper laser level could be transferred to it. The popu- 
lation of the lower laser level would also in turn be quickly depleted by 
thermal collisions. 

4. 1.2. 4 Specific Metallic Vapor Laser Systems 
A. The Na 2 Metallic Vapor Laser 
A. 1 Photon Absorption 

The Einstein coefficient for transitions from the lower level 1 to the 
upper level 3 (Figure 10) of the Na 2 laser is A^ - 107 s-l.< U) The 
absorption cross section is large, a a - A^ 3 / 8 * , ci /^^3 - 3x10 cm for 
absorption around 500 nm. The absorption length L then turns out to be 
10' 3 /p (cm) where p is the pressure of Na 2 vapor in torr. For reasonable 
active lengths, low pressures are necessary. The estimate of L is 
consistent with lasing in supersonic beams of Na 2 with an active length of 
order 1 mm. However, the fraction of the vapor which is in the form of Na 2 
is about 0.04 by number, with the remainder in the form of Na at around 
700K where the laser operated. ( 12 ) 
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Figure 10. Energy level diagram for a solar pumped dimer molecule 
* The vibrational levels only are shown. A molecule 
raised to level v 1 of the upper electronic A state can 
descend to the (v'-l), (v»-2)..., levels by collisions 
(continuous arrows) (and radiation-dotted arrows for 
heteronuc lear molecules) and remain in the upper elec- 
tronic state. Otherwise, it may radiate a photon 
(wave arrow) or be quenched and then drop to the lower 
electronic X state. 
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A. 2 Choice of the Lasing Levels 

The lower laser level 2 (Figure 10) must have a population of less 
than half that of the ground level 1 so that it can be adequately 
depleted. The reason is that the maximum population density achievable 
in level 3 with a strong pump source is n^ — . This means that 

(neglecting statistical weights) the density of level 1 at temperature T 
is divided into 2 between levels 3 and 1: n 3 ~ ” n^(T)/2 where n^(T) is 

the normal Boltzmann population. Population inversion is therefore 
possible only if n (T)/2 > n 2 (T).( n > At 800 K, the lowest value of the 
quantum number v" for level 2 is v" — 3. Lasing to the v M — 3 level would 
originate near the v' - 0 level in the upper electronic state. The low 
value of v * is helpful in utilizing most of the populations of states with 
higher v' , provided they can relax downward. Level 2 must be depleted by 
collisions thus requiring a buffer gas such as Argon. 

The closeness of the vibrational energy levels in both ground and 
electronically excited states raises the question of whether there could be 
lasing from more than one upper level to a corresponding lower level. The 
two transitions would have to be equal in frequency within the emission 
bandwidth of 3 GHz, but this is unlikely because the spacing of the vibra- 
tional levels is around 6x10^ Hz or 2000 times greater . 

A. 3 Requirement for Collisional Relaxation in the Upper Electronic State 

The collisional relaxation rate for levels near the upper laser level 
must exceed the spontaneous emission rate (neglecting quenching for the 
moment) , other wise the density of molecules in the upper level will not 

increase before threshold. If Argon collisions dominate, the requirement 
is: 

(Na 2 ) (Ar)a 3 v 3 > (Na 2 ) 7 ° 3 (4-18) 
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where the quantities in parenthesis are densities, a 3 the cross section, V3 
the relative velocity and 7° the spontaneous emission rate. With a 3 - 

5 x 10 ^ cm^, v^ - 7 x 10^ cm ^ and 7® - 1.4 x 10** s^ (Appendix D) 
the inequality is satisfied for Argon pressures greater than 11 torr. 

A. 4 Requirement for Stimulated Emission to Exceed Collisional 
Relaxation 

Lasing requires the stimulated emission rate to exceed the collisional 

relaxation rate out of v^ of the upper state. The stimulated emission 

depends on the inverted population density and emission cross sections. 

Assuming the upper level to be B^II U , with an angular momentum quantum 

number of 1, and the lower level to be X^S + with no angular momentum, then 

S 

the multiplicities for the upper and lower levels are g^ - 2 and g^ - 1 . 

The inverted populated density is: 

An - <Na r V) - 2(Na 2 ,v”) (4-19) 

where the quantities in parentheses are the densities in the upper and 

lower laser levels, respectively. It follows that (Na 2 , v^) must exceed 2 

4 2 2 

(Na„ , v 1 !) for lasing. The emission cross section a -A A on /47r cAA cm is 
2 i 0 e e 32 e 

evaluated at 1.3 x 10-12 cm2.(H) Here A e is the emission wavelength 525 
nm, A32 - 1 . 9 x lO? s “ ^ and the emission bandwidth Ai/ — 10^ Hz 
(corresponding to Doppler broadening at 800K) . The above value for o e 
agrees roughly with the quoted value 3 x 10"^^ cm^. The requirement 

that stimulated emission exceeds depletion by collisions of the upper level 
is : 

nc a AN > (Na 0 , v') (Ar) v Q (4-20) 

e l u J j 

where n is the photon density during lasing. For pressures corresponding 
to 0.04 torr of Na2 and 0.96 torr of Na and 60 torr of Argon, this is 
equivalent to requiring nAn/(Na2, v^) > 2 x 10^^(cm‘^) which should be 
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easily satisfied if nAn > 10 12 and less than 10 percent of the Na2 atoms 
were in the upper laser level. 

A. 5 Requirement for Collisional Relaxation in the Lower Laser Level 

For CW operation the total relaxation rate for the lower laser level 
must be greater than the spontaneous transition rate for a specific 
transition 3-2 or (Na^ v«) (Ar) > (Na 2> V) A^. With the above 

numbers and A 32 - 1.9 x 10 7 sj i ' 1 , the condition is equal to (Na 2> v^) / 

(Na^, v£) < 40, which is easily realized, 

A. 6 Probability of Filling the Upper Laser Level 

The energy level diagram of Na 2 (Figure 11)< 13 > shows numerous vibra- 
tional levels (of quantum numbers v' and v" in the upper and lower elec- 
tronic states) spaced closely together. At 800 K, there is a one percent 
probability of filling the v" - 10 vibrational state, and an absorption 
transition then populates an upper vibrational level with v' about 30. 

The energy must be transferred to the upper laser level without being 
dissipated from the molecule. A rough estimate of the probability of 
arriving at the upper lasing level v^ from a higher vibrational level v' 

will now be made. 

The molecule can drop from a higher to a lower vibrational level by 
collisions. A model by Landau and TellerO^) enables an estimate to be 
made of the number of collisions required for a transfer of vibrational 
into translational energy and the number is of order 1 for Na 2 ( 

The vibrational levels in the upper electronic state are about 120 cm 1 
apart corresponding to 170 K. If the gas temperature is 800 K, assuming 
the energy exchanged is equal to half the kinetic energy, then the jump in 
vibrational levels corresponds to AV - ± 2, and the transition can be 
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Figure 11. Energy diagram for showing the lower level <X) 1 E+ 
and the upper level There are 109 vibrationaf 

levels in the lower electronic state; only 45 are 
shown. There are 88 levels in the upper state; 50 are 
shown. The vertical dotted lines correspond to Franck 
Condon transitions upward from the v M * 0 level. 

The solid lines correspond to downward transitions 
from the v' - 0 level in the B*It state which end 
in the v" - 3 level. u 






either upwards or downwards. Such collisions can be due either to Na 2 
molecules or Argon atoms (<73 - 5 x lO'^- 5 ).^^^ Collisions of this kind 
leave the molecule still in a vibrational level in the upper electronic 
state . 

The rate equation for relaxing from level v* to level (v* ± Av' ) is: 

L_ (Na v') - (Na 2 ,v') (Na 2 ) 0 ^ + (Na 2 > v') (Ar)*^ (4-21) 

dt 

r 

where (Na 2 ,v') is the density of Na 2 molecules of vibrational quantum level 
v', and a - 10' 14 cm 2 and a 3 - 5 x 10' 15 cm 2 are the cross sections for 
collisions with Na 2 molecules of relative speed v X and Argon of relative 
speed V 3 . All collisions with Na are assumed to result in quenching, < 11 ) 
(see below) . 

The processes which cause the molecule to relax from the upper to the 
lower electronic state (apart from stimulated emission) are spontaneous 
emission of a photon (rate coefficient 7 ^ “ 1-^ x 10 s )(H) and 
quenching. Hence, the rate of spontaneous emission from the upper 
electronic state to the ground state is: 

jL (Na 2 ,v') - (Na 2 ,v') 7 3 (4-22) 

dt 

s 

Quenching can be due to Na 2 or Na. The quenching cross section for 
Na 2 , 04 is probably around 10 "^ cm^ , but the fraction of Na 2 is small. 

The quenching cross section of Argon is assumed to be negligible as it has 
no degrees of freedom. The cross section for Na is large ( 0 5 — 6 x 10 
cm2) (11) xhe reason is that the Na forms complexes with Na 2 which have 
internal degrees of freedom and can absorb the energy drop of 2 x 10 4 cm- 1 , 
before splitting up into an Na 2 molecule in the lower electronic level 
possessing vibrational and rotational energy and an Na atom, both with high 
kinetic energy. 
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(4-23) 


(Na^ , v' ) + Na ■+ Na^* -► Na 2 + Na + AE 
The quenching rate from vibrational quantum level v' is: 



(Na 2 , v' ) - (Na 2 ,v') (Na^cr^ + (Na 2 ,v' ) (Na)^ 


and V 2 is the relative speed of Na. 

The probability that a collision will leave the molecule in the upper 
electronic state is then: 


f - — (Na v') • I 

dt z dt 

r ( r 

The Na 2 laser experiments were performed at a total pressure of Na + 
Na 2 of about one torr with Argon pressures of many tens of torr. A plot of 
f vs. Argon pressure for these conditions for different percentages by 
number of Na 2 (from 2 to 20 percent) (Figure 12) shows that f increases 
with Argon pressure but is not very sensitive to the Na 2 fraction (and 
therefore also to o\ and < 74 ) . 

The experiments with laser pumped Na 2 vapor showed that the power 
output dropped to zero at about 70 torr of Argon, and the reason quoted was 
the increase of the relaxation rate of the upper laser level due to Argon 
collisions . Ul) it is possible this effect would not occur in a solar 
pumped laser, because the relaxation process would supply the upper laser 
level from levels above v' , as well as deplete it. Should pressures above 
70 torr be unworkable then f < 0.7 (Figure 11). On the other hand, if it 
were possible to use 300 torr of Argon then f » 0.9. The quenching by 

Argon is assumed small. (12) 

There are 88 vibrational levels in the B^-IIu state, and levels up to 
about n — 30 will be occupied. The probability of arriving at v' — 0 from 
a higher v' is next roughly estimated. The spacing of the vibrational 
levels is assumed constant at w e cm'V The collisions that transfer 


(Na v') + -JL (Na v') + (Na v') 

dt Z dt. I 

s Q I 


(4-24) 
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Figure 12. Plot of f, the fraction of molecules that stay in the 

upper level per collision vs Argon pressure for different 

fractions of Naj in a total pressure Na + Na2 of one 
torr. 
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vibrational to translational energy can raise or lower the molecule 

from vibrational level v' to v' ± Av. The molecule as it changes its 

vibrational levels through collisions, performs a one dimensional random 

walk in energy space with a step length Av'u e (cm - ^) and has to travel a 

distance (v' - v' )w to arrive at the upper lasing level v' . The 
u e e 17 ° u 

"diffusion coefficient is (Av' u e ) 2 y c where u c is the collision frequency. 

2 2 2 

The diffusion equation (v' - v') w - (An w ) v t, gives the number of 

u e e c 

collisions N - u c t as: 

N “ ( v ' - v;) 2 / (Av') 2 (4-25) 

The equation is approximate as it assumes a constant step length, but is 
sufficient for illustration. 

If a fraction f of the molecules stays in the upper electronic state 
at each collision, then in N collisions, the fraction remaining is on the 
average : 

F - r* (4-26) 

e.g., if for Na 2 f - 0.7, Av' - 2, v' - v' - 10, then F - 1.3 x 10'^. But 
if f - 0.9, then F - 0.39. 

Assuming all the vibrational levels above v^ are filled with equal 

probability (not too inaccurate an assumption because only levels near v' 

contribute) , then the contribution to the level v' from all levels is 

u 

GO 

s - l f [(v' - v ' ) /Dv ' ] 2 (4-27) 

v'-v'+l 

u 

The summation converges rapidly with increasing v' especially if f is 
small. For example, if f - 0.7, Av' * 2 then S - 2.47, but if f - 0.9 then 
S - 4.96. This is equivalent to saying that the first 2 (for f - 0.7) and 
5 (for f - 0.9) levels, respectively, are the only ones that contribute to 
the upper lasing level. 
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If levels above n — 5 are essentially unusable for lasing, then Figure 
11 shows the absorption bandwidth would correspond to 15 ran, the region 
between 476.6 to 490.6 nm. The solar efficiency v s for this narrow region 
is less than two percent. 

A. 7 Overall Efficiency of Na2 Vapor Laser 

The absorption efficiency r; A can be made to approach unity with very 
low Na 2 pressure if a path length of 1 cm is used. A rough estimate of the 
fraction of molecules of vibrational levels v' < 5 that arrive at the V - 
0 level can be made using Equation 4-27 which gives a value of about 0.5 if 
f - 0.7. We may call this a kinetic efficiency r? K . The quantum effi- 
ciency »jq here approaches unity. Thus the overall efficiency 
r, - nslPJIMQ is about 1 P ercent - If it: ls P° ssible to ^crease the Argon 
pressure sufficiently to raise f to 0.9, then the overall solar effi- 
ciency would approach 2 percent. 

B. The Li 2 Metallic Vapor Laser 

The temperature required to operate the Li 2 laser was reported at 
1000 k( U )- - easily attainable with large solar collectors. The laser was 
also reported as having the lowest threshold power. The criteria of 
Section 3 should be satisfied for Li 2 as well as Na 2 , as the coefficients 
are probably not very different. 

B.l Probability of Filling the Upper Laser Level 

The critical item in the efficiency calculation is the probability 
F - fN (Eq_ 4-26). The energy level diagram for the X 2^ and the B 11^ 
states of Li 2 ( 13 ) shows the spacing of the vibrational levels is greater 
than Na 2 . At 1000 K there would be one percent of the atoms in the v" - 7 
level and transitions up to about v' - 30 are possible. 
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The spacing of the vibrational levels in the upper state is about 250 

cm‘l • 360 K or approximately half the gas temperature, and hence Av' here 

2 

is about 1. Hence, N - (v' - v^) . The cross section for relaxation from 
vibrational level v' by collision is the same order of magnitude for Li2 as 
for Na2 (Section 4. 1.2.4 C below). The rates and cross sections for Li 2 
are not known, but should not be greatly different from Na 2 « 

B.2 Efficiency of the Li2 Vapor Laser 

At 1000 K vibrational level v" — 7 will have one percent probability 
of being occupied and the absorption bandwidth is from 400 to 490 nm. The 
cross sections and rate coefficients for Li2 are not known, but it is 
likely that the overall efficiency would be of the same order as that of 
Na 2 . 


C. Probability of Vibrational to Translational Energy Transfer 
The number of collisions Z required to deexcite a vibrational level 
n to a vibrational level n ± An is given by(l^) 


Z - Z Z Z _ (4-28) 

o osc tr 7 

where Z Q is a geometrical factor of order 3, Z osc is dependent on the 
vibration of the molecule and Zj- r or the translational energy of the 
impacting particle. Then for a molecule B C Struck by an atom A: 


osc 


2m b m c 

2 2 
m b + m c 


m a + m b + M c 1 


(4-29) 


m a 


2 * 


where the M's are the respective masses. The quantities $ and 6' are 
characteristic temperatures defined by 


k 9 - hi/ * he (4-30) 

where k, h, c have their usual meaning and w e corresponds to the 
vibrational energy in cm"^. 
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The temperature 6 ' is defined by 

k #' - 16 M v 2 L 2 (4-31) 

where M is the reduced mass , v - v e c the vibrational frequency , and L an 
interaction distance usually taken as 0.2 to 0.3A. 


The quantity Z^ r defined 



and T is the gas temperature in degrees Kelvin. 

For the Bln u state of Na 2 with u e - 117.3, L - 0.2A, then 9 - 168. 3K, 

9' - 2.316 x 10 4 K. 

For T - 800 K, then Z osc -7, Z tr - 3.5 x 10' 2 which gives Z - 0.7 or 
on average all collisions with Argon deexcited the Na2 • 

For the B 1 ^ state of Li 2 with w e - 255.5, L - 0.2A then 9 - 367. 9K 
and 9' - 1.552 x 10 4 K. Again taking Argon and T - 1000 K, then Z osc - 
2.9, and Z tr - 8.5 x 10 2 . Hence, Z « 0.7 and approximately every collision 

deexcites the Li2* 

4. 1.2. 5 Summary of Non-Dissociative Lasers 

A class of non-dissociative metallic vapor dimer lasers was investi- 
gated as potential solar energy converters. The high temperatures required 
for vaporization cause vibrational levels up to about v" - 10 to be occu- 
pied in the lower electronic state of Na 2 . This allows broad band 
absorption essential for high efficiency and yields a wide range of levels 
(up to v' - 30) to be occupied in the upper electronic state. The upper 
laser level is at or near v^ - 0, and the probability of the molecule 
relaxing to v^ from higher v* is small because it has to do so in a number 
of collisions, for each of which there is the possibility of being quenched 
to the lower electronic state. Inclusion of Argon is essential to increase 
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the collisional relaxation rate to overcome quenching. In laser-pumped 
experiments , the upper limit of Argon pressure in Na 2 was about 70 torr, in 
which case the ratio of relaxation to total (including quenching) col- 
lisions would be 0.7. This means that for Na 2 , essentially all levels 
above v — 5 would contribute little to the lasing power, and that portion 
of the solar spectrum which caused transitions to v' > 5 was essentially 
unused. For levels with v' < 5 approximately half the collisions resulted 
in relaxation to level v^. Even so, the overall solar efficiency was 1 
percent, a number which compares favorably with 1.2 percent for an IBr 
solar pumped laser. ( 12 ) If Argon pressures greater than 70 torr are pos- 
sible, the efficiency would be increased. The numbers for Li2 are not 
known, but it is probably its efficiency would not be greatly different. 

In addition to low efficiency, calculations by W.E. Meador of NASA/Langley 
have suggested that the threshold pumping levels would be much too high for 
practical use. 

4.2 Lasers Which Absorb All the Solar Spectrum: 

Indirect or Blackbody Lasers 

Indirect solar pumped lasers are often called blackbody lasers. The 
concept is that the total sun solar spectrum will be captured and changed 
into heat. A high temperature cavity will be created and thermal radiation 
from the blackbody cavity will be used as the pumping radiation for the 
lasers. Whatever bandwidth is absorbed by the lasing medium, the "hole" 
left in the spectrum, is then continually replenished by the blackbody 
radiation in the cavity. Thus, heat energy at wavelengths other than the 
absorbing wavelength Is transferred into this part of the spectrum, and the 
whole spectrum is utilized. 

Blackbody lasers have several advantages. First, in space, with 
appropriate thermal energy storage, they can operate in the earth's shadow 
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throughout an entire orbit. Second, they do not depend on the spectrum of 
the incident light, but only on the total quantity of energy which is 
degraded into thermal energy. Thus, chemical, electrical or nuclear energy 
can be used. Third, because the whole solar spectrum is utilized, the size 
of the collector-concentrator is reduced. 

There are two types of blackbody lasers, the cavity type and the gas 
transfer type. In a cavity laser, the laser is enclosed within the black- 
body cavity. Thus the lasant is subjected directly to the thermal radi- 
ation pumping. In the transfer laser, two gases are involved. The first 
of these, the transfer gas, flows through the blackbody cavity and becomes 
vibrationally excited by the blackbody radiation. Then the transfer gas is 
mixed with the lasant. Energy transfer occurs and the lasant emits radi- 
ation in the laser cavity. In 1982, the first in-cavity blackbody C0 2 
laser was achieved by Professor Walter Christiansen at the University of 
Washington. In 1984, he achieved blackbody lasing with N 2 0 as the lasant. 
Also, in 1984, a model for comparison with the experiment was developed at 
the University of Washington to explain the results of some of the recent 
experimental work. ( 16 > The author also developed a theoretical model for 
the in- cavity C0 2 laser in 1984. 

In 1982, Langley developed its own transfer laser program starting 
with N 2 into C0 2 . In 1985, the emphasis was shifted to other transfer- 
gas/lasant combinations particularly C0/C0 2 and C0/N 2 0. Also in 1985, 
mathematical models of the C0/C0 2 and N 2 /C0 2 lasers were initiated by the 
author at Langley, but the results were inconclusive as both the values and 
the variation of the rate coefficients with temperature were not well 
known. The portion of Figure 2 on page 10 starting with gas dynamic lasers 
was therefore studied only in a very preliminary fashion and no concrete 
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results were obtained. The work on both types of blackbody lasers was 
dropped in 1988 when the author pointed out that there was a serious 
limitation to all blackbody lasers. The limitations arose because the 
energy density in a cavity, determined by Planck's radiation formula, is 
limited. To get sufficient energy flow for high power would require the 
surface area of the lasing medium to be too large for practical use in high 
power lasers of outputs of order 100 KW to 1 MW. However the theoretical 
work carried out on direct blackbody pumped lasers could be useful in 
future investigations, especially of lower power lasers. 

In Section 4.2.1, some general considerations for designing cavity 
blackbody lasers are outlined, and a calculation is made of the dimensions 
of the surface area of the laser medium for a 1MW laser. One factor in 
estimates of the power absorbed from blackbody radiation is that for 
certain geometries, e.g., slab or cylinder, illuminated from the outside 
the oblique rays travel along much greater path lengths than those normal 
to the surface. A method of calculating the effect was published in Optics 
Communications and is shown as Appendix F. An actual CO 2 cavity laser 

is studied theoretically in Appendix G. Both steady and non- steady state 
outputs are illustrated and rough estimates of efficiency are given. The 
limitation on blackbody lasers due to the energy density and power flow is 
discussed in Section 4.2.2. Section 4.2.2, F gives overall conclusions for 
cavity blackbody lasers. 
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4.2.1 General Considerations f o r Designing C»vttY gjackbody Lakeys** 


A. Introduction 

The purpose of this study is to review briefly the physical processes 
which occur in a blackbody cavity laser and show how the characteristics of 
the lasing medium affect the efficiency when used as a solar energy 
converter. The characteristics of an "ideal" lasant will then be 
described. The blackbody concept^®) has been outlined by Dr. Walter 
Christiansen in the previous talk. 

B. Overall Considerations of Lasing Media 

In general, high laser efficiency requires a lasing medium of the 
molecular rather than atomic variety. The reason is that the laser effi- 
ciency depends on the "quantum efficiency" tjq . 

,q - (E u - ( ‘- 33) 

where E u and are the energies of the upper and lower laser levels, 
respectively, Each absorbed photon gives up energy E u , while the emitted 
photon possesses energy E u - Ej., and to recover the absorbed energy ijq 
should be high. Figure 13 compares atomic with molecular systems. For the 
atomic system shown »jq - 1/6. Also shown is a range of energies for molec- 
ular dissociation leading to an excited atom X* which causing lasing. If 
5eV were needed for dissociation, then as the energy of X* is usually 0.1 
leV, then qq < 0.2. Molecular levels are a few tenths of an eV above 
ground, but although both E u and E x are now quite low, rj Q is much higher, 
typically around 0.5 as shown. For C0 2 , which we shall take as an example 
for the rest of this talk, tjq - 0.43. 


**Talk given in the Workshop on Lasant Materials for Black -body Pumped 
Lasers at NASA Langley, September 1985. 
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c. The Blackbody Cavity 

The intensity of radiation in a blackbody cavity is given by Planck's 
radiation law and the radiation flux, E(A, T b ) per unit wavelength vs. A 
for different cavity temperatures T b are given in Figure 14. Also shown 
are the location of the absorption line for C0 2 at 4.256 pa, and a possible 
dissociation wavelength at 0.5 M m (2eV for Br 2 ). The peak wavelength A m 
and blackbody temperature T b are connected by Wien's displacement law: 

A m T b - 0.29 (cm K) . An ideal lasant should absorb at a wavelength A a which 
approaches A m so that the increased radiation intensity would give higher 
laser output and also E u and »7 q would be higher. 

The blackbody cavity would re -radiate out through the orifice where 
the Sun's radiation entered (Prevost's theory of exchanges), leading to a 

cavity efficiency rj c :^ 18 > 19 ^ 

r, c - l-4(T b /5785) 4 (4-34) 

where 5785 is the Sun's temperature. A plot of tj c vs. T b (Fig. 15) indi- 
cates that T b should not be much greater than 2000 K, a value also predi- 
cated from material considerations. 


D. The "Source Term" for the Upper Laser Level 

The source term is the rate of pumping the upper laser level. It 
depends on the energy flux in the cavity, E(A a , T b )dA a (W-cnT 2 ) where A a is 
the absorption wavelength and dA a the absorption bandwidth. This energy 
flux enters the laser medium and the fraction F per unit volume absorbed 


between a depth Xi and X 2 below the surface is 

F - (exp(-(C0 2 ) o & Xj) - exp(-(C0 2 ) o >/< X 2 ' X i> (4-35) 

where (C0 2 ) is the density of the absorbing medium, and a a its absorption 
cross section: 


a - A 4 A_ /4* 2 c dA 
a a Gu a 


(4-36) 


49 



CAVITY EFFICIENCY V VS Bug- ROPY TEMP 



Figure 15 
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where Aq u is the Einstein coefficient for the transition ground to upper 
level, and c is the velocity of light. 

The number of CO 2 molecules in the upper laser level (the asymmetric 
stretching mode, C 02 ( 001 )) produced per cc per sec is 

S " E(A a’V dA a F A a /hc (4 ' 37) 

where h is Planck's constant and dividing by hc/A a converts the energy 
flux absorbed into frequency of events producing the upper laser level per 
cc . 


The characteristics of the lasing material determining S are A a , dA a 
and Ag u , but S also depends on the density of the absorbing medium. The 
absorbing wavelength A a should be chosen to approach A m (Fig. 14) . The 
value of o a is not too important, as S depends on (CO 2 ) <7 a --and a low a a 
can always be compensated for by increased CO 2 pressure (similarity law). 
Too high a value of the product (CO 2 ) a a results in all the radiation being 
absorbed at the outer edges of the lasing medium near the walls, and , X 2 
and pressure should be chosen to maximize S. The absorption bandwidth dA a 
should be high to increase E(A a )dA a , but high dA a results in a lower a a . 

The bandwidth dA a for CO 2 absorbing at 4.256 /im is actually the sum of 
the bandwidths of about 30 closely spaced rotational levels near the asym- 
metric stretching mode vibrational level (Fig. 16). Each rotational level 
has a Doppler bandwidth dA^: 


^ . 7.2 x ^ , 

A M 


1/2 


(4-37) 


a 

where Tq Is the CO 2 gas temperature (=360 K) and M the molecular weight 
of 44. The Doppler widths could be increased with an ideal lasant of low M 
operating at high gas temperature. Increasing the bandwidths by pressure 
broadening might not be feasible if increased collisions depopulated the 
upper level. 
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The spread of rotational lines from 10.1 to 10.8 /im in Figure 16 
increases with gas temperature. The effective dA a would increase if there 
were more lines, or the spacing of the rotational levels could be 
decreased. The spacing in frequency is 

2 

At/ - h/4*r I (4-38) 

where h is Planck's constant and I the moment of inertia of the mole- 
cule. The ideal lasant should also have a low moment of inertia. 

E. Kinetics of the CO 2 Blackbody Cavity Laser 

The lasant must not be at high temperature or the lower level will be 
filled as will be shown shortly, and as helium has excellent heat conduc- 
tivity, it is usually introduced. A buffer gas such as Argon may also be 
introduced, which if possible should enhance losses from the lower laser 
level without affecting the upper level. The ideal lasant may therefore be 
a mixture of gases, making the choice still more complex. In addition to 
making the source term S as high as possible, an ideal lasant should have 
a high density of upper level states N u , and a low density of lower level 
states N^, to give as high an inversion population (N u - N^) as possible. 
(Fig. 17) A specific example of the processes determining N u and in a 
C02 - He-Ar mixture is illustrated in Figure 18. 

The ideal lasant should have low cross sections for collisions that 
deexcite N u , i.e. low rate coefficients k^, k2, k 3 , and a low Einstein 
coefficient Aq u , from the upper level to ground. A low value of A q u means 
a metastable upper state. Even more important, the filling of the lower 
level (Figs. 17, 18) from the upper state should be small, as this has a 
bigger effect on N u - N]_. (Fig. 18- -ky, kq, k^2 # ^13 should be small.) 
However, collisions which remove the lower level (k^o, ^n) are helpful. 
Losses of all levels to the tube wall occur by gaseous diffusion, which is 
important at low gas pressures (few torr total pressure). 
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Figure 18. Flow Diagram for C0 2 Black Body Laser 



F. Temperature Criterion 

Experimentally, it was found that if the C0 2 lasant became hot, lasing 
was inhibited. Most of the energy absorbed is deposited as heat in the 
lasant and must be removed by diffusion to the walls. In gaseous mixtures, 
the coefficient of heat conduction is independent of the total pressure, 
but depends on the gas ratios, with helium the best heat conductor. 

Assuming planar geometry the temperature of the lasant in the center of the 
laser can be calculated, and for C0 2 -He-Ar mixtures, values around 360-400 
K were obtained. 

Increasing the gas temperature Tg fills the lower laser level by 
Boltzmann statistics, and Ej. is only 0.165 eV above ground for C0 2 . In the 
steady state the fraction of C0 2 molecules with energies this high is 
exp(-E 1 /kT G ) - exp( -1972/Tg) or 0.5 percent at Tg - 360 K. Collisions with 
molecules whose energy is above can result in the symmetric stretching 
mode oscillation- -the C0 2 (100) lower lasing level. Only one in several 
thousand collisions accomplishes this, but nevertheless at elevated gas 
temperatures, the density of the lower laser level Ni state can exceed that 
of the upper laser level N u , and lasing will cease. The lasant must there- 
fore be cooled and the cooling requirement is yet another reason why an 
important characteristic of an ideal lasant would be to have E 1 as high as 
possible above ground. 

G. Characteristics of an Ideal Lasant 

The lasing medium X is usually mixed with a cooling agent Y, (usually 
He) and a buffer gas Z which preferable deexcites the lower level. All 
three components interact in the kinetics of the laser. The ideal lasant 
should satisfy the following criteria: (1) Number of atoms in the mole- 
cule: the lasing medium should have at least three atoms so that there are 
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several energy levels to provide an upper level, and a lower level which 
cannot be the ground level. (2) Absorption cross section a a : it should 
have an absorption cross section o a > 0 at an appropriate absorption wave- 
length, A a . The cross section need not be large as increasing the pressure 
can compensate for small a a . (3) Absorption bandwidth dA a : this should be 

as large as possible (although a a oc (dA a )’^--see 2 above), as a greater 
portion of the blackbody spectrum is now utilized. An ideal lasant would 
have large dA a if it had a small molecular weight M, and a large moment of 
inertia I. (4) The laser energy levels E u , E]_: E u should be high so A a is 

small, and approaching the peak of the blackbody curves, where the radi- 
ation density per unit wavelength is higher and raising rj q - (E u -E^)/E u . 

The lower laser level, E^ , should be high so that the level is not filled 
by temperature effects. However, both E u and must give a reasonable ^q. 
(5) Pumping the upper laser level: absorption in the lasant X at a density 

o 

(X) cm“ J , depends on the correct ((X) o a d) , where d is the average depth 
(between X]_ and X2) and d should be small, to enable adequate heat 
conduction to the walls. (6) Losses from the upper laser level: the ideal 

lasant X should have low losses from its upper level X u , meaning collisions 
of the type X u + X,Y,Z should not deexcite (see Fig. 18). Spontaneous 
emission from X u should also be small, or Einstein coefficient A u g should 
be small, i.e., a metastable upper level is desirable. (7) Stimulated 
emission cross section a e : should be large for low threshold, high gain, 

and high power: 

a e " X e A uL /4 * 2cdA e (4 ’ 39) 

Here A e , the lasing wavelength is fixed by E u -Ei, A uL the Einstein coeffi- 
cient for spontaneous emission should be large and dA e the emission 
bandwidth should be small. (8) Filling the lower laser level: the rate 
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should be small. The two main contributions are collisions from the upper 

level, requiring rate coefficients such as k 7 , k 12 , M3 < Fi S- 18 > to be 
small, and filling by temperature effects should be small meaning should 
be high, and the gas temperature T G low. (9) Losses from the lower laser 
level: should be large so that the density of the lower level states (X L ) 
is small and (Xu) - (X L ) » 0. Hence collisions of the type X L + X, Y, Z 
should have large deexcitation cross sections (large rate coefficients K 4 , 
K 5 , K 6 in Fig. 18). It is also helpful to have other vibrational modes at 
energies in near resonance with El, which provide pathways to deplete the 
lower laser level, e.g., for C0 2 , we have (100)^(020)^(010)^(000). 

H. Scaling of a Blackbody Cavity C0 2 Laser to 1 MW 
In scaling a C0 2 blackbody cavity laser to 1 MW there is a serious 
limitation due to Planck’s Radiation Law, which leads to excessively large 
surface areas for the laser. If a flat "box- type" laser is assumed, the 
surface areas of the solar collector, and heat radiator are considerably 

smaller . 

H.l Power Density in a Blackbody Cavity 

The C0 2 laser absorbs at a wavelength A a - 4.256 /.m, with a bandwidth 
dA - 4 x lO'Vn. Planck's radiation formula then gives the power entering 
unit area of the laser surface as 

E ( x )« - 107 *1^: — (w-»' 2 ) < 4 - 40 > 

a a exp(3381/T b ) - 1 

where T b is the cavity temperature. For T b - 1500K, and 3000K, respec- 
tively, E(A a )dA a - 1.26 and 5.14 mW per square cm, values which are very 
small. Large power outputs therefore require a large laser surface A L . 
There is no concentrating of input flux inside the blackbody cavity, as 
with parabolic collectors. 
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We define a device efficiency ij as P out /Pin where P in is the absorbed 
power to produce C02<001) the upper laser level, and P out is the output 
power assumed to be 1 MW. 

H.2 Calculation of the Dimensions of a 1 MW Laser 

The model is shown in Figure 19. A parabolic collector of area A c 
collects the sun’s radiation of 1.4 KWm" 2 and the total power P^, is 
directed into the cavity, which has negligible losses. All of P c is 
assumed to be input power P^ n to the laser of surface area A^,. Although 
the cavity is at temperature T b , the temperature of the gas T G is assumed 
much lower - 400K. A cooling system removes heat to a radiator of area A r , 
at temperature T r which emits power P r by Stefan's Law: P r - A r £aT r 4 where 
c the emissivity is assumed unity, a is Stefan's constant, and T r - Tq. As 
ij is usually much less than 1, P r = P^ n . With these assumptions the areas 
of the collector A c and radiator surface A r can be calculated as shown in 
Fig. 20. 

In this simplified treatment A c , A^, A r are proportional to P ou t/*?, 
only A l depends on Planck's formula. Assuming cavity temperatures T b - 
1500K and 3000K, and efficiencies of 1% and 10%, the values of the areas 
are given in Figure 21. Also given is a representative linear dimension L 
where L 2 - area, to give a feeling for the size of these components. The 
radiator may to some extent be folded, while the laser area, Al, can be 
shared by the two sides of the laser. However, a folded geometry for Al 
would be difficult because the surface has to be cooled as shown. Even in 
the extreme limit with T b - 3000K and rf - 1, A L is still 20000 m 2 if P out - 
1 MW. 
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I. Conclusions 


The characteristics of an ideal blackbody cavity lasant have been 
outlined. The choice of an ideal lasant is a complex process depending on 

a large number of factors, including the choice of a cooling medium and a 
buffer gas. 

Planck's radiation law limits the power input per unit area into a CO 2 
blackbody cavity laser making the surface area for high powered lasers 
excessively large. It is suggested that an alternative application might 
be small 1 W lasers for communication and surveillance, because It would be 
easy to maintain the cavity temperatures in synchronous orbits where 72 
minutes each day are spent in the earth's shadow. 

4*2.2 A Limit ation in Scaling up a Blackbody Laser 

A. Summary 

A blackbody pumped CO 2 + buffer gas laser is assumed cylindrical of 
length L and radius r. The pumping rate is determined by the rate of 
arrival of radiation in a bandwidth 0 . Ifi at the laser surface according to 
Planck s radiation formula. The power input per unit volume however is 
proportional to 1/r. A maximum value of r is determined from experi- 
mental results for threshold at 1500 K. Increased power can only be 
achieved by increasing L which becomes impossibly large. The limitation 
can be circumvented by a box type laser. 

B. Introduction 

The purpose of this section is to illustrate a limitation on the pump 
power input into a blackbody laser if the geometry is cylindrical. The 
input energy flux is limited by Planck's radiation formula, and there is no 
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means of concentrating it, as would be done by parabolic solar reflectors. 
The pump power per unit volume is estimated from the experimental results 
of Insuik et al.( 16 ' 20 ) To attain a given power input per unit volume, it 
will be shown that there is an upper limit to the radius. Hence in scaling 
up to large power output, a cylindrical laser turns out to be impossibly 
long. A method of overcoming the limitation using a box laser geometry is 

given. 


C. Absorption of Blackbody Radiation 

Experiments have shown that for a C0 2 blackbody laser, with a buffer 
gas, absorption occurs at 4.256 n with a bandwidth of 0.1 m- ( 18) The lasing 
medium in the form of a cylinder of radius r and length L receives 
energy from the blackbody cavity at a rate per unit surface area given by 
Planck's radiation formula. 


E(A) dA - 


3.74 x 10 


-12 


.5, 1.439 

A (exp 


AT 


dA 

- 1 ) 


-2 

watts cm 


(4-41) 


where A is the wavelength, dA the absorption bandwidth both in cm and T 
the temperature in degrees Kelvin. For the above wavelength and 


bandwidth : 


E( A) dA 


2.68 


, 3381 
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T 


-2 

watts cm 

-1 


(4-42) 


a function of temperature only. 

The power entering the laser medium is then E(A)»2xrLdA, but to obtain 
the power P entering unit volume on the average we must divide by 

or 

P - 2E( A)dA/r (4-43) 

Hence, increasing r, decreases P, and there is a limitation on r, 
which has a maximum value given by equation (4-43) when P is just 
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su ^^ c ^ ent f° r threshold. To scale up for high power output, (for a given 

efficiency) the only remaining alternative is to increase L. 

Some illustrative values of E(A)dA and P are shown in Table III, 

where r is assumed to be 1 cm. 

Table III. Values of E(A)dA, the power per square cm 
and P, the power per cc vs . temperature 


Blackbody 

E(A)dA 

P 

temp K 

watts cm’ 2 

watts cm'^ 

1000 

0.094 

0.188 

1100 

0.130 

0.260 

1200 

0.170 

0.340 

1500 

0.314 

0.628 

2000 

0.606 

1.212 

3000 

1.285 

2.569 


D, Estimate of Threshold Input Power 

An approximate estimate of the pump power per cc to attain threshold 
can be obtained from the experimental results of Insuik et al.( 2 °) Their 
experiment was done in a sapphire tube of internal diameter 8 mm, outside 
diameter of 9.5 mm, with 50 cm of its length in a hot zone. The medium was 
co 2 with He and Ar as buffer gases. Lasing occurred with the blackbody 
temperatures ranging from about 1100 to 1500 K. 

By varying the pressure, it was found that a drop off in power 
occurred around 7 torr total pressure when there was 21% CO 2 , suggesting 
that the pressure and path length were just great enough to absorb all the 
input radiation. Assuming the absorption length L* to be of order 0.8 cm, 
this leads to an absorption cross section « 3 x 10’ 17 cm 2 and a relation 
between the path length and CO 2 pressure: 

L' p a 1 (cm torr) (4-44) 

where L - 2 r is in cm, and p in torr for maximum output. Neither the 
pressure of CO 2 nor r can be increased indefinitely in scaling to higher 
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power, otherwise the radiation will not penetrate to the center of the 
tube. If Lp is too small, the radiation will not be absorbed. 

A very rough estimate of the pump power per cm3 for threshold can be 
made. A number of experiments < 2 °) showed that maximum output power 
occurred over a narrow range of total pressures P t ; 2.5 < P t < 12.5, and 
C0 2 pressure p; 0.5<p<2.5 torr. A laser power output vs. temperature 
curve for P t - 7 torr, p - 1 torr showed that threshold occurred around 
1100° and the power output rose linearly to a value of 4 mW at 1500 K. 
Assuming 1100 K to be the threshold temperature for r - 0.4 cm, then 
equation (4-43) gives the minimum power per cm 3 for threshold to be P min - 
0.65 watts cm* 3. 

Now P min may not be a constant value if r is varied; increasing r 

decreases wall losses. Assuming it is correct in order of magnitude, then 

the maximum value of r is obtained from equation (4-43) 

r - 2E(A)dA/P , (4 ' 45) 

max ^in 

With P mln - 0.65 watts cm-3, a.,, at T - 1100* x m - 0.4 cm, and at 1500‘C 
r max “ 1 cnu 

E. Power Output and Laser Dimensions 

If the overall efficiency of the laser is t/, then the power output P Q 
is 

p - n 2irrL E(A)dA (4-46) 

o 

and increases with rL for a given cavity temperature T. To increase P 0 , r 

should approach r max , determined for given T from Equation (4-46) . Then 

given L can be calculated; e.g. for P Q - 10KW, T - 1500 K, , - 0.1, 

r - 1 cm then L - 500 m for a laser tube of 1 cm radius. The length L 
i max 9 

is impossibly large. 
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The difficulty can be overcome by a box laser arrangement (Fig. 22) 
when the medium is in the form of a flat plate. The output power is now 

P 0 - 2 t) xy E(A)dA (4-47) 

the factor 2 arising as it is pumped from two sides. Using the same 
numbers as above, then xy — 16 m^ , or x — y — 4m, much more reasonable 
values . 

F. Conclusions 

In scaling up a blackbody cavity laser to high powers, there is an 
inherent limitation because the energy flow into the laser medium is 
limited by the Planck radiation formula. Concentrating it by parabolic 
reflectors does not seem possible. In addition, the power input per unit 
volume varies inversely as the radius if the medium is cylindrical, and the 
value of the radius is limited by a threshold condition. For high power, 
the length L must therefore be increased to a value impossibly large. 

The limitation can be overcome by using a box laser arrangement, which 
essentially offers a greater area for pumping in a much more compact 
arrangement . 


4.3 Dye Lasers as Solar Energy Converters 
In 1982, a preliminary survey of Dye Lasers for solar power conversion 
was made. The reason for the survey was that the efficiencies of the 
systems considered hitherto were found to be very low. In brief, the 
previous investigations had first considered different absorbing (Br 2 ) and 
lasing media (CO 2 , HF , HCN, H 2 O) and the overall efficiency for all combi- 
nations was less than 0 . 5 %. (21, Appendix C) Investigation of an IBr laser 
which uses the same medium for absorption and lasing and where the 
"transfer efficiency" was largely eliminated showed the final efficiency 
1.2% or 1 ess.( 15 > Appendix D) 
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The new types considered were (1) liquid dye lasers and (2) vapor dye 

laS 6IS ♦ 

A. 3.1 Liquid Dye Lasers as Solar Energy Converters 

An experimental investigation of liquid dye lasers was carried out by 
M.D. Williams of the National Aeronautics and Space Administration in 1982. 

In general, such lasers have high gain and operate at room tempera- 
ture under giant pulse pumping. For solar pumping in the steady state a 
number of difficulties arise. 

Lasing in organic dyes differs from one material to the other, but in 
general, the molecules possess a pair of electrons which possess a certain 
freedom of motion within the molecule. The motion of this electron pair 
determines the electronic configuration of the molecule. The lowest elec- 
tronic configuration Sq is a singlet state (spins opposite). There are 
excited singlet (S lt S 2 , . . . ) and triplet (T^ , T 2 , . . . ) states. Transitions 
between singlet states are "allowed” by the spectroscopic selection rules. 
They give rise to intense absorption and emission spectra. Transitions 
between different triplet states are also allowed, but singlet-triplet 
transitions, known as inter-system crossings , are forbidden. 

The presence of molecules in one of the metastable triplet levels may 
be detrimental for laser action in several ways. First, the molecules 
stored in a triplet level are eliminated from participating in the laser 
cycle. Second, when some of the lowest triplet levels are populated, the 
material may become absorbent for the laser radiation because of triplet - 
triplet transitions that overlap the laser line. 

It is possible to quench the triplet states, for example, by satu- 
rating the solution with oxygen, and pulses of 500 ps duration have been 
achieved by Schmidt < 22 ) in air saturated rhodamine 6G solution. Commerical 
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cw dye lasers are now available, but they have to be pumped by ultraviolet 
lasers. An excellent compendium on cw Dye Lasers is given by Snavely 
and Schafer(^) . 

Another problem is the inducement of optical inhomogeneities 
(differences in reproactive index) caused by the heating and non-uniform 
excitation of the solution. 

In conclusion, the quenching of triplet states, the presence of 
optical inhomogenieties , and the problem of cooling the liquid at high 
power levels to prevent boiling would raise severe problems for cw solar 
pumped liquid dye lasers. 

4.3.2 Vapor Phase Dye Lasers 

The use of vapors is advantageous for solar pumping because of struc- 
tural simplicity and uniformity of medium. Recently, a paper by Basov et 
al.(25) reported experiments with vapor phase dye lasers. The results 
enable some tentative conclusions to be drawn on whether they would be 
feasible for solar pumping. 

The pumping was performed with the third harmonic of a neodymium laser 
at 355 ran. This wavelength lies on the low wavelength side of the solar 
spectrum peak which is at 500 ran, and the number of photons per unit wave- 
length is about half that at the peak. The absorption bandwidths of dyes 
in liquids are typically of order 100 ran, and if this is also true in the 
vapor phase, then the fraction of photons in the solar spectrum absorbed 
the solar efficiency- -would be somewhere in the region of 3%. 

Values of efficiencies of the vapor lasers were quoted, but these 
numbers did not include a "solar efficiency." The highest "efficiency" for 
a vapor phase laser was for the material POPOP, which is 
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p-bis [2 - (5-phenyloxalyl) ] benzene and was 23%. Inclusion of the solar 
efficiency shown gives an overall efficiency of the order of 1%. An 
efficiency of over twice the above value was quoted for POPOP in liquid 
ether solution, a liquid dye laser. It would be unlikely that temperatures 

sufficiently low to keep the ether in the liquid state would be possible in 
a solar laser. 

As steady state working would be desirable for solar lasing, the 
threshold pumping intensity is important. The lowest value quoted was 20 
kW cm' 2 which would correspond to the sun's radiation being concentrated 
1.4 x 10 5 times. However, the length L of the active region of the laser 
was only 2.4 cm. If it is true, as in the IBr laser, that the product CL, 
where C is the concentration factor, is roughly constant, then C would 
be about 3 x 10 3 for L - lm. 

In conclusion, the threshold pumping intensity for such a vapor laser 
would be within reasonable bounds, and its overall efficiency would be 
comparable to IBr. It is noted that the experiments were performed with 
ether as the buffer gas at 36 atmospheres. 

4.3.3 Conclusions on Dye Lasers as Solar Energy Converters 

It was concluded that liquid dye lasers would probably not be 
feasible. Vapor phase dye lasers offer some possibility although their 
overall efficiency would not be much greater than IBr at around 1%. 


5 . OVERALL CONCLUSIONS 

Solar pumped lasers used as energy converters have many important 
applications for NASA's space program. They can be used to beam power from 


one space vehicle in a mission to other more maneuverable vehicles, or else 
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from orbiting vehicles to rovers on the surface of planets. Low powered 
solar pumped lasers (in space or terrestrial) can also be used for 
communication and surveillance without the need for auxiliary power 
sources. An important criterion from the space program's point of view is 
high output power/weight (Section 3, Appendix A). The main weight woul 
that of the heat radiator, and its size would be reduced if the overall 
efficiency of the laser would be maximized, and also if the laser could run 
at high temperature. Studies of the far field patterns indicate so far 
that if the power outputs are high, then the area of cross section of the 
lasers will be large, giving rise to modes at right angles to the axis with 
resulting spread of the beams. This aspect needs further study. 

The mechanisms of the lasers and the various families of solar pumped 
lasers are shown in Fig. 2 (p. 10). The class of lasers utilizing only a 
band of the solar spectrum were studied in detail. The dissociative types 
were divided into two, namely where one material absorbs and hands over 
energy to another which lases and second where the same material absorbs 
and lases. A detailed study of the first type in the form of Br 2 -C0 2 -He 
(absorber, lasant, coolant) Br 2 -H 2 0-He, Br 2 -HCN-He , and I 2 -HF-He showed 
that in all cases, even assuming complete absorption, the overall 
efficiencies were below 5 x 10 ^ . (Appendix C) For the second type--IBr-- 
the estimated efficiency was 1 . 2% . (Appendix D) 

In the case of the non- dissociative lasers, the studies first concen- 
trated on some physical aspects. First, it was shown that stimulated 
emission in Na 2 occurred at somewhat different wavelengths to regions of 
self absorption, a necessary condition for lasing (Section 4. 1.2.1) 

Second, a simple theory was obtained to estimate the ratio of dimers to 
monomers in metallic vapors as a function of temperature (Section 4. 1.2. 2). 
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Third, the question of metastability in N 2 and CO while the vibrational 
levels in Na 2 are short lived was explained by estimating the spring 
constants of the molecules and relating these to the spacing of the levels. 
For Na 2 , the spacing is of order the thermal energy, for N 2 and C0 2 the 
spacing is 15 times greater-hence thermal collisions can deactivate the 
higher vibrational levels in Na 2 , but not N 2 and CO 2 . Fourth, a new 
hypothesis on the electronic energy transfer in Na 2 based on a suggestion 
by N. Jalufka was proposed. (Appendix E) It is based on the blackbody 
radiation in the oven at 750°C causing stimulated emission between the 
electronic bands B and 2 2 . The results have been submitted to the 

o 

Journal of Molecular Spectroscopy. 

Specific metallic vapor phase laser systems in the case of Na 2 and Li 2 
were studied and the overall solar efficiencies were estimated at around 
1%. However it was pointed out by W.E. Meador that the threshold poten- 
tials for these materials as solar pumped lasers would be far too high for 
practical use and it seems unlikely that they would be useful. 

The blackbody lasers have the advantage of thermal storage and can 
operate in the earth's shadow. The size of the collector is reduced as the 
whole solar spectrum is absorbed. A physical problem addressed was the 
effect of isotropy of the radiation for certain geometries of absorber and 
a paper was published in Optics Communications (1?) (Appendix F) . A 
theoretical study of a solar pumped blackbody CO 2 laser outlined the 
physical mechanisms and estimates were made of output power as blackbody 
and gas temperature. The results depended critically on the values of the 
rate coefficients which were not well known, and the agreement with the 
experimental results of Insuik and Christiansen (IEEE J. Quantum 
Electronics QE-20 622, (1984)) was poor. However, the author pointed out 
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that there was a limitation in scaling up blackbody lasers because the 
energy density in the cavity was United by Planck's radiation formula, and 
for high powers, the area of the absorbing surface of the l.sant would have 

to be impossibly large for space applications. 

A brief period was spent studying liquid and vapor phase dye lasers, 

and it was concluded that they would not be feasible. 

The overall picture obtained is that it is possible to obtain lasing 
by solar pumping, and indeed, experiments on IBr, on the perfluoroalkyl 
halides and on black-body lasers have borne this out. So far, for energy 
conversion, the conversion efficiencies which in turn effect the output 
power to weight ratio, seem too low for immediate application and further 
study is needed. The question of far field patterns which was put in 
abeyance also needs further study. 
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Criteria for the Evaluation of Laser 
Solar Energy Converter Systems 


W.L. Harries’ 

Old Dominion University, Norfolk, Virginia 

I. Introduction 

O NE concept for collecting solar energy is to use large solar 
collectors on orbit.ng space stations, which then transmn 
the energy as laser beams. Direct conversion by solar pumped 
gas lasers has already been considered and the ^er conver- 
® on efficiency disscussed. 1 ! However, a figure of mem. of 
more im portance for space missions, is the output power pe 

UI1 An idealized energy conversion system is shown in Fig. '■ 
Tht .ol.r r.ylc. “diced • p.r.bollc nfltajr »M« 
a required output power, its area is determined by the l“ er ^ 

fluency. The power not converted int0,h * ’^'^L^which^ 
Hissmated by a heat radiator emitting by Stefan s Law whicn, 
as the converter efficiency is usually low, will contribute the 
major part of the total weight of the syst * m , d = , 0 

The overall efficiency of the laser can be decomposed into 
the product of several efficiencies, which affect the collector 
IXZror ir, different ..»• The .rad™.,., m decent 
in Sec. II and the energy flow is discussed mSecJIL he 
dependencies of weight on area are then assumed, which deter- 
mine the weight for a given output power. Some comparisons 
will be shown of different laser systems. 

H. Efficiency of Solar Pumped Lasers 

The overall efficiency of solar pumped lasers can . be sub- 
divided into the product of four efficiencies. First, M .abor- 
tion occurs over a bandwidth X, to X*. the fraction of the solar 
spectrum used, or “solar utilization efficiency, is n s . 


Ts 


= p 2 +(X)dX/ *(X)dX 


III. Energy Flow in Solar Pumped Lasers 

The energy flow for the idealized model in Fig. I is shown in 
Fig 1 The solar power falling on the collector of area r is 
p fhe mirror of reflectivity r concentrates the flux, and a 
fraction rP reaches the first transparent wall containing the 
lasing medium. However, a fraction (l -r)P has been ab- 
sorbed bv the mirror and has to be radiated. The first con- 
tainer wall is assumed to have a transmittance r, , so the pow 
arriving at the lasing medium is r,,rP„ and the power absorbed 
bv the wall is (1 - r, )rP , , which also has to be radiated. 

^e power arriving at the laser medium is distributed in 
wavelength over the solar spectrum; only a fraction, is- van i be 
of .n» only e, .. absorbed. The po.ee .hoehed 
p is ncT),r,/-P, and the output power P 0 is 1* • Vkhl J* 
part lost due to the lasing process. ( l - ix-lp )f*4 • has 10 
radiated The power not absorbed by the lasing medium is 
r ^l’i^ and reaches the second wall of transmittance 
r Then r U -irtJttP, is transmitted through it. but an 

? m ou»U I - 1.) »,* absorbed " “ 

radiated. . ■ D 

The total power that has to be radiated is r, 

p =P j(i -r) + r( 1 - r,) + rriij s t?d4(l “Wq) 

+• j ( 1 - Vs^A H l ~ r :M * 

As the output power P Q and input power P } are related by 
Po = rTxy]Pn the relation between the radiated and output 

power is 

( 2 ) 


where + (A) is the solar radiance in photons per square meter 
bC ThT absorption efficiency V U that ^on °f photons 

tion. (AO the density of absorbers, and d the thickness of gas. 

Of the photons absorbed, only a fraction. r, K . end up pro 
ducing a lasing transition. The kinetic effiaen^ nx is deter- 
mined by the detailed kinetics of the laser and, for example, it 
takes accounts of losses in the upper laser level due to quench- 
ing etc.*- 3 Finally, the quantum efficiency is th * 
the energy of the emitted photon to the average energy of the 

ab T°hTdtvucr efficiency would then be iid = VaViOQ‘ al3d the 
overall^solar efficiency ,-,,*.'"9. which is usuaily an 
order of magnitude less than r^ D . 

— ^ a 1QR4- Technical Note received April 22, 1985. 

cSZ Aeronautics and Astronautics. 

Inc., 1985. All rights reserved. 

•Professor, Department of Physics. 


Pr = KP 0 /?\ 


where 


K= ( 1/rrjUl -r) +r(l - r,) + rr x ri 5 r^ A (\ 
- rr^l — Tj) (1 — ) 1 


(3) 


IV. Output Power to Weight Ratio 

, * a f i t u KWm 4 , and if the output 

The solar irradiance / 0 is 1.4 Nwm , 

power is P Q (kW), then the area of the collector A c ( 

A c =Po/UoVrrf <4) 

The radiator emits by Stefan’s Law and its area is 

A r = P,/otT* (5) 

where <r is Stefan’s constant, e the emissivity, and , T the 


parabolic 

rse lector of 

REFLECTIVITY 

r , AREA Ag 5 q^AR . 

RADIANCE t 


laser output 
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TRANSMITTANCE T, 
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radiator of 
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Table 1 

Comparison of various lasing materials in solar nnm^rf i...~ 


Material 


Vk *lQ 

V 

7\K 

P 0 /W 

0/a = 1 $/ a — jo 

0H1GIMAL PAGE IS 

IBr 

C 3 F 7 I 

0.12 1 

10~ 2 1 

0.57 0.18 

l 0.22 

1.2 x 10 “ 2 

2.2 x 10 ^ 

423 

873 

1.2x 10'-* 4.4 x JO 7 " 1 

2.69x I0’ 3 2.48 x 10” 3 

OF POOR QUALITY 


46»orb«j in [7 , IP " ' \ 

y '' Mirror _« ftADUTCO 




JiK- 2 Energy now chan for a solar pumped laser. 


ing temperature of the laser and that A, is just sufficient to 
radiate the power required. 

Assumptions must now be made about the dependence of 
ihe weights of the collector H / , and radiator W on their 
respective areas A, and A r ; here, both are assumed propor- 
tional to their areas. H'.-a/t,. W t .0A„ where a and 0 
(Kgm -) are constants of proportionality. (The latter assump- 
non may only be approximate as. for example, heat pipes may 
be needed to conduct to the radiator.) 

The total weight W of the laser, collector, and radiator is 
W f + W t . + iY t , where W L is the weight of the laser. It is 
not unreasonable to assume W , « 3P, , and, jf so, the out- 
put power to weight ration P Q /W is 



k „ on in *'* *'*** fur different Icmperaiurrs 

showing effect of varying 0 /*. Both jt and y are dimension lew, while 
a and ^ are In units of kg/m’ 2 . Valuev awumed are r-Q.n 
») « U b> 0/« - |0, an extreme cave. 


p o /w = i,/[ (a/!.4rr t ) + (/3AVo«D ]; W L < W (6) 

Assuming that a and which are determined by structure 
design, are constants for all lasers, then Eq. (6) enables a com- 
parison to be made for different lasing materials, provided the 
various efficiencies at the working temperatures are known. 
For fixed or, we may also make comparisons in terms of 0/a. 
The dependence on the efficiencies is P 0 /W=f{ x ,y) t where 
A - 'M i and y-r \^ a plot of P 0 /W on an xy plane is 
shown in Fig. 3, for r=0.98, 7, = 7 , =0.9. (The value 0.9 is 
achievable if the walls are either KCs or Csl. 4 ) The values of 
(i/a are taken as I and 10, the latter being an extreme case for 
illustrative purposes. The range of a* and y is only 0-0.3, cor- 
responding to actual efficiencies in practice; *=0 9 and 
300 < Ts 1500 K. * 

For both cases P a /W is either constant or increases 
monotomcally with both jr and y. The effect of 0/a is large at 
low temperatures, (comparing Fig. 3a and 3b) and here P Q / 
fFal/0. At high temperatures, the T* dependence quickly 
makes the area required for the radiator small, and P Q /W is 
independent of both T and /3/a, provided T is high 

A comparison of IBr and C 3 F 7 J lasers is shown in Table I. 
It is assumed that the gas pressure and depth are sufficiently 
high, if.., * 1 . The kinetic efficiency of C 3 F 7 I is assumed as 1. 
The value of P a /W t is calculated from Eq. (6), assuming 
/■-0.98, r, =r 2 =0.9, a = 1 . The temperatures are approxi- 
mate working temperatures. Table 1 shows that the IBr laser 
has a power/weight ratio several times higher than the C\F 7 I 
laser for 0/a =1 and 10. 


V. Conclusion 

Assuming that a laser solar energy converter has a radiation 
collector and heat emitter whose weights are proportional to 
their areas, and that the weight of the laser is negligible in 
comparison, the output power per unit weight can be ex- 
pressed in terms of the efficiencies and working temperatures 
o the system. This ratio seems to be several times higher for 
Br a C i f ‘7 | lascr because if* is greater for the former, 
although the working temperature is lower. However future 
converters should operate both at high efficiencies and’at high 
temperatures. * 
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angley Research Center, Hampton, Virginia. The author 
acknowledges many useful discussions with Dr. J H Lee of 
NASA Langley, who first pointed out the significance of the 
power-to-weight ratio. 
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FAR FIELD PATTERNS OF HIGH-POWERED LASERS 

By 

Wynford L. Harries* 

I . INTRODUCTION 

The overall purpose of this investigation is to study the feasibility of 
using lasers as solar energy converters, with special emphasis on their role 
in space missions. Examples would be sending power from one space vehicle to 
others, or from an orbiting vehicle to vehicles on the ground 1 . The criteria 
for lasers as solar-energy converters have been given already^ , and high 
efficiency, and high output power per unit weight are important. Apart from 
the above criteria of efficiency, power/weight, and low threshold, another 
problem for high power lasers is the far field pattern, especially where the 
geometry of the lasing medium is large. Emphasis in Langley Research Center 
at the moment is on iodine lasers, and on solid state lasers such as coupled- 
diode arrays^ >^. The theoretical studies of diode arrays have concluded that 
reasonable far field patterns are feasible - 80% collection efficiency into a 
3m diameter dish placed 5 x 10^ Km away. 

In the period 16 January to 15 June 1989, a preliminary study has been 
made of obtaining far field patterns for a continuous lasing medium in the 
form of a right circular cylinder pumped radially. Thus the method of 
pumping from the outside is entirely different to the diode array. For 
either a master oscillator-power amplifier (MOPA) system, or an oscillator, 
we consider the last pass through the medium. The effect of amplification is 
included, and most important is the effect of amplification which varies with 
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radius. The medium could be solid, liquid or gas. The entering wavefront 
amplitudes are assumed planar or Gaussian with constant phase, tailored for 
best mode pattern. Although the discussion assumes a steady state CW laser, 
the concepts here can in general be applied to steady or pulsed lasers. 

In sections II and III, the minimum size of the laser is discussed for a 
given output power, and the difficulties arising from large size. The advan- 
tages of a master oscillator-power amplifier are given in section IV, and 
section V describes further difficulties with large amplifiers. Section VI 
describes the method of analysis, which is based on the method of Fox and Li 5 
for plane parallel resonators. The differences in our treatment are empha- 
sized and preliminary results are shown. Sections VII and VIII give 
conclusions and future work required. 


II MINIMUM CROSS-SECTION FOR A HIGH POWERED LASER/AMPLIFIER 
The power P from a laser of cross-section A is given by 

P - h n hi/c T A (1) 

2 

where n is the density of the lasing photons, h is Planck's constant, c 

the velocity of light and T the transmission factor of the output mirror. 

12 

High values of P require A to be large for a given n, e.g. if n * 10 as in 
a gas laser with hw - 0.44eV and T - 0.05 a power of 1 MW would require A - 
0.7m 2 with a radius - 0.5m. On the other hand, if a one pass amplifier were 
used with T - 1, then for the same n, A would be reduced 20 times and the 
radius would now be about 0.11 m. 

Again the intensity I of the lasing photons (I - n hvc) must be less 
than the saturation intensity I s - h^/2ar where o is the stimulated emission 
cross-section, and r the spontaneous lifetime for the lasing medium 7 . 
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The power output 


P — AI < AI S — A hi/ 
2 or 


( 2 ) 


Again we require A to be large, for large power. 

(Equations 1 and 2 are almost equivalent as — 


c AN where AN is the 


population difference and AN ~ n near saturation.) 


HI. MODES ARISING WITH LARGE RADII 

We consider a laser in form of a cylinder of length L and radius r. The 
larger these values, the greater the number of modes possible because the 
frequency spacing between two consecutive modes, either along the axis or 
perpendicular to it would be 8 


v Al - u - 
n+1 n 


2L' 


(3) 


where L' can be either L axially, or 2r In the plane of cross section. (For 
solid and liquid lasers the refractive index must be included.) These modes 
will have wave numbers k nL axially, k mr radially, and hence, the waves can 
add to give a resultant wave number k^ m given by 8 


k 2 - k 2 + k 2 
run nL mr 


(4) 


Thus the number of modes can be large, and picking one mode only, difficult. 
This difficulty can be overcome by using a master oscillator-power amplifier 
system. 


IV. ADVANTAGES OF A MASTER- OSCILLATOR POWER AMPLIFIER SYSTEM 
The use of a laser amplifier-oscillator combination offers an advantage 
over a stronger laser-oscillator that is the same as in conventional elec- 


tronics. It is possible to design and control the low-power oscillator more 
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precisely than an oscillator of large power -generating capability. A low- 
level laser oscillator may be built with confinement of excitation to rela- 
tively few modes. The amplifier that follows the oscillator need not be 
highly selective, but it must provide amplification in the frequency range 
containing the oscillator output. 

The gain of a laser amplifier is limited because a very long amplifier 
or one very highly excited, may turn itself into a pulse generator. There- 
fore, when a laser system with a very large gain is required, it must be 
built using several amplifiers in series, and the passage of radiation in the 

reverse direction must be inhibited by isolators. 

Without the isolators, if the gain of one stage were G and the 
reflectivity of the window of the next stage y, then instability would occur 
if Gy > 1. A series of amplifiers would require isolator stages and beam 

expanders . 

Also the problems of threshold and amplification of high power are 
relegated to separate units. It can be advantageous to use unstable mirror 
systems which reverse the radiation back through the same amplifier, but here 
our discussion centers on an axisymmetrical cylindrical one pass amplifier, 
the last stage of the series. It is assumed that a plane or Gaussian 
wavefront enters this amplifier and that adequate phase conjugation has been 
attained prior to this final amplifier* - 10 . The sides are considered trans- 
parent to enable penetration of the pumping radiation. 


V. MODE PATTERNS WITH AMPLIFIERS OF LARGE RADIUS 


5.1 Oblique Rays 

If there were no amplification in a medium in the form of a cylinder 
(Fig. 1), then a beam of light entering on the left would be expected to 
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Figure 1. Medium in the form of a cylinder with wavefront entering 
from left. (A) Amplification a constant, ray CD gets 
amplified more than AB as longer path length; (B) a varies 
with radius; ray EF gets amplified less than GH, while 
GJ travles over regions of varying a. 





continue to travel as a straight beam if r » A (Ray Optics) . The beam 
profile could be found by Huygen's Principle, where each point at entry is 
assumed to emit a spherical wavefront and the wavefronts interfere. With 
amplification however, there may be changes to the beam pattern. 

First, with constant gain a, if r/L is appreciable, ray CD will be 
amplified (Fig. la) more than ray AB. as CD passes through a greater length 
of medium. Obliqueness is neglected in the usual axial approximation. 

If a varied with radius, and were greatest near the outside for example, 
then (Fig. lb) the amplitude of ray EF would be less than GH, while ray GJ 
could be intermediate. The value of a would depend on a number of factors, 
one of which is how the pumping frequency is absorbed. 

5 . 2 Absorption of Pum ping Frequency 

The absorption efficiency „ A is the fraction of the pumping radiation in 
the absorbed bandwidth that is deposited in the medium and should be high. 2 
For one pass of pumping radiation, we require that the absorption length A a 
be « 2r, but if so the profile of the intensity of the pumping radiation 
could be much higher at the outer surface of the cylinder than on axis 
radiation. The gain a - o s (N u - N X ) , where a s is the stimulated emission 
cross-section, depends on N u and N X , the densities of the upper and lower 
laser levels. Now the rate of production of N u depends on the pump inten- 
sity, the rates of stimulated emission (proportional to the laser photon 
density), spontaneous emission, the quenching rate and the rate of diffusion 
of upper level states in the medium. If the pump intensity were small on 
axis, then N u and a could be small there also. 

5.3 Multipass Pumping 

A reduction in the variation of N u with radius would be achieved by 
sending the pumping radiation through the medium several times by reflection 
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with the absorption small (low pressure gas) during each pass (Fig. 2). 
However for a number of passes, the beam at each reflection has to pass 
through the outside surface twice, where some loss occurs, and be reflected 
with reflection losses. If the fraction remaining after each reflection is R 
and there are n passes, the absorption efficiency or fraction of radiation 
deposited in the gas is 

00 

rf a - (1 - exp(-yd)) £ R^' 1 ^ exp( - 7 d(n-l) ) ( 5) 

n-1 

where y is the absorption coefficient to create the upper level (y - Na a , 

N is the density of absorbers in the medium, a a the cross-section for absorp- 
tion), and d is one path length through the medium; d =* 2r. A plot of rj a 
versus number of passes n for different R is shown in (Fig. 3) for 7 d - 
0.01 (i.e. absorption length 100 d) . It can be seen even for R - 99% that 
only half the energy is deposited in the gas, however many passes. To 
increase rj a then 7 d must approach unity - in which case then the pumping 
radiation will be mostly absorbed just within the outside surface of the 
cylinder. 

5 . 4 Saturation 

The gain a(r) will also be reduced on axis if a wavefront with a 
Gaussian amplitude profile, with its maximum on axis is transmitted through 
the medium, and conditions are near saturation on axis. Then neglecting 
other factors, a will be reduced to ay/2 on axis, where a u is the unsaturated 
gain. 8 Evidently, the effect of gain varying with radius must be 
investigated. 


VI. METHOD OF ANALYSIS 

To find the field patterns of a cylindrical laser, the method of Fox and 
Li was applied. 5 Their original work investigated a plane-parallel resonator 



UMPING 

DIATION 


; 2. Miltipass absorption where pumping 

radiation is reflected through medii 
several times. Mirror losses occur 
at A and D, surface losses at B 
and C such as absorption in a con 
taining vessel wall if the medium 
is a gas. 
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with the radiation reflected many times, and showed how modes would arise. 

Here we consider one pass only through an amplifier where the initial ampli- 
tude is either planar or a Gaussian, produced by the previous amplifier. 
Thereafter, it proceeds onwards without amplification to form a far field 
pattern. The treatment in this report so far is for transmission through the 
amplifier medium only. The present method differs from the original in 
several ways. Here the exact Kirchoff integral is used, and not the paraxial 
treatment, as large radii would ensure offaxis rays. Amplification in the 
medium is also included. 

The electric field distribution is assumed known at plane 1 (Fig. 4) 
and is represented by a quantity U^X]^) where U L consists of two parts 

U r i and U^i, real and imaginary. The quantity yf^i + u n 1S then the 
amplitude, and the phase is tan' 1 (U U /U r i). The amplitude squared is pro- 
portional to the energy density, and when integrated over the surface is 
proportional to the total power. To obtain the distribution on plane 2, the 

Kirchoff integral is used . 5 

IL (x. , y. )exp(ikr) (l+cosfl)dx 1 dy 1 

1 1 1 - (6) 

R 

where R is the distance between Pi and P 2 . 9 is the angle P^ makes with the 
normal to the surface at Pi, A is the wavelength of the radiation; k the wave 
number and dxidyi is the element of surface around Pi- By writing U X - U rl + 
iUn, the quantities U r2 (x 2 ,y2). U i2 (x 2 ,y 2 ) are obtained for each point P 2 in 
plane 2. The integral is evaluated over all surface 1 for each point P 2 in 
plane 2. It is claimed the integral cannot be performed analytically and 
integration was performed element by element, by computer. 


U 2 (x 2 ,y 2 ) 


- 

2A J J 
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Figure 4. Method of solving the Kirchoff Integral. 
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If the values for each point in plane 2 are stored, a reiteration pro- 
cess can then find the pattern for plane 3 etc., as was done by Fox and Li. 5 

So far, amplification has not been introduced. 

6.1 Criterion of Step Length 

A critical factor in the calculation is the step length dx 1( d yi and the 
number of steps n-2a/d X1 in one dimension. The number of steps required 
depends on the Fresnel number 

a 2 (7) 

N F “ ~ 

LA 

where (Fig. 4) L is Z 2 -Zi, and A is the wavelength. The Fresnel number is 
the number of diffraction fringes in a width 2a a distance L away. To get 
accuracy at plane 2 there should be about 10 points per fringe, and therefore 

a criterion for n is 

n ^ 10 N p (8) 

Initial runs were made with 20 x 20 steps and therefore should be accurate 

for Np ^ 2. Higher values of Np caused instabilities in the program 
evidenced by very large fluctuations of U r2 ,U i2 and indicating values of 

power in plane 2 greater than plane 1, which is inadmissible. 

For A „ 10 -6 m a value of Np - 2 indicates a laser that is too long and 
too narrow, e.g. if a - 0.1 m then L - 5000 m. More reasonable numbers: 
a = 0.1 m, L-lm give Np - 10 4 and then n > 10 5 . With n 2 increments in 
plane 1 feeding into n 2 increments in plane 2. the total number of additions 
would be n 4 , impossibly large, and the programs too long. 

Two methods were used to reduce computer time, and were tried with 
20 and 200 x 200 runs. The first consisted in storing values only for the 
square ABC (Fig. 4), and then using these values to fill in the mirror images 
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of the rest of plane 2, as axisymmetry could be assumed. However this only 
cut down the time to 1/4. A better and faster method could be used for 
cylinders where it was assuming U r ,Up - 0 outside the dotted circle D in 
plane 1. Then values were filled in only for the 10 elements in the radius 
AB of plane 2. The program then filled in all the units of plane 2 by giving 
them values equal to those units in AB which were the same radius from the 

center B. The number of additions was now n 3 /2, and the program lengths were 
more reasonable. 

6.2 Initial Results with 20 x 20 Arrays 

The first task was to see if runs with only 20 x 20 in 200 x 200 steps 
gave correct results and to recover the results of Fox and Li. Initial runs 
confirmed that a planar amplitude distribution in plane 1 produced Gaussian- 
like distributions with subsequent passes as expected. 3 Results were 
obtained for TM00 modes where the initial amplitude was constant over plane 1 
(Fig. 5a- 1) or was finite only inside circle D of Fig. 4 (Fig. 5a-3). 

Results for a TM01 mode were obtained when the sign of U was changed halfway 
across plane D1 (Fig. 5b). Only the results of the first pass are included 
in Fig. 5 although up to five passes were recorded. For small values of Np 
the program runs correctly. 

6 . 3 Collection Efficiency 

The plots of intensity vs position can give a measure of collecting 
efficiency or the fraction of energy emitted at plane 1 which is received by 
plane 2. The energy emitted at plane 1 is E 1 - r i i (x i ,y i )dx i ,dy i’ ° ver 
the area Ap, calculated by summation over the elements dx^y^^ and that 

received at plane 2 is likewise E 2 - 1^ ,y 2 )dx 2> dy 2 . If gain G 
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for different cross-sections 


Figure 5(a). Results after one pass 

(1) Initial intensity 

(2) Intensity after one pass 

(3) Initial intensity, circular cross-section 

(4) After one pass 
No amplification. N 


= 0.465. Mode is TM00 . 



Figure 5(b). Results after one pass for a TM01 Mode. 

(1) Intensity was assumed constant over initial plane 

(2) U r was assumed to have opposite phase over 
right side of plane 

(3) Intensity after one pass 



occurs in the medium, then a measure of the "collecting efficiency" over an 
area A 2 would be E 2 /GE^ . This quantity will be used in section 6.4. 

6.4 Spherical Mirrors and Focussing 

Analytical methods are possible for plotting the mode patterns of 
confocal unstable resonators using Hemite- Gaussian functions 11 , and computer 
runs gave good results. However the inclusion of gain, and gain varying with 
radius resulted in functions which were not tractable (incomplete Bessel 
functions), and this analytic method was not pursued further. 

The effect of spherical mirrors or lenses can be included into the 
computer program by assuming that the rays emanating from plane 1 (Fig. 4) do 
not originate in the (x,y) plane of plane 1 (Huygen's principle) but from the 
surface of a sphere or paraboloid resting above and tangentive to it 
(Fig. 6). If 0 is the origin, then at (x 1>y i), distance r from the origin, 

(r Jx* + y^) the ray from B is assumed to start not at zj. - 0, but z F - Az. 

Assume a plane wave approaches along front BD. Then for the wave to focus at 
C, the path lengths BC and DO + 0C must be equal. Now OC - f the focal 
length and therefore 

z - Az - f! 

1 4f 

and equation (9) applies to paraboloidal mirrors. 

The program ran successfully for N F < 4 (few fringes mean large spread) 
for various values of f, and focussing was evident; the effect creased as N F 
increased from 0.5 to 4. The method may be applied later for higher N F , 
where the rays would be more beam-like (ray optics). 

6 . 5 Variation of Gain with Radius 

An estimate of the gain c as a function of radius based on laser 
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Method of treating a lens or curved 
surface of wavefront. Rays from D 
at radius r are assumed to start 
from B, where BE = Az. 



To simulate such variations we 


parameters has not yet been carried out. 
assume a is a parabolic function of r. 

The total gain of the laser over a distance L, is G- exp(aL) , where a is 

the gain per unit length. This is the gain in intensity I, which is 
proportional to uj + uf . We assume that to attain G in intensity that U r 
and Ui were both amplified equally by yfZ We now define a gain a' per unit 
length, which is such that exp (a'L) - >/T and use o' to operate on U r . U t . 
It follows o' - o/2, and variations in o and o' with radius are connected 
linearly. We assume that o' (r) can be a symmetrical parabola about the 

axis: 

a'(r) - (1 + 

where p is a constant. The sign of p determines whether the gain is a maximum 
( £ < 0) or minimum ( p > 0) on axis. A maximum on axis might result if the 
effect of focussing the pumping radiation on axis outweighed other consider 
ations. Here we shall consider o' has a minimum value on axis consistent 

with the arguments of section 5.2 and 5.4. 

Consider a maximum gain G max at the edge of a cylinder and a minimum 

gain G m i n at r - 0. The corresponding values of a' are 


a' - In 
max 



a' . - In 

mm 



( 11 ) 


By equation 10 also, 
and can be evaluated 


a max “ S min (1 + ^max 5 


P 



( 13 ) 


It Is convenient to use G max and G min when inserting numerical values. 

Next consider a cylindrical laser with gain a function of radius 
(Fig. 7) and axis at (x^.y^) - (0,0). We calculate the total gain 

G - J a(r) • dl from point (x^ y r to (x^ z,,) . For any point 


at height z the x,y coordinates are 


X£ - x^ , 


Z + X, 


(14) 


y 2 * y l 


z + y 


(15) 


At height z, between z and z + dz the gain a’ is 
a'(z) - a' ( 1 + p (x 2 + y 2 ) ) 


(16) 


and amplitude U r (or U^) is increased by dU r in traversing path dl (Fig. 7) 
dU r - U r a' (z) dl (17) 

But dl/dz - constant and therefore 

dU r - dl U r a'(z) dz 
dz 

Expressing x and y in (Eq. 16) in terms of z, (Eq. 18) can be 
integrated, and leads to an effective a' over the path AB: 


(18) 


“eff " °o V (X 2 ’ X l )2 + (y 2 * y l )2 + (z 2 ' z l )2 


£ 

3 


1 + - (Xj + X 2 X 1 


+ x 2 + 


y?y 


i + y ?> 


(19) 


The total gain g AB between A and B for U r , Uj is then 
g AB - ex p a ; ff R 
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where R Is the distance between A and B. Thus when performing the 
integration of the Kirchoff integral over plane Z L for any point in plane Z 2 , 
for each element at (xi^) at A, its contribution to U r2 , U i2 at point B is 
multiplied by g^g which in turn depends on the coordinates of A and B. 

The function in (Eq. 19) can be checked in the following ways: first, 

it holds for all Z 2 , and if Z 2 - 0 or B is in plane A but not coincident with 

A, we still get gain, as the radiation still traverses a distance AB. 

Second, if x 2 - X]_, y 2 - yj on the beam is parallel to the axis the gain 

reduces to exp (1 - /9r 2 )Z 2 where r is constant. Third, if 0 - 0 and 

a is independent of r, then the gain is exp a'Zo, So the case a' - 

o L 

constant independent of r is included. 

The intensity patterns for a laser where the gain a was 20 at the 
outside and center (Fig. 8-2) show considerably more intensity after amplifi- 
cation than the case for a gain of 20 at the outside and 2 at the center, as 
expected (Fig. 8-4). The question arises that for a parabolic variation of 
a which of the gains G max or G min decides the collection efficiency (Sec. 
6.3). A series of runs were made where G max - 20 and G min varied from 20 to 
2, and the collection efficiency plotted versus G min /G max (Fig. 9). Qualita- 
tively, it can be seen that the power output is approximately proportional to 
^min> although G max was kept constant at 20. 

The results so far have been for plane waves, of constant amplitude and 
phase over the initial wave front. Results were also obtained for waves with 
Gaussian amplitude profiles for similar Fresnel numbers. Gaussian profiles 
are more applicable to an actual case where previous laser/amplifiers would 
have already formed such a profile, with constant phase, and it was assumed 
U r - Uj. over plane 1. The expression used was 

U r - exp ( - (r -r Q ) 2 /Ar 2 ) ( 20 ) 
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( 2 ) 


fife 


Results with a a parabolic function of radius with total 

gain G = 20 at edge of cylinder and total gain G on axi 
max 

(1) Initial intensity (2) Intensity after one pass with 

G • = G = 20 (3) With G min = 10 (4) WiCh G min = 2 ‘ 

min max min 

As G is reduced, intensity on axis falls, but at the edge 

min 

stays almost constant. In (4), the intensity rises at the 
s . N = 0 .465 . 


where r 0 was the coordinate of the axis, and Ar was the radius where U r drops 
to e In our case, for plane Z x enclosing 20 x 20 linear units, Ar was 
chosen as 5 units, and the emitting and collection surface on planes 1 and 2, 
respectively were confined to circles of radii 8 units. Using circles rather 
than the square surfaces did little to alter the results for Gaussian 
profiles when Ar was much less than r. 

The effect of varying G min /G max on the collecting efficiency was 
repeated and the results agreed well with Fig. 9, except that their values 
were about 10% lower due to the smaller collecting area on plane 2. Qualita- 
tively, exactly the same dependence on G min /G max was obtained. 

The variation in total intensity versus G min can be explained by the 
spreading of the beam of radiation. In Fig. 8-2, the intensity at the edge 
of plane 2 is only a small fraction of that on axis; in Fig. 8-4, the inten- 
sity at the edge is as great as on axis, showing the beam was spreading and 
the mode structure modified. These results were for a laser with Np - 0.465, 
and radius of 0.8 mm, length 1 m, with wavelength 1.37 pm. 

VII. CONCLUSIONS 

The emission of large power from a laser would best be accomplished by a 
master oscillator-power amplifier system where the amplifier has a large 
cross-section. For one pass in such an amplifier, the method of solving the 
Kirchoff integral gave results which worked successfully, but only so far for 
small Fresnel numbers - i.e. amplifiers with a diameter to length ratio much 
smaller than envisaged for high power. The preliminary results show that for 
a cylindrical amplifying medium pumped from the outside, difficulties would 
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arise in attaining an efficient absorption efficiency and at the same time 
collecting all the power at the far end of the medium, due to deterioration 
of the mode pattern. 


VIII. FUTURE WORK 

Future work should include: 

1. Extending the above treatment for one pass by increasing the cross- 
section of plane 2. 

2. Extending the above treatment to include onward transmission of the 
radiation after leaving the laser to a plane many laser lengths 
away to find the far field pattern and collection efficiency. 

3. Estimating the gain coefficient for a cylindrical medium pumped 
from the outside which includes the laser paramenters. 

Including the effect of saturation of a. 

5. Converting the programs to include higher numbers of steps, thus to 
treat higher Fresnel numbers, and lasers with larger diameter to 
length ratios. 
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— The possibility of using solar-pumped lasers as solar energy converters is exam- 
ined. The absorbing media considered are halogens or halogen compounds, which are dis- 
sociated to yield excited atoms, which then hand over energy to a molecular lasing medium 
Estimates of the temperature effects for a Biy-COj-He system with He as the cooling gas are 
?**?• " lgh temf * ratures , can cause the lower energy levels of the CO, laser transition to be 
filled. The inverted populations are calculated and lasing should be possible. However the 
efficiency is less than I0" 1 * * * * * . Examination of other halogen-molecular lasant combinations 
(where the rate coefficients are known) indicate efficiencies in all cases of less than 5 x to -3 . 


I. INTRODUCTION 

One concept for collecting solar energy is to use large solar collectors on orbiting 

space stations, which then transmit the energy as laser beams. If the solar energy 

could be converted directly into laser radiation, the inefficiencies in converting the 
energy through different transducers could be avoided. Gas lasers would be advanta- 
geous because they provide a uniform medium, structural problems are reduced, and 
the volume could be large. Here the concept of solar-pumped gas lasers used as 
energy converters is examined. 

The difficulties lie in the fact that the major fraction of the solar radiation lies in 
the longwave (visible) region of the spectrum with a peak at about 2 eV. An efficient 
absorber [1] must be broadband, [2] must absorb near the peak of the solar spectrum, 
and [3] the excited state must be at an energy level sufficiently great to yield reason- 
able quantum efficiency. The halogens can be dissociated at these energies, resulting 
in one atom being in an electronically excited state due to an electron spin flip. The 
excited levels are at a few tenths of an eV, comparable with the levels of molecular 
lasants. Only the halogens and halogen compounds are considered here as absorbers, 
although of course other materials may bq_ feasible. 

Two possibilities are open once the energy has been absorbed: either the absorb- 
ing medium lases or it hands over the energy to a different medium. The case where 
the absorber and lasing medium are the same will be treated elsewhere. In view of 
the high efficiencies of CO* lasers, a BRg-COg-He laser is considered as a first 
example with the helium acting as a cooling medium. 

Recently Gordiets, Gudzenko and Panchenko (1) published a theoretical analysis 
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Fig I Potential curves of the Br. molecule and the absorption coefficient The photons 
corresponding to the shaded area cause dissociation. Those that cause the electron 
raised to the 3 n„- level below AB create vibrational states and hence increase the tempera- 
“ ture. The potential curves were taken from Ref. 3. 


of such a solar-pumped gas laser. The physical mechanisms assumed here are largely 
the same, but our analysis differs in several respects. First, extra nonlinear terms are 
included to take account of the depletion of the particles, otherwise at low pressures 
the above theory predicts more excited states than there were neutral particles in the 
first place. Second, the energy flow is discussed more fully, and estimates are shown 
of the temperatures of the medium at high concentrations of solar radiance. T e 
temperature decides the density of the lower laser level. Third, solutions are ob- 
tained for the inversion density in terms of a wide range of parameters, including the 
gas composition and pressures. Fourth, an estimate of the efficiency is given. 

The physical mechanisms are discussed in Sec. 2, and estimates of densities of 
Br* and the upper laser level are obtained in Sec. 2.1. Section 2.2 considers the 
temperature rise of the laser medium and the population of the lower laser level. The 
inversion population is then displayed graphically as a function of gas composition. 
The efficiency estimates are presented in Sec. 3. Other combinations of halogen 
absorber-molecular lasants are then discussed in Sec. 4 and their overall efficiencies 

estimated. 


II. PHYSICAL MECHANISMS OF A Biy-CC^-He LASER 

The frequency of dissociations per bromine molecule is D = <t>AXcr, where d> is the 
solar flux (photons ernes' 1 per A), AX the absorption bandwidth, and <x the cross 
section at the photodissociation peak. The cross section vanes with wavelength X 
and also with the gas temperature T. The absorptivity of Br 2 vs X and T has been 
studied by Passchier, Christian, and Gregory (2), who showed good agreement be- 
tween theory and experiment. The theory showed that increasing the temperature 
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Fig. 2. Energy flow diagram of a Br^-COr-He solar-pumped laser and for an absorber M. a lasing 

medium Y . and a cooling gas Z. 


from 300 to 800 K reduced the peak value of a at 4200 A by 26 %, but the IK increased 
correspondingly and hence the temperature had little effect on D. 

The absorbed photons cause an electronic transition to a higher upper energy level 
(3) (Fig. 1). and in Br 2 there are two potential curves associated with such transitions 
in the visible. The transitions obey the Franck-Condon principle, and only those 
electrons which arrive in the upper level above the line AB can cause dissociation, 
Br, + hv — » Br* + Br, corresponding to only the shaded part of the absorptivity 
curve (3). Below AB , the molecule ends in an electronically excited vibrational state. 
The transitions into the 3 n, u state all result in dissociation, but the cross section is 
probably lower. It will be assumed only a fraction F of the absorbed photons creates 
Br*. while the remainder causes heating of the gas, with F somewhere between 0.5 
and 0.9. 

It is convenient to summarize the reactions in an energy flow diagram (see Fig. 2). 
The Br* can hand over energy to C0 2 to raise it to an asymmetric oscillation mode 
(001) at about 0.3 eV: CO* + Br* — *■ C0 2 (001) + Br + A£„ with A£, — 0.2 eV, and a 
rate coefficient k x ** 6 x 10 -12 cm 3 s~‘ (4). The Br* can also be deactivated by 
collisions with C0 2 which end up in levels other than (001) [rate coefficient k t = 10~ u 
cm 3 s"‘ (5)], and by collisions with Br 2 [/fc 3 * 4.7 x 10 -13 cm 3 s -1 (6)]. The rate 
coefficient for deactivation by He is negligible. The C0 2 (001) can be deactivated by 
collisions with C0 2 [£ 4 ** 10 -M cm 2 s -1 (7)], He [yfc 5 * 10~ u cm 3 s -1 (7), Br and Br 2 [k t 
= 8 x 10~ 15 cm 3 s~‘ (8)]; however, it will be seen that Br << Br 2 under usual 
conditions. The C0 2 (001) can then lase into the lower level, C0 2 (100). 


II. I . Densities of Br * and the Upper Laser Level 


The rate equation for Br* is 


dBr* 

dt 


= CF<t>AA<r(Br 2 - Br*) - 


Br*[(* t + k t ) CO, + & 3 (Br 2 - Br*)), 


( 1 ) 
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u crki^r radiation is concentrated. If a focusing 
where C is the number of times ^ ^ multiple passes, then C = C 'x. Equa- 

mirror concentrates C times an the expression Br 2 - Br* is written 

tion (1) differs from Ref. I in that, ins e B *. p depletion of Br, is small for 

to include the effect of depleting the Br 2 . The ettect ot ucp 

C< 100 and pressures <1 Ton. 

The rate equation for the C0 2 (001) - * ooi IS 


= t.BfiCO, - s m ) - .v„, lucoo - ' 

d/ 


( 2 ) 


The first parentheses inclu, de the elfe* of « ^ 

differences to the results when the * ^ bm nQt inc|uded in the final expression, 

included: the latter was men ion ^ resuUs for the inversion popula- 

Its inclusion with *« = 8 * lu cm 
tion by about an order of magnitude. 

Using the steady-state solution of Br from Eq. 1. 


Br" 


•Voo, = (kt + k^fa + k^aUk I + Br*/(3.5 x 10- pa) 


(3) 


where the gas composition is descry the ratios « - C Ogr and . = He/Br, 
p is the pressure of Br 2 in Ton; the term 3.5 x 10 pa 
ensures Nom<C 0 2 . 


II. 2. Temperature of the Medium _ 

The lower lasing level is assumed to . be ; 
above ground level. If the temperature of the gases ts T. then oy 

N m = N^exp (-£,«/*!) - «CO, - N«)«p(- » • 

, t x ; n k. The depletion of C0 2 if SI W i becomes 

:; p rec,a;,e B «lTr;^gh CO Xnce is included. Clearly r must be hep, as low as 

^ Maftrst consideration ^ 

average wall transparency ts r. t en tQ thc wa jj come s from the gases. 

Z&Z M then 

perfectly conducted to a ^^“by Stefan's law: 

imately the same temperature <T„) as the wah. tnen r 

(5) 


1/4 

T = 400 j (in degrees Kelvin) . 


v / n "» if n is small. The area A r could be made equal 

30 ^^^^^ ^tiwever, Th^hea^has^t^^^ccmduct^ 

eff r ^s^r^^r^erin, the energy 
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deposited and then assuming it is earned to the walls by gaseous conduction. The 
amount of absorbed energy that goes into heating is obtained from Fig. 2 by consider- 
ing the branching ratios. The fraction of the absorbed Quanta iof average energy 6) 
that ends up in producing C0 2 (001) is Fk t a/[{k l + k 2 )a + * 3 ], and hence the fraction 
contributing to heating (F‘) is 

c< — | _ Fk x a 

iky k 2 )a -+■ k 3 ' ^ 

For example, if F = 0.5. then F' lies between 1 (a-»0) and 0.81 (a->*). 

The heat absorbed per sec per cm 3 is dQ/df. 

^ = CQlkaF'i [Br 2 ] = 10 ~KF'p (Worn' 3 ) (7) 

as <t>.lAcr = 1.15 x 10 -2 (s' 1 ) and i = 2.7 eV for Br 2 . 

The coefficient of heat conduction of a single gas is given by k = kcnC v l 3, where A 
is the mean free path, c the particle velocity, n the gas density, and C v the specific 
heat of the gas at constant_volume. The conduction is independent of the pressure 
(An = constant) and V T. The conduction mechanism is a transport of the hot 
particles with a diffusion coefficient Xc/3; hence, for a mixture of gases 1, 2, 3, where 
1 = He, 2 = C0 2 , 3 = Br 2 , the resultant coefficient k is 


K 


X I 


1 + a 


m 2 C t ., 2 
m x C,, t 


+ b 


r ™ r c r . 3 
V m, c t ,, 


1 + a 


Zl 
cr , 


+ 


b 


Zl 

<?i 



( 8 ) 


where m n = molecular weight, C v , n = specific heat at constant volume, <x„ = collision 
cross section, and 5 is defined as shown. The coefficient is in terms of as He has 
the highest conductivity and is regarded as the cooling gas. 

Assuming a one-dimensional heat flow out from the center of the gas to the 
enclosing side wall: dQ/dt = #cA, ( T - T w /)l, where A, = unit cross section, / = 
distance, then a rough estimate of T at the center is obtained: 


TOO 


5.85 x 10- 


, I C *pl \ 

( f— ) + 7 3' 2 

V K,A t ) 

VS/ 1 * J 


2/3 


(9) 


The numerical constant assumes k, = 3.27 x 10~ 4 cal s -1 cm -1 K _l (9). Apart from 
the one-dimensional approximation, Eq. 9 may also not be accurate if the radiation 
produces considerable dissociation. However, it can be seen that T is a function of 
(xpl/A t ), or a similarity law is obeyed (note that T w contains C'). Increasing p is 
equivalent to increasing either x or l/A x . The gas ratios are described by ( F'/S ). 

Computer plots of T at pressures of 1 and 3.6 Torr on an a-b plane are shown in 
Fig. 3 for C = 100, F = 0.5, T w - 300 K. A flat laser is assumed of depth d = 2 cm, so 
/ =* d/2 and A/l=* 1 for unit area. Then Fig. 3a corresponds to p( Br 2 ) = 0. 17 Torr, and 
Fig. 3b corresponds to 3.7 Torr. The small increase in pressure causes a dramatic 
increase in 7 because more radiation is trapped, yet the conduction does not increase 
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Fig. 3. Plots of gas temperature 77 K) on an a-b plane: a - CO,/Br„ b= _He/Br„ C - 100. F ’0.5. 
7" = 300 K; if ,4// = 6, then (a) corresponds to p = I Torr, (b) to p - -0 Torr. if 4 

corresponds to 0. 17 Torr, (b) to p = 3.3 Torr. 


accordingly. If A/l = 6 (for example. A/l = 2ir for a cylinder, if / 
one-dimensional argument is then inaccurate), then Fig. 3a wou 


radius, but the 
correspond to 


l Torr and Fig. 3b to p - 20 Torr. 

For any given p, the variations of T with a and b are mostly due to the quantity 
1/5, or variations in the heat conductivity of the mixture, which can vary over a 
range of 6 to l, whereas F' varies only from 0.81 to 1. Evidently overheating may 
curtail high-pressure operation of the laser system. 


11.3. Inversion Population 

Evaluation of T yields N 100 by Eq. 4, and the inverted population N = N«ot ~ 
can be estimated from Eqs. 3 and 4. Computer plots of Non a-b planes for different 
pressures are shown in Fig. 4. The ratios a and b vary over the range 10' to 10« and 
N is on a linear scale. Where N is negative, the plots indicate zero. The values 

assumed are C — 100, F — 0.5, A/l = 1. T w — 300 K. 

When the effect of Br 2 deexcitating the CO, is omitted (the k A Br, term), the values 
of N are an order of magnitude greater. Also, failure to include the depletion terms 
results in ATs too high for low p, and low a: for example, if p - 0.001 Tw, and a 
lO" 3 , b = 10- J , the estimated N Ml - 10 13 , whereas there are only 3.5 x 10 molecules 
of CO, present. The depletion terms result in the flattening of the graphs at low p and 

l0 *The drop hi Na t Upvalues of a is the effect of the lower level fdling up because 
of temperature rise for increased CO, content (Fig. 3b). The temperature rise gets 
rapidly worse as p increases and results in the lower level exceeding the upper level 
at about 1 Torr for A/l = 1 . On recalculating N with A/l = 6, the curves for p - 10 , 
10 -2 and 10 -1 Torr are indistinguishable from the above, consistent with the ow 
temperature rise at lower p (Fig. 3). However, the values in the 1-Torr plot are about 
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five times higher, and at 3 Torr there is still an inverted population. All these results 

are ob tain ed assuming the wall temperature is 300 K. 

At T w - 400 K, the curves for 0.001 , 0.01 , and 0. 1 Torr become identical to those 
in Fig. 4, but the 1-Torr curves show maximum N falling from 3.5 x 10 1 * to 3.5 x 10“ 
cm -3 . At T w = 700 K, the 0.001 and 0.01 Torr curves are still unchanged, but the 
0. 1-Torr curve now drops to half the above values, and the 1-Torr curves show no 
inverted population. It seems that if the wall temperature is not kept low, then lower 
gas pressures will be required. 

Reduction of N at high values of b for all p is due to the He deexcitmg the N ooi 
level, yet the heat conduction does not increase with the pressure. 

For both A/l = 1 and 6, the highest N values are obtained for T w = 300 to 400 K, 
C = 100 at p - 0.1 Torr with a - 10- and 10- < b < 10'*, and are about 10“ cm* . 
An inversion population of this value would result in very high gain, maldng practical 
solar-pumped systems superradiant. Hence, if necessary, p, a, and b could be 
changed from their values at optimum N to improve the efficiency. 


C-8 


^ W. L. Harries and J. W W ils*>n 

There is a region of parameters where laser action is possible despite high gas and 
wall temperatures. Plots of .V at b = l and 10- show that at b - l. N is drastic y 
reduced Ts I, is increased, but when b = 10 - there is still an inverted l population 
even at T u = 1000 K. (However, this result should be treated with caution as tem- 
peratures this high can introduce rotational effects which are i not inc J h J gh 

lemnerature region occurs where the concentration of He and CO, is small. At nign 
erpe a u es H as a coolant is unnecessary: Eq. 8 converts to the resultant con- 

: ::: z 0 ^ ga se S - o. and wu h v* ** 

adequate. Again, as ^ 0 the deexcitation of CO, (001) by He becomes unimportant. 


III. EFFICIENCY AND POWER OUTPUT OF THE LASER 

///./. Efficiency 

The overall efficiency can be subdivided into the product of several 
First, if absorption occurs over a bandwidth X, to x„ the (iractton of the solar spec 
trum used or “ solar utilization efficiency" is 17 ,: 


■q, = p<t>(x)dX/ j X <t>(x)dX = 0.18 


( 10 ) 


'"if the radiation passes through the gases x times by reHection andthe^abso^^ 
small, the fraction absorbed or -absorpnon efficiency b, <B ^ Iem £ ratur( , 
of C' and obeys a similarity law, r) A = r) A (xpd). However, p is am y 
rise considerations and x is limited by the reflectivity of the reflectmg wall (-99 %) to 

SSHSsS- 

“kinetic efficiency,’ *)*• Assuming F 0.5, , ’ d © 19 as a — ». 

10-3 cm a s - l; then „ - 0 as * - 0, r, k = 0.18 for u - 1, and the 

Jaw - 

fra^n' tttS&rtJ* 

(or slightly more, depending on F) is PK»bte. It is ^ ^ is sufficien t 
efficiencies are achieved only in those regions of Fig. 4 where me ga 

for steady lasing. 


Ill 2. Power Output of the Laser 
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a 


Fig. 5. Effictency of a Br^COj-He laser at *Br f ) = 0. 1 Torr vs a = C<VBr 2 for different numbers of 

passes x. 


power. Assume A, faces a collector of area A c so A c /A, = C \ If the solar flux is 0. 14 
W cm- 2 , the power absorbed per cm 3 is 0. and with tj. = 0. 18, n A — 1 , then 

the power absorbed P* = 2.6 x 10" 2 C' (Wcm" 3 ). The laser power emitted is 
Pab'nicVQf and at a ^ 0. 1, tj* = 0.18, t} q = 0.039 the power emitted is 1.8 x 10 _4 C' 
(Wcm 3 ) = 18 kWm -3 for C' = 100. A laser 1 m 2 and 2 -cm thick should emit up to 
360 W and would require a coUector of 200 m 2 . A 10-kW laser would have a volume 
of 0.56 m 3 and would require a collector of 5.5 x 10 3 m 2 . Such a system is at least 
comparable with other possible solar-pumped laser systems (12). 


IV. OTHER SOLAR LASER MATERIALS 

As the efficiency of a Br^-CO*— He laser seemed low, other combinations of mate- 
rials were considered. The difficulty of obtaining rate coefficients especially for 
collisions between unlike components led to estimates of the laser performance 
based on a series of “figures of merit.” The method assumes the partial pressures of 
the gases are equal — an oversimplification. Even so, it turns out that the results are 
useful for indicating poor performance, which usually is the case where one rate 
coefficient far exceeds the others, and the others can be neglected. The figures of 
merit can only be used in comparing one system with another. 
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The energy flow chart for a laser consisting of molecular absorber M , a molecular 
laslm id a coolant Z is included in Fig. 2. Material M is a compound containing 
he halogen / and dissociates to yield the excited specie X*. The excitation energy 
oT£t‘hen hLed over .o .he upper level r v of .he lasing med.um r. The lower 
lasing level is Y L . The rate coefficients are defined as shown. 

The absorption coefficients « (mole- /cm-) for the halogens and halogen com- 
Dounds (10) are plotted vs wavelength and compared with the so ar ra 1 

nr 1 A" 1 ) < 1 1) in Fig. 6. It can be seen that the three compounds of the perfluoral- 
kW halide aZ have high absorption in .he ultraviolet and >ha. -vera, ma.enab 

cL absorb up to 6000 A. The absorbing power of a smglej M 
a (cm J ) = 1.66 x 10-* 1 e and D is averaged over wavelength X. The o(\) can be 
exoressed by a single or the sum of several Gaussian curves vs X. and AX is the fu 

width of the curve at the rms value of <r. Values of D are » vcni " t ^ sta ^^ and 
Table 1. Only a fraction F of the absorbed photons will yield the excited state X*, and 

F ^ rm h etr S of the density of excited halogen atoms^ 

“ > i, Sen “ is essentially determined by D and and a figure of men. for 
absorber storage is defined: 


F, = Dfk, . 


(ID 


A low value of F, indicates that the absorbing material M is deexcitmg X* by colli- 
sions, and a high value of D does not necessarily guarantee lasing. 
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TABLE 1 

characteristics of the halogens and halogen compounds \ = wavelength fob 

THE PEAK ABSORPTION. Ax = BANDWIDTH ,t = ABSORPTION CROSS SRTinin*^ 
WHERE . - SOLAR RADIANCE. - DEEXCITaTION RATE COEFFICIENT. 51"™ 6 Fret RE Of 

M C D IT r - n i cunr-l rx fa r' a o , . r IUJ L, I\C KJ r 


Gas 

A ,, 

(Ai 

Aa 

(A) 

<7„ 

(cm*) 

D 


Ti 

(‘for Br*) 

Dk, 

I. 

5000 

640 

l.70x I0 " 

8.4 X 10--’ 


3.6x 10" 

2.33 < 10* 

Br. 

4889 

378 

1 50x io"* 

4.33 X 10- ' 





4141 

360 

3.9x |0- |a 

7. 18 x 10-' 

1.15x10-’ 

4.7 x |0 " 

2.45 x 10'" 

Cl, 

3283 

332 

1 . 19x I0' 1 * 

1 .23 x io- 1 




(Cl 

4804 

284 

8.42x 10- io 

1.80x10-' , 





4425 

549 

1 .66 x 10“ 19 

6.03 x )0- J 





3249 

254 

1 .56 x 10--° 

1.18x10-' 

• 7.2x10"' 

— 

— 


2440 

258 

2.29 x i0- |a 

4.58x I0~ 4 ^ 




IBr 

5072 

291 

2.82x |0-'* 

6.55 x 10“ 3 | 





4773 

508 

4.78 x 10"'* 

1.79x10-- 

• 2.5x|0- 2 

* 1 .0 x 10' 1 - 

2.5 x 10'° 


2682 

277 

1.31 x 10-'* 

4.74x I0" 4 J 




BrCI 

4555 

404 

4. 17x IO-™ 

I . I7x 10- * 





3736 

354 

2.01 x 10-'* 

3 29x |0- 3 

4.5 x 10" 3 

*2.9x IO” 14 

1 .55 x 10" 


2278 

261 

3.54x 10-*’° 

4.7! x 10 5 




CFJ 

2650 

210 

5.40x 10-'* 

1.37x IO' 3 


6x 10“ 57 

2.28x 10" 

C,F : I 

2745 

270 

6.20x |0-'* 

2.45 x 10- 3 


1 x IO" 17 

2.45 x 10" 

C\F*I 

2880 

280 

5.60x 10-'* 

2.88x IO" 3 


— 

— 


The “transfer efficiency” is a measure of the fraction of collisions between A”* and 
y that result in the formation of the upper level Y t . Actually the fraction should be 
k x Y(k 3 M + k z Y), but here it is simplified to the number F T , defined by 

F t = kj(k t + k 3 ) s kjk 3 . (12) 

The expression is useful if k 3 M >> k z Y, and also assumes the densities Y and M are 
of the same order. 

The “upper laser level storage'’ is a measure of the density in the upper laser level 
Y c . The rate of production of Y v is proportional to k x X* Y, and the rate of loss is 
proportional to k<M + k s Y + k*Z. The density Y v is large when the quantity F u is 
large, defined by 

F u = ki/(k A + k s + k t ) ~ k x /k< . (13) 

The value of the individual figures of merit lies in that they can pinpoint the quality 
of performance in one particular reaction or region of the energy flow chart. 

Numerical values for D and the rate coefficients for the systems Br^COj-He, 
Br r -H 2 0-He, Br^HCN-He, and I*-HF-He are shown in Table 2, which also gives 
the figures of merit and the quantum efficiencies. 
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Fig. 7. Overall solar efficiency v of various laser combinations as a function of a = M/Y i- las- 

ant/absorber). Complete absorption is assumed. 


Despite the very large cross section of I, for photoabsorption, the high rate of 
deexcitation of I* by I, (*J reduces the F, value to one-tenth that for Br„ implying 
that for equal conditions the value of I* is an order of magnitude less than Br* 

The coefficient k , is much less for I^HF than any of the Br, systems, implying 
poor transfer of energy to HF, ( F T ), and also low retention in the upper level because 
of deexcitation (FJ. The highest value of is that of Brj-HjO, which also has the 
highest transfer efficiency F T . Unfortunately, it has a very high deexcitation rate of 
the tipper laser level (low FJ. The Br r HCN-He system has values of F T and F u 
which are intermediate between the combinations of CO* and H,0. 

The overall efficiency of the laser can be expressed as tj = ai) k r \ Q , where it is 
assumed that complete absorption occurs. The constant a = 1.4 x 10— 8 £>e/<r and for 
Br* »s 1*32 x 10- and for I, is 1.4 x 10- units. The value of F, the fraction of excited 
species formed, is taken as 0.5. Plots of 17 vs a, the ratio M/Y{ Fig. 7), show 77 » 5 x 
10 3 for Br*-H,0-He and Br,-HCN-He, provided a > 10— and for Ir-HF-He if a > 
10*. Unfortunately, in the latter case the high value of a will probably result in the 
gases becoming overheated, thus filling the lower laser level. The BrnCOt-He laser 
has an efficiency of around 10— for a > 0. 1 . 
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V. CONCLUSIONS 


A solar-pumped Br 2 -C(VHe laser should lase provided the radiation is concen- 
trated sufficiently (> 100 times) and the pressure of the gases is low. High pressures 
cause overheating, which results in the lower laser level be.ng filled. A suitable 
choice of gas pressures for optimum inverted population would be Br 2 — 0.1 Torr, 
CO = 10~ 3 Torr He = 10~ 3 Tort: but for greater efficiency a better combination 
would be Br> = 3 Torr; C0 2 , 10' 3 to 10' 2 Torr; He, 10- to 10— Torr (Fig. 4). 

As a power source, the laser is limited by low efficiency. The efficiency is ^the 
product of the fraction of the solar spectrum used, the fraction of radiation absorbe , 
a kinetic efficiency depending on the branching ratios, and a quantum efficiency. The 
absorption efficiency can be improved by increasing ( xpd ). The pressure can be 
increased to about 3 Torr if the cooling is enhanced by introducing fins into the gas 
for good heat conduction to the walls. Alternatively, a number of flat "box lasers 
could be placed on top of each other to absorb different sections of the spectrum m 
sequence, with each producing its own output beam. With complete absorption, a 
10-kW Brz-COz-He laser would require a collector approximately 70 x 70 m, and 
would have a volume of about 0.5 m 3 ; such dimensions would be feasible for space 
applications. The radiator has to be large to reduce the temperature, {A e lA r -* D- “ 
the temperature is high, lasing still seems possible at low concentrations of He and 
C0 2 , but at further reduced efficiency. 

Other systems using halogen absorbers were also examined; namely, Brj-UUr- 
He B r2 -H 2 0-He, Br,-HCN-He, and Ir-HF-He. Assuming complete absorption the 
efficiencies were all below 5 x 10' 3 . A 10-kW laser with this efficiency would require 
a collector of several thousand m 2 and a radiator of comparable size. 
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LIST OF SYMBOLS 


The following represent the densities (cm -3 ) of the respective materials: 
Br>, Br, C0 2 , CO 2 (00l) = fy*,,; (001) is the vibrational state; 

Br*. X* excited atoms 

V/ molecular absorber containing the halogen atom X 

Y lasing medium; K t upper, Y L lower level 

Z cooling medium 

,V inverted population 

Density ratios: a = C0 2 /Br 2 ; b = He/Br 2 
pressure of Br 2 (Torr) 

number of times focusing mirror concentrates 
number of passes through lasing medium 
C'x — effective concentration of solar radiance 
solar radiance (watts per unit area per Angstrom) 
bandwidth of absorption 

cross section for absorption by molecule (cm - *) 

Planck's constant 
Boltzmann's constant 
frequency of photon (s -1 ) 
average energy of quanta 
excess kinetic energy released in absorption 
energy of C0 2 in ( 100) state 
Rate coefficients: 

for Br* handing over energy to CO 2 (001) (cm 3 s -1 ) 
deactivation of Br* by C0 2 which ends up in levels 
other than the upper laser level 
deactivation of Br* to Br by Br 2 
deactivation of CO 2 (001) by C0 2 , He, and Br and Br*, 
respectively 

fraction of absorbed photons resulting in Br* 
fraction of absorbed photons contributing to heating 
areas of collector, heat radiator, and subtended by laser 
unit area of cross section of laser 
depth of laser 


P 

C' 

x 

C 

<t>(A) 

AA 

<T 

h 

k 

V 

t 

A£, 

£100 


^3 

* 4 . *». 

F 

F' 

Af, A,, 

A, 

l 

T)i. 

Va 

T 

t 

Q 

K 

A 

c 
n 

m n 
C r .n 

s 


Vk, 

V 


solar, absorption, kinetic, quantum, and overall efficiencies 

transparency of container wall 

emissivity of heat radiator 

quantity of heat deposited 

effective heat conductivity of mixture of gases 

heat conductivity of nth gas 

mean free path of gas molecule 

mean gas particle velocity 

gas density (cm -3 ) 

molecular weight of nth gas 

specific heat at constant volume of nth gas 

constant representing the effects of three gases on 

heat conductivity (Eq. 8) 

gas temperature 

wall temperature 
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Abstract — The possibility of using an IBr laser as a solar energy converter is examined 
theoretically, and reasons for its choice are given. Broadband absorption results in dissocia- 
tion with the formation of excited Br* atoms, some of which then lase to the ground state Br, 
The ground state is depopulated by three-body recombination and, more importantly, by 
exchange reactions which more than compensate for the high quenching in heteronuclear 
halogen systems. Kinetic modeling indicates lasing is possible in the pulsed mode and 
possibly in the steady state with a cooled gas flow system. Temperature effects are dis- 
cussed. The efficiency of the laser approaches 1.2 9c at optical thicknesses large enough for 
complete absorption of the photons. 


I. INTRODUCTION 

The concept of collecting solar radiation in large mirrors on orbiting space stations, 
and then transmitting the energy via laser beams has been considered previously 
(1-3). The efficiency of the system is expected to be highest if the laser could be 
directly pumped by the solar radiation. 

The criteria for an efficient solar pumped laser are as follows: 

• There must be broadband absorption; 

• Peak absorption should occur near the peak of the solar spectrum; 

• High quantum yield into a long-lived (metastable) state which serves as the upper 
laser level; 

• In general, quenching of the excited state should be small, but, as will be pointed 
out later, this condition is alleviated if 

• The lower level is rapidly depopulated to maintain inversion; 

• The upper and lower levels must be sufficiently separated to yield a reasonable 
quantum efficiency; 

• The process must be reversible; if not the components must be reconstituted by 
flow methods. 

Gas lasers are advantageous because of uniformity of medium and because size is 
not a limitation. Two classes of laser can be considered: [1] where the absorbing me- 
dium is distinct from the lasant, and [2] where one material performs both functions. 
The first solar pumped laser examined theoretically was type [1] above, namely a 
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Br -CO, -He mixture (2-3). The Br 2 acted as a broadband absorber, the C0 2 was the 
lasing medium, and the He acted as a coolant. This laser was an electron, c-to- 
vibrational energy transfer laser, and the low -transfer efficacy . as well as ab- 
sorption efficiency resulted in an overall efficiency of less than OAWc 

Higher overall efficiency might be attained if the transfer efficiency could b 

eliminated, as in lasers of type [2] above. An example would ^ 
a molecule to yield an excited atom, which then lases to the ground state, 
population of the lower level might be removed by chemical processes. High pres- 
sure working should be possible which would enable more efficient absorp ion 
The objectives of this paper are to study theoretically the IBr solar pumped laser 
as In example of class [2] to understand the essential lasing features by determining 
“he dominating reactions, and to estimate the efficiency for power conversion. The 
criteria that have been listed for the selection of candidate lasants will also be 
assessed. In particular the question of whether steady-state lasing is possible will 

^TsThe solar radiance must be concentrated many times, and since the IBr absorp- 
tion near the peak of the solar spectrum, the laser will qutckly heat up 

and lasing inhibited (by mechanisms to be discussed), unless cooling is pro 1 ^ 

the theoretical study it is therefore assumed that the lasant is maintained around roo 
« ^p^ being .o examine ,he po.emial for lasing under favorable 

operating conditions. Methods for accomplishing cooling are proposed in Sec IX. 

WhiTe th“s udy was in progress, an experimental investigation of an IBr la e 
Dumped by a xenon lamp was performed by L. Zapata (4). The objective of that 
experiment was to demonstrate lasing, and not necessarily to provide definitive 
Quantitative data Hence, only limited comparisons could be made, but the expe 
ZZTsu's were useful in assessing the effects of excessive heating ,n the actual 

experimental environment. 


II. CHOICE OF IBr 

Broadband absorption is essential for high solar r effic I^y ' 
compounds ty^i diatom* homp y 

«• - rao,e£Ules 

° f hH^^^s, rated for type «. using a -on 

can absorb near the peak of the solar spectmm 

W m n K b0 ‘ h f ramL a ^ h a*“nd I* have energies of about 0.1. 0.2, 0.44. and 
The excited atoms F , Cl , Br , ana i ^ ^ values high enough to give 

acc^ptabl^quantunf efficiencies. Lasing ^V^popu- 

between the rate of reduction of the excited species X by quen^i^ £ 

r^e'^ 

since*" t^^nmm™ f^ieMylnwea^ 1 wilhlhe atomic number of the lasing atom, it 
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appears that IBr should be the best laser candidate of type [2]. It absorbs near the 
solar peak with a high probability of dissociating into I * Br*. 


III. PHYSICAL MECHANISMS 

The processes occurring in IBr are assumed to be photodissociation with the 
formation of excited and ground state atoms, quenching of the excited atom re- 
combination and exchange reactions. The species are assumed to be eight . all, 
namely, IBr. I. Br, L, Br>, Br*. and I*. and photons created by spontaneous and 
stimulated emission. Vibrational excitations are neglected at room temperature be- 
cause of very rapid relaxation. The reactions and rate constants are summarized in 
Table 1. 

The absorption of photons results in reactions 1 through 6. The photodissociation 
rates 5 are equal to Cd>(A)AAcr„(A)(yV), where C is the number of times the solar 
radiation is concentrated, <t>(A)AA is the number of photons arriving per unit area per 
second (8), cr„(A) the absorption cross-section in question, and hY) the number den- 
sity of absorbers. As the ratio Br*/Br initially produced by the photodissociation of 
IBr is critical to attaining an inverted population, care must be taken in evaluating 5i 
and 5 4 in Table 1 . The potential curves for the A'CI 1 '), A 3 n,, and B 3 n 0 ^ levels for IBr 
were plotted by computer on the energy diagram (Fig. I) using data from Huber and 
Herzberg (9). The B' 3 n + is included, and is drawn in approximately. The lowest 
vibrational level (v = 0) and the Franck-Condon transitions are shown. The horizon- 
tal dashed lines B and C are based on the transitions from the lower to the upper A 3 n, 
and £ 3 n 0 levels and indicate the peaks and widths of the absorption curves. The 
absorption cross-section for IBr as a function of wavelength has been measured 
(10-11), and the function can be represented by three Gaussians whose peaks are 
at 268. 477. and 507 nm, respectively. These are plotted on the left. The peak of 
Gaussian F coincides almost exactly with the line C. but the peak of the Gaussian G 
is displaced from B. We believe the Gaussians to be accurate as well as the asymp- 
totes .4 and D. 


The integral 


f d>(A)cr „(A) dA 

is a measure of what fraction of absorption events are caused by transitions to the 
level in question and <r a (A) is the respective Gaussian. Thus Gaussian H produces Br. 
but can be neglected as </>(A) is small here (8). The part of Gaussian F above 
asymptote A produces Br*; the part below A does not produce dissociation. The 
£' 3 n,7 and BTIj" curves cross and possibly absorption into the latter could result in 
Br. However, the translational energy at the crossing is sufficiently high and the 
nature of the crossing is such that the probability is near unity that all parts of the 
absorption curve above A in Fig. 1 correspond to dissociation into I + Br*, as 
confirmed experimentally (13). The Gaussian G, corresponding to absorption into 
the A 3 n, level produces I + Br. Performing the integration of d>(A)cr a (A) dA then 
yielded the fractions of absorption events resulting in the production of Br*, Br. and 
I (Table 1). The total absorption rate for IBr was found to be C x 1.25 x 10 IS (IBr), 
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TABLE l 

r 1ST OF REACTIONS 



Scarce 


Quenching 


4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


hv * IBr 


hv U 


hv -+* Br 


Recombination 
(2 body) 

Recombination 

(3 body) 


Br* 

Br 


l +■ 

I - 

Br 

Br 


Br* 

Br 


Br 


IBr — Br + IBr 

- l 2 Br - l. 

-r Br, — Br ^ Br 2 
+ I — Br + I 

- Br - 2Br 
l* + IBr — I - IBr 

+ 1 2 —• ► I + L 

- Br, — I - Br, 


15 l, + Br, ^ 2lBr 


Exchange 
<2 body) 


16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


Br* + Br - IBr 

— Br, + IBr 
Br + Br ^ IBr 

— Br, + IBr 
Br* + Br + Br 
Br + Br * Br, 
Br* +■ l + IBr 
Br ^ l ^ IBr - 
I* + I + IBr - 
+ IBr 

I* + Br + IBr 
I + I + IBr — 
+ IBr 

Br + IBr — ► I 
Br + I 2 


5„ = 

Br*: 

Br: 

I: - 

5,i = 

l*. - 
I: - 
S:u : 

Br*: 

Br: 

0. 
Qz 
Q> 
Q* 
Q * 
Q 6 
Qr 
Q« 

K; 


C x i.25 < 10 l5 ( IBrl 

- 5, = 0.7045,, 

- 5 4 = 0.259S,, 

S- = 0.%35 n 

: C * 2.5 x 10 15 (L) 

■ S 2 = 0.25,i 
S* = 0.485,, 

= C x 1.25 x l0 15 (Br : ) 
: - 5, = 0.60 55, , 

- 5r = 1 . 38 1 5n 

= I x 10“ 12 (cm 3 s' 1 ) 

= 1.86 x 10“ 12 
= 4.7 x 10“ 13 
= 1.8 x I0“ u 
= 1.8 x !0“ li or 0 
= 6 x 10“ u 

= 3.5 x ir u 

= 5.6 X ir 11 


- I x 10' 13 
-(1.6 x I0" 16 ) 


c x = 4 X 10" 32 (cm 6 s"‘) 


9. 10. 11, 13 
14. 15. 16 
17. 18, 19 
20 . 21 


22, 23 
22 
22 
14 

24 

24, 25 
24. 26 

31 


Assumed 
equal to C { 
Assumed equal 



c t = 3 X 10- 30 

to C 

- 2 Br 2 

f, = 4 x 1 (T 31 

27 

2Br 2 

c 4 = 3 X 1CT 30 

27 

— 2IBr 

f s = 1 x l(T 3J 

28 

- 2IBr 

II 

X 

o 

k 

28 

- U 

n 

ii 

X 

© 

k 

28 

— 2lBr 

C. = 3 x KT” 

28 

I, 

c t = 3 x ur 30 

28 

+ Br 2 

£, = 3.5 x nr" (cm 3 s~') 

29 

• + I 

£, = £, or 0 

29 


where C is the number of times the solar radiation is concentrated, of which a 
fraction 0.704 resulted in Br*. 0.259 in Br and 0.963 in I. Similar estimates for I 2 and 

Br *Fo7aU l the" 'inititf pressures of IBr considered, it was assumed that 4% by number 
density of l 2 and Br 2 were present at 300 °K, according to the law of acUom 

Thereafter, ’all the densities varied with time. The ^nd, TaCtTo to” 

IBr pressure for low pressures (and absorption lengths of cm), 

essentially all the photons are absorbed. ..... c ,u» 

The degeneracy factor for the upper laser level Br* is 2, while that o 
laser level Br( 2 P 3(2 ) is 4. The population inversion SN for lasing is the 
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R C Rn g s t . ) 

IBr LOWER LEVEL: X1SIG+ 

UPPER LEVELS: B3PI0+ P3PIi 


Fig. 1. Energy level diagram for IBr. The potential curves AVI*). A 3 U { . fl 3 ITv are Morse functions 
generated from data of Huber and Herzberg. The repulsive curve £' 3 n„* is approximate. The Gaussian 
absorption curves are plotted from data of Seery and Britton. 


A.Y = Br* - Br/2 . (1) 

The degeneracy factors in this instance help to reduce the threshold for lasing. 

The quenching of Br* and I* is given by reactions 7-14 in Table 1 ; the quenching 
of Br* and I* by Br* and I* was neglected as well as the quenching of Br* by Br and 
I* by I and Br. Computer runs in which these latter coefficients were arbitrarily 
assigned values equal to Q 4 , instead of zero, showed no great difference in the 
results. The large rate coefficient for the quenching of Br* by I is due to the elec- 
tronic to translational energy transfer resulting from the “crossing’ of (IBr)* poten- 
tial energy curves (see Fig. 1). The process has been called the inverse predissocia- 
tion mechanism (14). 

The two body recombination of I 2 and Br.. (item 15) is the only reversible reaction 
included. At room temperature, by the law of mass action (assuming no photodis- 
sociation), the concentration of I 2 and Br.. in IBr is 0.04. It then follows that if the 
forward reaction 2IBr — ► I, + Br 2 has a rate coefficient K 7 , then the reverse coeffi- 
cient K h is (0.04) 2 x K 7 . 

Three-body recombinations are listed and it is seen that the reactions involving I* 
and Br* are less likely than those involving I and Br because of the difficulty of 
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e£,u,?onS£s?»n^ 


getting rid of the excitation energy. Some of the values are uncertain and are based 

° n The«change reactions that depopulate the lower laser level and J""* 
possible are listed as items 25 and 26. Experimental measurements of react on 25 
tave been repotted by Clyne and Craze, (27) who deduce a high rate coefficient E, 

3 5 x erf s- . We have been unable to find a rate coefficient for reaction 26^ but 
computer rans for E, = 0 and ft = E, did no, differ greatly because the concentratton 

of I 2 is much less than that of IBr. 


IV. RATE EQUATIONS 


In the following equations the terms in 
atomic and molecular species; n is the density of photons 


parenthesis represent the densities of 
The quantity A is the 
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Einstein coefficient for spontaneous emission, while T is the rate of stimulated emis- 
sion. The rate equations for the eight species present are 


— = -5,,(IBr) - C 3 (I)(Br*)< IBr) - QlMBrXlBr) 

+ C(I*)(Br)(IBr) - £,(Br>(IBr) + £,(Br)(I 2 ) 

* 2AVI.)(Br 2 ) - 2£ 8 <IBr) 2 

^ =5 7 + S, 2 S 5 - C 5 (I)(Br*)(IBr) 

- C s (I)(Br)(IBr) - C 7 (I*)(I)<IBr) 

- 2C 9 (I)(I)<IBr) + £,(Br)(IBr) + £ 2 (Br)(I,) 

^-=- 5 21 (L) - C^I*)(I)(IBr) + C 9 (I)(I)(IBr) - £.,<Br)(I 2 ) 


( 2 ) 


( 3 ) 

( 4 ) 


52l<Br - ) + c i<Br*)<Br)(IBr) + C 2 (Br)(Br)(IBr) 
i- Cj(Br*)(Br)(Br 2 ) + C 4 (Br)(Br)(Br.) 

* £,( Br)( IBr) - KAMBr,) + £ 8 (IBr) 2 (5) 

rlf 

° =5, + S 3 - C,(Br*)(Br)(IBr) - C :i (Br*)(Br)(Br 2 ) 

- C 5 (Br*)(I)(IBr) - (?,(Br*)(IBr) - 0 2 (Br*)(I 2 ) 


- £> 3 (Br*)(Br 2 ) - £> 4 (Br*)(I) - Q 5 (Br*)(Br) - A(Br*) - T (6) 

d( d Br) = S 3 + S 4 + 2 S 6 + £>,(Br*)(IBr) + £ 2 (Br*)(I,) + 0 3 <Br*)(Br 2 ) 

+ £ 4 (Br*)(I) + 0 s (Br*)(Br) - C,(Br*)<Br)(IBr) - 2C 2 (Br)(Br)(IBr) 

- C 3 (Br*(Br)(Br 2 ) - 2C 4 (Br)(Br)(Br 2 ) - Cg(Br)(I)(IBr) - C„<Br)(I*)<IBr) 

- £,(Br)(IBr) - £ 4 (Br)(I 2 ) + A(Br*) + T (7) 

d( d I * ) = S 2 ~ Q4 < - (? 7 <I*)(IBr) - £?g(I*)(Br 2 } 

- C^n(IXIBr) - C 8 (I*)(Br)(IBr)' (8) 

^-=V + GA(Br*) <9) 

dr t c 


The quantity G in Eq. 9 is the geometrical factor 2 r\jL}, where r b is the radius of 
the laser beam and L the length of the laser; for isotropic spontaneous emission only 
this fraction of photons is confined within the laser cavity. The last term represents 
the loss through the end mirrors of reflectivites r, and r,, and r c is the lifetime of an 
average stimulated emission photon travelling parallel to the axis; 

t c = -T/clnlr/j), 


D-8 



Iqq W. L. Harries and W. E. Vendor 

f .-u t u*ro r = s X 10* s. The quantity T, the stimulated 
rmisTionrate. Ise^ to na e clS where ex, is the stimulated emission cross-section. 
Then 


(j. 


= \TA/47T 2 AX ( 


where x is the emitted wavelength. A l - 11 is the Einstein coefficient for s P»"' a “' 

frnm t Ua imner laser level and AX. is the emission bandwidth. With AX, 
ous emission from the upper laser ievei. <*»« , _ , x 10 -n cm 2 

taken as the bandwidth for Doppler broadening at 300 K, then <r e - 1 b x tu c 

Fouations "through 9 were solved by computer, and are compared with the 

q t i ^ c,,itQ of Zanata in Figs 2a and b. In the experiment, the light intensity 

was measure/ with a pho.ode.ec.or, and varied approximately 


as 


sin^ 


27 Tt 

2 x if 


where t is in seconds. The theoretical source terms included this function. 


IV. COMPARISON WITH EXPERIMENT 

Ju^t— 

K3 =er * ' St 

ST- reflector was no.Terf^y earned 

radianc^^Thr vLues^The equivalent ^^^Kes^^he 'sp>ctral 

2000 and 40,000 for the vanous runs. There were also a 

distribution of the xenon lamp and the solar spectrum ““Z- 1Br , aser pumped by 
Figure 2 compares the measured light output toma 1-m J„ 

a xenon discharge lamp powered from a capa i ’ j n t h e 

of Eas 2 through 9. The IBr pressure was 3 torr and C was aoo 

experiment, which are the values assumed sim iiar. Both 

The overall shape of the experimental and emulated ^»e d P ^ ^ ^ ^ of 

start about the same time, and consist o a s H> oscillations in the theoretical 

spike to tail is higher in the calculated version, constanl of the 

ZZ pJXdetector chcuhshowed'h was high enough to reduce the high frequency 
‘“ZTm^f^perimenta. puls, for 

shorter than the calculated values by a tempera- 

to the lasant being heated near the axis. Rou * . f ral hun dred "K in 50 (is. 

s«*. a value too high to have been 
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caused by photodissociation (it would require C 2= 10 5 ). High temperatures could 
also increase the ratio of reverse to forward rates of the exchange reactions, thus 
tending to neutralize the mechanisms necessary to depopulate the lower laser level 
(30). Also, the quenching reaction Br* + I — Br - I is enhanced in the heated gas 
because of the greater number of I atoms. The calculations assumed no heating. 

On changing the pressure p both experiment and theory showed the output laser 
pulse duration was roughly independent oip for pumping pulses of 200 /us duration. 
Experimentally the output light signal amplitude increased withp up to about 20 torr, 
where it reached a maximum and then decreased: asp increased, absorption in the 
outer layers decreased the pump power on axis. The theory showed the amplitude 
proportional top because constant pumping energy density was assumed throughout 
the gas. 

Both experiment and theory showed that increasing C increased the laser ampli- 
tude proportionally. 

There was also agreement between the measured and calculated output power of 
the laser. In an experiment with a tube length of l m and 2.22 cm radius, and an 
output mirror of 95 % transmittance, the measured output was approximately 2 kW 
for an IBr pressure of 4 torr, with C = 10 4 (4). 

The calculated peak output power P is 

P = cnhvTnr 2 . * 10) 

Our results showed n = 1.2 x 10 12 cm' 13 for C = I x 10 4 , corresponding to 1.7 kW 
which is in reasonable agreement. 

V. INTERPRETATION 

By varying the parameters in the program one at a time, it was established that the 
dominant mechanisms were quenching of Br* by IBr and I and the exchange reaction 
Br + IBr - I + Br 2 . In particular, the program indicated lasing would not occur if the 
exchange reaction was removed. Three-body recombinations made little difference to 
the shape of the light output for pumping pulses of 200 /lis duration. Fortunately the 
rate coefficients most accurately known turned out to be dominant ones. However, 
for computer runs of longer duration pulses, it was found essential to include all the 
terms. 

Evidently the exchange reactions can overcome quenching, and high quenching 
cross-sections per se may not be a valid selection criterion for rejecting otherwise 
promising gases. 


VI. THRESHOLD CONDITION 

The threshold condition is r,r 2 exp(2a:L) = 1, where r^are the reflectivities of the 
laser mirrors, a = is the gain per unit length, and L is the length of the laser 

(29). The threshold condition is contained in Eq. 9, and occurs when 

T > t)/t c - GA(Br*) =* n/T c . 

The two relations are equivalent as can be seen from the definition of r c and require 
AN • L 2= ln( l/r,r 2 )/2cr e = constant . 
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Id8 


r 



L (cm) 


Fig. 3. Plot of lowest C for lasing vs laser length L. 


As AIV « C it follows .ha. .he Uphold value for C " l^andMcm' 

runs were performed for successively °^ er va J^ s compared with experimental 

and .he threshold values of C estimated^ The values ^re compa M the 

walttll' ZSZ ; =ry and -penmen. show C 
rrf-i Thus C can be reduced by increasing L. 


VII. EFFICIENCY OF THE LASER 

Estimates of .he efficiency can 

For a laser of length L collec.mg w.dth w 0Jt L. If the 

solar power into the laser is 0. ... x 10 - '°n/CE, where n (the 

transmission factor r is 0.05. the effic.ency ,s n - 3-93^ 0 f c, 

l: increase „ provtded 

overheating can be avoided. 
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Fig. 4. Efficiency of IBr solar pumped laser plotted vs length L at different pressures, and concentration 
factors C. The reflectivities are /•, = I, r, - 0.95. respectively. 


An understanding of why the efficiency is low can be obtained if the overall 
efficiency is considered to be the product of four efficiencies ( 3 ): tj = VsVaVkVu- Here 
17,, the solar efficiency depends on the absorption bandwidth of the total solar radi- 
ance. and its value for IBr is 0.12. The absorption efficiency 7 j a depends on the IBr 
gas pressure and depth of absorption d: for low pressures t)„ = cr„( IBrW = 10 ~ 2 pd. 
where p is the pressure in torr. d the depth in cm. The expression is valid up to 
pd = 50; for higher pd most of the photons are absorbed and T) a — * I. The kinetic effi- 
ciency T) k depends on the number of absorption events producing Br* and also on the 
competition between stimulated emission and quenching. Approximately 0.7 Br* 
atoms are produced for every absorbed photon (Table 1). The fraction of Br* atoms 
which yield a stimulated photon can be obtained by examining Eq. 6 , and comparing 
the value of T with the rest of the loss processes, as all eight variables vs time were 
known. The ratio T /(all loss processes) was 0.82, when tj was a maximum, the major 
competitors with T being quenching by IBr and I. The quantum efficiency is r) a = 
0.44 eV/2.5 eV = 0.18. Hence for pd < 50 torr-cm, 17 = 6 x 10 -4 at 5 torr, d = 1 cm, 
agreeing with 5 x 10 -4 in Fig. 4. With pd > 50 torr-cm, r) a — < ► 1 and the maximum 
efficiency would be about 12 x 10~ 2 . 

Pulsed working would further reduce the efficiency because of the duty cycle, 
unless there was adequate storage of energy when the laser was not emitting. 


VIII. POSSIBILITY OF CONTINUOUS LASING 
The computer runs showed that IBr was depleted, while the quenchers I 2 , Br 2 , and 
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I grew with time. For 200 ns pumping pulses, the photon density n had a maximum 

value just prior to the peak input. . , , n _ n to 

The input light pulse was then assumed to be proportional to sin 10 up to 
SO ms and thereafter to be constant. Forp = 3 torr the calculations indicated a pulse 
of 350 ms. and for? = 13 torr the duration was of order 10 ms with C - 5000_Plots of 
the other variables indicated a gradual depletion of IBr and a growth of I. The long 
time scale suggests that a gas flow system would be feasible for steady -state working 
provided the gas temperature was kept within limits. 


IX. TEMPERATURE EFFECTS 

As the absorbed photons have an average energy of several eV, and as the ^ heat 
conduction of IBr is low. very high temperatures can occur. If the pressure p of IBr is 
greater than 50 torr in a vessel l cm deep, then the temperature of the lasing gas can 
be roughly estimated assuming all the photon energy is deposited in the gas: if p <50 
torr then roughly a fraction 10 ~'-p is deposited. The heat produced travels an average 
distance r to the walls of the containing vessel by gaseous conduction. The walls are 
assumed at a temperature 7 U and are thermally connected to a radiator into space. 
Above about 1000 °K. the gas would also act as a blackbody radiator, stimates 
fhow thauf , = 0.5 cm. then for C = 2000 andp = 5 torr the gas temperature would 

bC If f the^ pressure ^s increased (a requirement for high absorption efficiency) still 
higher temperatures are possible, as heat conduction in gases is independent of pres- 
sure Cooling can be enhanced by reducing .v by inserting fingers in the gas. or by 
admitting a noble gas with a high heat conductivity and a low quenching cross 
section. Assuming perfect conduction from the walls to a radiator the same size as 
the collector then the radiator temperature would be somewhat less than 300 K. 
High temperatures can have large effects on quenching, recombination and disso- 

ciation rates, and on exchange reactions. 

To our knowledge experimental data on the effect of temperature on the quench- 
ing of Br* by IBr or I, are not available. The possible analogous case o quenc ing 
I* bv L was reported by Kartazaef er al. (31), whose measurements showed the 
quench ng rates at 1000 °K were 20 times lower than at 300 °K. If a similar depen - 
See occurred for the quenching coefficient of Br* by IBr or L, the high temperatures 

would prove advantageous in this respect. 

Turner and Rapagnani (26) have reported measurements of the recombination 
ratlstol T 7 ’ t- 21,, where* = U * The coefficient ts 1000 tunes 

smaller at 1000 °K than at 300 °K and is analogous to C 6 used above where I recom 

bined with Br and the third body was IBr. . . . e , Rr 

Hish temperatures cause thermal dissociation of IBr wtth the format, on of I. Br_ 
and Br whkh have large cross sections for the quenching of Br The depletion o 
*e lower level by exchange reaettons would also be reduced as the increase ml 
would enhance the reverse exchange reaction. On balance then, heattng would seem 
to be very detrimental to IBr solar laser operation. 


X. CONCLUSIONS 

The time varying solutions of the species rate Eqs. 2 through 9 show that inversion 
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is possible near the start of a pumping pulse, but continuous lasing may not occur. 
The reason is that early in the pulse the Br created from Br* by quenching and 
stimulated emission is sufficiently removed by exchange reactions to maintain in- 
version. Eventually , however, the IBr decreases and the Br density increases to cut 
off the lasing. 

The dominating mechanisms are the source terms, quenching by IBr and I. two 
body exchange reactions and stimulated emission. With 200 ps pumping pulses the 
behavior of the light output can be represented using only these reactions. However, 
if long pulses are to be simulated, then all the reactions must be included. Our 
computer results showed that for 3 torr of IBr, steady depletion occurred and the 
pulse lasted for about 350 /^s with C = 5000. L =1 m. At 13 torr the pulses w ere many 
tens of ms for the same conditions. 

The exchange reactions can overcome the effect of high quenching cross sections, 
and high quenching need not necessarily indicate that a material is unsuitable for 
solar pumped lasing. 

Heating will reduce quenching. It will also cause dissociation and depletion of IBr 
with the formation of I. l>. Br and Br, all of which quench Br*. The exchange 
reaction rates will be reduced causing the population inversion to disappear, thus 
killing the laser. 

If the path length L is small, requiring the solar concentration C to be high for 
threshold to be reached, then the pressure has to be kept low to prevent heating ( <s50 
torr). The overall efficiency tj is then l0 -4 p (torr) or 0. 19£ at 10 torr. Such a 10 kW 
laser would require a collector of 85 x 85 m 2 . If L is large, so that overheating does 
not occur, and p is high enough to absorb all the photons, then tj could reach 1 .2 9? — 
a 10 kW laser collector would be 30 x 30 m 2 . The above efficiencies can be compared 
to the calculated efficiencies of a Br>-CO>-He laser and a C : ,F r I laser, each of which 
was around 0. 1 7c. 

The effective length L can be increased by multiple passes in a flat "box type - ' 
laser. If Z. = 10 m, then lasing should be possible withC — 300. The temperature rise 
would then be 700 °K even if p = 50 torr, with a depth of 0.5 cm, T w = 300 and no 
admixture of a cooling gas. The receiving area of the laser would be 400/300 - 1.4 m 2 
and an effective length L of 10 m could be achieved in about 8 passes. 

The possibility of steady-state working seems feasible if a continuous flow system 
is used. 
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ELECTRONIC ENERGY TRANSFER IN MOLECULAR SODIUM 


(by N.W. Jalufka and W.L. Harries, 
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ABSTRACT 


Emission from Che a1j£ ■ X l Ig band of Na 2 has been observed following 
excitation of the B% u state by the *88 n» line of an Argon ton laser. The 
populating mechanism of the A state appears to be a transfer from B 1 n u to 
the 2 Big state followed by radiative transition to the A 1 !^. The experiment 
was performed with no buffer gas in the sodium cell and at a pressure where 
the time between collisions of Ba(3 S ) and Na 2 long compared to the 

radiative lifetime of Since collision, can not be responsible for the 

Bln u * 2 % transfer the results suggest a different mechanism. We suggest 
that the mechanism is stimulated emission from B 1 n u to the lower vibrational 
levels of 2 x Ig with the stimulating radiation field provided by the black- 
body emission of the cell walls. 
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I . INTRODUCTION 


Electronic energy transfer is known to occur in alkali metal dimers from 
the observation of A-X band fluorescence in both Li2 and Na2 following laser 
excitation to the B state ( 1-3 ) . The recently discovered 2 1 Sg electronic 
energy level of Na 2 (4) has been identified as playing a role in the transfer 
Q) . The transfer of energy from B^*n u to 2^-Eg is generally explained in 
terms of collisionally induced transitions ( 2. 5. 6. 7 ) . While there is no 
doubt that the transfer of energy from B^I^ to 2^Eg and subsequent 
radiative decay to A^-E^ is responsible for the observed A^E^ - X^Eg fluores- 
cence Q) it is not clear that collisions are always responsible for the 
transfer. We have observed A^E^ -+ X^Eg fluorescence following excitation of 
Bln u by the 488.0 nm line of an Argon ion laser under conditions of no buffer 
gas and densities such that the time between collisions for Na(^S) - 
Na 2 (Bln u ) were long compared to the radiative lifetime of Bln u . The 
radiative lifetime of Bln u is known to be » 8 ns (£) so that whatever 
mechanism transfers energy from B^II U to 2*Eg must have a rate comparable to 
the rate of radiative decay of the B^II U state. In our experiment this can 
not be collisional processes. Direct radiative transfer by spontaneous 
emission can also be ruled out due to the small energy separation and the 
dependence which makes the spontaneous transition probability small. 

However, this transition probability is not zero and we suggest that 
radiative transition from B^IIu to 2^Eg by stimulated emission with the 
blackbody radiation of the cell walls providing the radiation field is 
possible. We report in this paper a detailed study of the B^TI U - 2^-Eg 
transfer based on this possibility. 
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II. experimental details 

The experimental arrangement consisted of a stainless steel cell con- 
taining the sodium vapor. The cell was fitted with four viewing ports ar- 
ranged in the shape of a cross. The viewing ports were sealed with sapphire 
windows and a stainless steel cylinder was welded perpendicular to the plane 
of the cross and served as a reservoir. The reservoir and the cross were 
encased in separate stainless steel blocks. Each block could be heated by 
four commercial cartridge heating elements (300 watts each). The temperature 
of each block was controlled separately and temperatures were measured with 
chromel/alumel thermocouples. This cell design allows the reservoir 
maintained at a fixed temperature which controls the total vapor pressure in 
the cell while the cell itself may be maintained at a different temperature 
which controls the Na/Na2 ratio. 

Fluorescence from the excited Na 2 molecules was collected at right 
angles to the laser beam and focused on the entrance slit of a 0.3 meter 
monochromator . The monochromator was equipped with a 1200 line/mm grating 
and had a first order resolution of 0.10 nm. A photomultiplier tube (S-20 
response) was placed at the exit slit of the monochromator. The output of 
the photomultiplier tube was recorded on a strip chart recorder. Calibration 
of the wavelength drive was accomplished by use of a low pressure mercury 

lamp . 

The laser had a power output of 10 mWatts at 488.0 nm which excited the 

B ln u (V - 6, J - 43) level. Spectral scans of the B 1 n u - X x zJ fluorescence 
provided verification of the excitation of b1r u (v - 6, J - 43) as the spec- 
trum was dominated by transitions from this level. 
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III. EXPERIMENTAL RESULTS 


The system was operated with no buffer gas and at a reservoir and cell 
temperature of 723 K. At this temperature the Na atom density is 4.4 x 
10^ cm -3 and the Na 2 molecular density is 2.6 x 10^ cnT^. The 488.0 nm 
line of the argon ion laser corresponds to the X^Eg (v - 3, J - 43) -+ B^n u 
(v - 6, J - 43) transition (£) . The observed fluorescence consisted of the 
B -► X band as well as transitions in the region 600 nm-750 nm. While our 
resolution was inadequate to allow identification of the individual transi- 
tions the overall characteristics were in agreement with the observation 
of Astill et al. (2) who identified the feature as A^-E^ - X^-Eg emission. It 

appears, therefore, that B^-IEu -> 2^-Eg transfer occurs even with no buffer gas 
present and at densities where the time between Na-Na2 collisions is long 
compared to the radiative lifetime of the B^II U state. 

IV. CALCULATION OF COLLISION TIME 

The time between collisions is taken as the reciprocal of the collision 

frequency Z - a n v where a - collision cross-section, n — number density of 
collision particles and v - R.M.S. speed between the particles. The 
collision cross section was taken to be a - *(r)^ where r - 1.5 x 10*^ cm, 
i.e. 1/2 the inter-nuclear spacing. For sodium vapor at 723 K the frequency 
of collisions of an Na2 molecule with Na atoms is r - 1/Z - 3 x 10“ ^ sec. 
This is two orders of magnitude longer than the radiative lifetime of 8 nsec 
for B x n u and therefore collisional rates are not competitive with the 
radiative rates. 
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V. CALCULATIONS OF THE BLACKBODY RADIATION INTENSITY 

We treat the cell as a blackbody cavity. The energy density of radia- 
tion in the wavelength range A to A + dA is given by Planck's radiation 
formula p(A,T) dA - SirhcA'S (exp(hc/AkT-l) ^dA ergs-cm^ where the symbols 
have their usual meaning. The peak of the blackbody radiation field is given 
by Wien's law, TA max - 0.28978 cm deg and for T - 723 K, A max - 4.008 pm. 

The energy difference between B^ (v-6) and 2^(v-0) is 2492 cm'l from the 
energy level values of Kusch and Hessel (2) for and Effantin et al. (4) 

for 2 1 Zg. This energy difference corresponds to a radiation wavelength A * 
4.013 n . The Doppler width at A ■ 4 /■ and T - 723K is AA - 0.835 nm and 
the blackbody radiation intensity in this range AA is p - 3.99 x 10-7 
.3 

ergs* cm 


VI. ESTIMATES OF TRANSITION RATES 

1. Possible Transitions 

The experimental results indicated that absorption of the 488 nm Argon 
ion laser light in sodium occurred from the transition X(v-3) - B(v-6) (path 
#1, Fig. 1). Fluorescence occurred at wavelengths of 730 nm or greater, and 
the emission apparently came from the transition A -* X (path #4, Fig. ) 

The dotted path #5 from B^II U to A 1 ^ is forbidden (u-*u) . However, the 
pathways ^1%, and subsequent path 2^ - A^ are allowed. We 
consider the rate of transfer for these last two transitions, especially 
the transition for " curve- cross ing, " B 1 !^ - 2 1 2 g . 


The transition B^II U - 2% would be competitive only if the stimulated 
emission rates were as great as the spontaneous emission rates B^-X 1 ^ 
We first estimate the ratio R of the stimulated emission rate W (B 1 n u -2 1 E g ) 
with the spontaneous emission rate (B^II u -*X^Sg) . 
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2. Preliminary Estimate of Ratio of 
Stimulated to Spontaneous Rates 

As is well known, the rate of stimulated emission is given by W - pB, 
where p is the energy density of the black body radiation and B is the 
Einstein coefficient. The value of p was obtained from Planck's radiation 
formula and depended on the wavelength A ^ determined by the transition 
from level v^ to level v 2 and the absorption bandwidth. The latter was 
obtained by assuming absorption occurred over a width of order 100 rotational 
lines, so AA was 100 times the Doppler bandwidth of one line. Pressure 
broadening was negligible under the conditions here. 

The Einstein coefficient B is (10) : 


B 


— JT 





( 1 ) 


where iii- t is the matrix element for the electron dipole moment of the 

I l,v 1 ,v 2 l r 

transition B^n u ->2^Sg from level v^ to level v 2 , c a is the dielectric 

constant for free space, and ft is Planck's constant/2ir . 

The matrix elements are of the form |/i| - e J ^ r where 

are the wavefunctions of the two respective states, e is electronic charge, 
and r the length of the dipole. According to the Franck-Condon principle r 
varies only slightly during the transition and we may treat r as a constant 
and take it outside the integral. Then IjJ can be expressed as 



e r 



( 2 ) 


where F- is the Franck-Condon factor for the transition vi ,vo. The 

l,v 1 ,v 2 z 

rate of transitions from the v^ - 6 to a level v 2 in the 2^*2g state is then 


E-8 



(3) 


W 


V 1 V 2 


- x(er) p (A 


V 1’ V 2 


,T) F? /3 £ h 2 

1 , V 1’ V 2 ° 


V 1 “ 6 


The total rate of all transitions from - 6 to all the levels in the 2 1 E g 
state below in energy (corresponding to an upper value of v 2 - v n ) is 


n 


W - l 


w 


- 6 


v 2 -0 


V 1 ,V 2 


w 


The level corresponding to v n in the 2 1 Z g state must lie below the v - 6 
level for the B state for stimulated emission, and its value is v n - 24. 

Next, we consider spontaneous emission, B^IIu-OC^Zg with a rate 
coefficient A given by (10) : 


2 

8 ir 2 , V- , v 

A - — ; - 6 (5) 

3 h v x l 

with the same nomenclature as Eq. (1); |M2| is for the transition B 2 g 

and A a the wavelength corresponding to the transition between the vibrational 
quantum state, B( V1 ) to X(v 3 ) . However, the energy difference is about 2.6 
eV for v X - 6 and A a - 0.5 nm, a value 10 times smaller than the peak 

wavelength for black body radiation. 

Again, in the electronic transition, r is assumed constant, and is the 
same value as in Eq. (3). The total rate of spontaneous emission is then 

( 10 ) 
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( 6 ) 


A 


8 » 2 (er ) 2 
3 ft c • A 3 


OO 

2 




v- - 


6 


where F 2 jV ^,V 3 is the Franck-Condon factor for the transition v^, V 3 . 
Experimentally, the spontaneous lifetime of the B^n u state is known to be 
8 ns, so A - 1.25 x 10® s'^. 

Equations (5) and ( 6 ) can be evaluated if r and the Franck-Condon fac- 
tors are known, but first it is instructive to take the ratio R - W/A - 
pB/A: 



Here, as both transitions started from the B^n u state, r cancels. The 

Franck-Condon factors are finite for only specific values of A . To get 

V 1 » v 2 

a rough order of magnitude value for R, we assume the contributions from the 
Franck-Condon factors in the numerator and denominator are about equal and 
cancel, and define a quantity R' - p (A v ^» A^/4ft. The quantity R' 

essentially reflects the contribution from the radiation density. Assuming 
A a - 0.5 pm and constant, R' is plotted as a continuous function of A 
for different temperatures (Fig. 2). It can be seen that for ^v^,V 2 > ^ M m 
that R' « 10 6 , but if A Vl V2 « 0.5p then R' « 1, indicating that here R 
would be much less than 1. Hence, stimulated emission should predominate at 
around 4pm, but not around 0.5 pm. 
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3. Calculation of Transfer Rates 
We confirmed this conclusion by calculating the stimulated and 
spontaneous rates based on estimates of the Franck-Condon factors. The 
Franck-Condon factors were calculated in two ways: first Morse functions 

(11) were fitted to the potential curves of Jeung (11) ; second the same 
curves were digitized using a computer interpolation technique. Values of w e 
were obtained from Huber and Herzberg (12) and Effatin, et al. (4). Values 
of w e x e were obtained from w e x e - w e 2 /4 D e where D e was read from the 
curves. Values of w e y e were taken as zero. For the second method, the value 
of w e from reference (4) gave poor results and w e was obtained by fitting a 
parabola to the curve near the equilibrium displacement to give a spring 
constant K and estimating w e assuming the known effective mass of the atom. 

The wave functions were calculated by first obtaining asymptotic 
solutions for large r^ the internuclear distance, and then reducing r£ in 
small increments and obtaining iterative solutions. They were well behaved 
for at least 20 levels and satisfied the critera of orthogonality, the 
vibrational sum rule, and displayed a Condon parabola on a v 1( v 2 plane. The 
Franck-Condon factors varied between 0 and 0.2 and agreed with Kusch and 
Hessel (2) in order for the B^-II U -* X^-Sg transition. 

The Franck-Condon factors obtained from the two methods differed by 
factors of up to 5 when comparing the transitions from vj_ - 6 in the upper 
level to individual levels in the lower state. For some values of v 2 the 
Morse curve results were larger, for others it was the reverse. The two sets 
of Franck-Condon factors were then used to calculate the stimulated emission 
rates, where again the answers were different for individual levels. 

However, when the rates were added for 21 lower levels, the differences 
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averaged out and gave total transition rates in good agreement as shown 
below. 


4. Stimulated Emission Rates 

The value of W was obtained using the value of B from Eq. (1) where fix 
was taken as (er) Fi Vl , V2 with F 1 the v « lue for the transition B^n u -*2^Lg: 



( 8 ) 


The value of r was obtained from the known value of A, Eq. (5), where 1^2^ “ 
er F2, and A a - 0.5 then r “ 1.5A, a reasonable value, half the equilibri- 

um intemuclear distance . 

A plot of W vs X , the wavelength of the oven black-body radiation 

v l> v 2 

is shown (Fig. 3), for values of v X - 5, 6, 7, using the Morse curve numbers 
and shows a series of "lines" occurring near A - 5 jim with values * 10 12 s' 1 . 
The addition of all the lines between 0 and 20 fi gave a total rate of several 

times 10 12 for the stimulated emission rate B 1 ^ -* ]Eg. The interpolation 
method gave individual values of W, (v^ — 6) for different V2 in quantitative 
agreement with Fig. 3, but, with individual values which showed differences 
of up to 5. However, when the rates for all the lines were added the net 
rates calculated by the two methods agreed within 20 percent. The spontane- 
ous emission rate from B 1 n u -*X 1 Sg is of order 10®, so the ratio R again 
approaches 10^ . 
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Stimulated emission does not occur from B^-IIu-^X^Zg. The 2 eV difference 
in the electronic levels calls for wavelengths of 0.5 m® where the value of p 
is about six orders of magnitude smaller at T - 450 C. 

In a similar manner the stimulated and spontaneous emission rates were 
calculated for the various pathways shown (Fig. 1) and are shown in Table 1. 
It is evident that the stimulated emission rate B 1 II u - > 2 1 Sg far exceeds any 
stimulated or spontaneous rate B - X. For the A - X transition the stimu- 
lated and spontaneous rates are comparable. 


VII. CONCLUSIONS 

The model we have assumed in the calculations is approximate only. It 
may have inaccuracies in the value of r and also in the estimates of the 
energy levels and hence the actual values of the wavefunctions and Franck - 
Condon factors. Nevertheless, the insensitivity of our estimated total rates 
to variations in parameters such as D e , and the differences of 5 orders of 
magnitude between the estimated stimulated and spontaneous emission rates 
(Table 1) suggest the pathways in Fig. 1 are valid. Moreover the suggested 
mechanism of stimulated emission due to the oven black body radiation seems 
valid for transfer from electronic energy states which overlap, and for 
transfer in "curve crossing" in sodium. 
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list of symbols 


Molecular Electronic Energy Levels 

X^g, A 1 ^, 2^, 

z Collision frequency 

g Collision cross section of molecule 

n Density of molecular sodium vapor 


V 

r 

A, dA, AA 

P 

T 

n 

h 

ft 

c 


VI, v 2 , v n 
^1, v lf V 2 

^2, v lt v 2 


e o 

^1. V»2' ^2* 
F 1.vi,v 2 


R.M.S. speed of molecules 
Internuclear space, or length of dipole 
Wavelength; increments 

Radiation density per unit wavelength 
Absolute temperature In K 
Pi 

Planck's constant 
Planck's constant/211 
Velocity of light 
Vibrational levels In molecule 

Matrix element for transition from vibrational level 

V 1 to v 2 f° r ® ** 

Matrix element for transition from vibrational level 
to v 2 for B -*■ X 

Dielectric constant for free space 
Wave functions 

Franck -Condon factor for tranaitlon fro, vibrational level 
VI to v 2 
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F 


2 


> v l> v 3 


B 

W 

V 1 > v 2 
W 

A 

R 


Franck Condon factor for transition from vibrational level 
V]_ to V3 

Einstein coefficient for stimulated emission 

Stimulated emission rate from vibrational level v^ to level V 2 
Sum of stimulated emission rates for all levels 
Einstein coefficient for spontaneous emission 
Wavelength for transition B(v^) to X (V3) 

Rate of stimulated to spontaneous emission - W/A - p B/A 
Electronic charge 
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Table 1. Transition# r.tas In Na 2 with Initial and final vibrational levels 
average energy difference i(eV) and average wavelength of 


transitions X. 


Initial 

State 

Transition 

Upper vb 
level vi 

Final vb 
level V 2 

A(eV) 

Max 

X 

(a*®) 

i 

i- 1 

>— * 

Spont . 
A(s) * 1 

B 1 n u 

B 1 n u -2 1 Sg 

6 

0-20 

0.3 

4 

8xl0 12 

2xl0 5 


B 1 n u -X 1 2g_ 

6 

0-35 

2.5 

0.5 

* 7xl0 2 

1 . 4x10® 

2 % 

2 1 Sg-A 1 2^ 

0-20 

0-35 

0.6 

2 

2xl0 13 1.6xl0 6 


A^-X 1 ^ 

0-35 

0-35 

1.7 

0.7 

25 4x10^ 

3.7xl0 7 
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Fig. 1 


Fig. 2 


Fig. 3 


Energy level diagram for Na£ according to Jeung Q) , with the 
transition pathways shown. The dotted curves are best fitted Morse 
functions . 


Plots of R' - pX^/bh vs X 

b! v 1» v 2 


for different oven temperatures 


T(°C). The experiment was carried out at 450 # C (723K). 


Plot of transition rates W from vibrational energy levels v^ 

(5,6,7) in B 1 ^ state of Na 2 to vibrational level V 2 in the 2^2g 

state vs X . The bandwidths (« 1.4x10-3 urn) are too narrow to 

vi,V2 

be evident. Oven temperature 723K. 


E-18 











APPENDIX F 


Volume 63. number 6 


OPTICS COMMUNICATIONS 


[ 5 September 


ABSORPTION OF ISOTROPIC RADIATION 

Wy nford L. HARRIES 

Department of Physics. Old Dominion University. Sorfolk. l.'-l 23508. LSA 


Received 4 May 1987 


Estimates of the absorption of isotropic radiation are given for an absorbing medium whose geometrical shape is either a plane- 
parallel sided slab or a cylinder. By including the directional variation of absorption path length, the results show appreciable 
differences to the ‘ unidirectional radiation" formulae. 


The absorption of isotropic radiation such as 
blackbody radiation, is of interest for estimates of 
rates of energy deposition into a medium either as 
heat or molecular vibrational energy. We assume the 
medium is uniform and isotropic and has a char- 
acteristic absorption length A. For narrow band 
absorption of the radiation into molecular vibra- 
tional energy, as in laser pumping [1], k=(na)~', 
where n is the density of absorbing centers and a the 
absorption cross-section which is single valued. The 
fraction of radiation absorbed along a path length / 
in crossing through the medium is then 
F= 1 -exp( — Ilk), and for unidirectional radiation, 

/ can be related easily to the medium geometry. 
However, for isotropic radiation, / varies with direc- 
tion. Here the total fractional absorption for two 
medium geometries, a flat lab and a cylinder, is esti- 
mated and comparisons made with the "unidirec- 
tional formula”. 

The radiation energy enters the medium from the 
black body cavity walls according to Planck’s radia- 
tion formula and the emission and absorption coef- 
ficients of the wall are assumed unity. If the medium 
is contained within a transparent vessel inside the 
cavity, reflections from the vessel surface may occur, 
which are not considered. 

For a plane parallel infinite slab of thickness?/ ( fig. 

1 ) let the energy enter from a small surface dS at PQ. 
A spherical coordinate sytem ( p , a, /?) has origin at 
O. It is convenient to work in terms of angle 
X = nf 2 -a. For isotropic emission the amount of 



P 0 Q 
— — dS — — 


Fig. I . Infinite slab geometry with energy entering through a small 
area d S at PQ. Spherical coordinates ip. a. ft) have origin at O. 
Isotropic radiation entering at angles less than ; is assumed 
completely absorbed within an error e. 

radiation in direction y is proportional to the open- 
ing, namely, dS cos a ( Lambert’s law ) = dS sin y and 
the fraction of solid angle between y and y + dy is 
cosy dy. The path length in the slab is OM = d/sin y 
and thus the fraction absorbed is 

f. d y = /c, sin y cosy 

x [ 1 — exp( —d/k sin y)] dy , (1) 

where is a constant. To integrate numerically we 
assume for y<y\ all the energy is absorbed within 
an error t, because of the long path lengths, i.e.. 
exp( - d/k sin y) = t, where 6 is a small quantity, cho- 
sen as Then 


and the fraction of radiation absorbed for 0 < y < y 

is/,: 


0 030-4018/87/S03.50 © Elsevier Science Publishers B.V. 
( North-Holland Physics Publishing Division) 
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. / 

t\ = k, i sin y cos y dy 


= \k, [d/k In ( 1 c )] : 


(3) 


The fraction absorbed for ; 


v'^y ^n/2 is/;'. 


/-=fc, 


sin y cos y 


x [ 1 — exp( -dlk sin >’)] dy 


( 4 ) 


Integration was performed using Simpson s Ru e. 
with 11 values of y between y and nil. The tota 
fraction absorbed is / 3 =_/l + h and / is a function o 
the ratio kid. so universal curs es can be plotted. The 
normalization constant k x can be evaluated by real- 
izing that ft'-k, sin y cosy dy is the total power 
through dS, or 0.5 The results were insensitive to 
t provided it was small, except where Aid-* 0. For 
j!d4 0-3 it is seen /,> l due to numerical errors. 
Table 1 gives /, versus kJ d and also shown is the 
function F, = 1 — exp( —d/k ). for comparison. The 
raU o /,/F, indicates the error introduced by using the 

“unidirectional formula . 

A. more common geometry for black-body cavity 
experiments is a cylinder of radius r with energy 
entering from the stdes ( fig. 2 ) . Let the cylinder _and 
the small area dS at L be tangential to plane AFGD 
at AD. The amount of radiation at angles 0 and 0 to 
plane AFGD is proportional to dS sin 0 sin 0. For 
section ABCD at angle 0 , path length LM is 
2rsin (9/sin 0. As before, we require exp( -LMM ) -« 
at o’ or 0'(0) = sin _l [2r sin 0/Aln(l/e )], and 
between 0 and 0 + d0 all the radiation for O<0<0' 
is absorbed, consisting of a fraction f 4 (d): 


f J (d) = k 2 I sin0 sin0d0 


— k 2 sin 0(1 -cos0") . 

Between O’ and nil the fraction absorbed is/ 5 (0): 



2 Energy enters through a small area d.S at L. Plane A FGDL 
angential to cylinder. Sect.on ABCD ts at angle d to AF and 
iorptton paths LM are at angle 0 to LD Isotropic rad.auo 
lering this sectton at <S<0 is assumed completely absorbed. 


jt/2 

f 5 (d) = k 2 J sin 0 sin 0 

O 

x[l-exp(-2rsin0/Asin0)]d0, ( 6 ) 

and f b {6) =f*+fi is the fraction absorbed between 0 
and 0 + d0. For Os£0<k/ 2 the fraction absorbed is 

/■ 

n/2 

/ 7= J/ 6 (0)d0. (7) 

0 

The value of k 2 is obtained by observing that the 
total power in through dS for 0. 0 between 0 and n/2 
is K 1 jo 2 !o 2 sin6sin0d6 d0 and K 2 =l. A double 
integration of / 7 over 0 and 0 was performed numer- 
ically using 1 1 values of 0 between 0 (0) and nil for 
each of 1 1 values of 0. Again / 7 is a function of kJr 
(table 1). Comparison with a '‘unidirectional cal- 
cuiation for path length 2r, F ; = [ 1 -exp( - 2r/k )], 


362 


F-3 



Volume 63. number 6 


OPTICS COMMUNICATIONS 


I 5 September 1 98' 


Table I 

Estimates of the fraction of isotropic radiation absorbed versus /Jd for a slab, and versus kjr for a cylinder: is the absorption length, d 

the thickness or the slab and r the radius of the cylinder. The functions J\. t- are the fractions alculated here. C and F : are the results of 
the "unidirectional radiation formulae" The ratios /y'/y and ./-'/ 7 ': are a measure of the errors in using the latter The values of ?\ for 
a d<. 0 . 3 should be ^ 1 . 


Slab geometry Cylinder 




F 

FJF X 

/Jr 

c 

F : 

t- F : 

0.10 

1.042 

1.000 

1.04 

0.10 

l 000 

1.000 

1.00 

o.:o 

1.020 

0.993 

1.03 

0.20 

0.999 

1.000 

1.00 

0.30 

1.002 

0.964 

1.04 

0.30 

0.997 

0.999 

1.00 

0 40 

0.978 

0.918 

1.07 

0.40 

0.994 

0.993 

1.00 

0 50 

0.948 

0.865 

1.10 

0.50 

0.990 

0.982 

1.01 

0.60 

0.916 

0.81 1 

1.13 

0.60 

0.984 

0.964 

1.02 

o.?o 

0.887 

0.760 

1.16 

0.70 

0.978 

0.943 

1.04 

0.80 

0.848 

0.713 

1.19 

0.80 

0.970 

0.918 

1.06 

0.90 

0.816 

0.671 

1.22 

0.90 

0.962 

0.892 

1.08 

1.00 

0.785 

0.532 

1.24 

1.00 

0.952 

0.865 

1.10 

1. 10 

0.756 

0.597 

1.27 

1. 10 

0.942 

0.838 

1.12 

1.20 

0.728 

0.565 

1.29 

1.20 

0.931 

0.81 l 

1.15 

1.30 

0.702 

0.537 

1.31 

1.30 

0.919 

0.785 

1.17 

1.40 

0.678 

0.510 

1.33 

1.40 

0.908 

0.760 

1.19 

1.50 

0.655 

0.487 

1.35 

1.50 

0.895 

0.736 

1.22 

1.60 

0.633 

0.465 

1.36 

1.60 

0.883 

0.714 

1.24 

1.70 

0.613 

0.445 

1.38 

1.70 

0.871 

0.692 

1.26 

1 80 

0.594 

0.426 

1.39 

1.80 

0.851 

0.671 

1.28 

1.90 

0.576 

0.409 

1.41 

1.90 

0.846 

0.65) 

1.30 

t j 

o 

o 

0.559 

0.393 

1.42 

2.00 

0.833 

0.632 

1.32 

3,00 

0.432 

0.283 

1.52 

3.00 

0.720 

0.487 

1.45 

4.00 

0.352 

0.221 

1.59 

4.00 

0.628 

0.394 

1.60 

5.00 

0.297 

0.181 

1.64 

5.00 

0.556 

0.330 

1.69 

6.00 

0.257 

0.154 

1.68 

6.00 

0.498 

0.284 

1.76 

7.00 

0.227 

0.133 

1.70 

7.00 

0.450 

0.249 

1.81 

8.00 

0.203 

0.118 

1.73 

8.00 

0.411 

0.221 

1.86 

9.00 

0.184 

0.105 

1.75 

9.00 

0.378 

0.199 

1.90 

10.00 

0.168 

0.095 

1.76 

10.00 

0.350 

0.181 

1.93 

1 1. 00 

0.154 

0.087 

1.78 

11.00 

0.326 

0.166 

1.96 

12.00 

0.143 

0.080 

1.79 

12.00 

0.305 

0.154 

1.99 

13.00 

0.133 

0.074 

1.80 

13.00 

0.287 

1.143 

2.01 

14.00 

0.124 

0.069 

1.81 

14.00 

0.270 

0.133 

2.03 

15.00 

0.1 17 

0.064 

1.81 

15.00 

0.256 

0.125 

2.05 

16.00 

0.1 10 

0.061 

1.82 

16.00 

0.243 

0.118 

2.07 

17.00 

0.104 

0.057 

1.83 

17.00 

0.231 

o.i n 

2.08 

18.00 

0.099 

0.054 

1.83 

18.00 

0.220 

0.105 

2.09 

19.00 

0.094 

0.051 

1.84 

19.00 

0.21 1 

0.010 

2.11 

20.00 

0.090 

0.049 

1.84 

20.00 

0.202 

0.095 

2.12 
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shows /? can be considerably greater than Fz. 

In conclusion, estimates of the absorption of iso- 
tropic radiation should take account ot the contri- 
bution of rays proceeding at angles which result in 
very long path lengths. The above procedure can be 
extended for broad band absorption if er is a known 
function of wavelength b.y reiterating over small 
wavelength ranges. 
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THEORETICAL STUDIES OF A SOLAR-PUMPED BLACK-BODY C0 2 LASER 

By 

Wynford L. Harries and Zeng-Shevan Fong 


SUMMARY 

A method of pumping a CO 2 laser by a hot cavity has been demonstrat- 
sd (ref. 1). The cavity, heated by solar radiation, should increase the 
efficiency of solar pumped lasers used for energy conversion. Kinetic mod- 
eling is used to examine the behavior of such a C0 2 laser. The kinetic 
equations are solved numerically vs. time and, in addition, steady state 
solutions are obtained analytically. The effect of gas heating filling the 
lower laser level is included. The output power and laser efficiency are 
obtained as functions of black body temperature and gas ratios (C0 2 -He-Ar) 
and pressures. The values are compared with experimental results. 

1. INTRODUCTION 

The possibility of using solar collectors on orbiting space vehicles 
and transmitting the energy to earth via laser beams has been discussed 
already (refs. 2, 3, 4). A critical factor for energy conversion is the 
laser efficiency. Oie suggestion for increasing the efficiency is to heat a 
black body cavity by solar radiation and then let the cavity pump the laser, 
according to the Planck radiation law. If the laser mediun absorbed a nar- 
row bandwidth from the spectrum then the black body cavity would immediately 
refill any "hole" in the distribution. The pimping radiation is thus 
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replenished and the potential efficiency should be much higher than if the 
laser is pumped directly. 

It is the purpose of this paper to examine the kinetics of a black body 
pumped laser. In addition to directly pumped C0 2 lasers, there are schemes 
whereby N 2 is pumped and the energy transferred to C0 2 , but these are not 
considered here. The direct C0 2 lasers usually have He and Ar as additives 
which perform the functions of deexcitation of the lower level and heat 
conduction; these are included. 

The physical mechanisms of the laser are discussed in Section 2. For 
the first time to our knowledge the effect of gas heating, which serves to 
fill the lower laser level and inhibit lasing, is included. The temperature 
rise is calculated, and an estimate of the effect of filling the lower laser 

level from Boltzmann statistics is made. 

In Section 3, time varying solutions of the kinetic equations show the 
behavior of the upper and lower laser level densities, and also the photon 
density and output power vs. time. The solutions show that steady state 
conditions can be achieved. 

Analytical solutions are obtained for the steady state in Section 4. 

The steady state output power is shown as a function of black body temper 
ature and gas constituent pressures. Estimates of the laser efficiency as a 

function of the variables are also given. 

Throughout the discussion it is assumed that the laser tube of 0.4 cm 
radius has a length of 50 cm heated in an oven at 1000-1500 K and contains 
COj +He+Ar at total pressures of 1 to 12 torr in varying ratios. The radi- 
ation from the oven is passed through nitrogen which is transparent to it 
and enters the laser tube. The purpose of the nitrogen, which is flowed, is 
to cool the laser medium. The measured exit temperature of the nitrogen was 
340 K (ref. 1), and this was assumed to be the wall temperature. The gas 


temperature is slightly higher than tMf and „ d1scussed ^ ^ ? 

The above ambers correspond to experimenta, conditions in reference 1 . ‘ 

2. PHYSICAL MECHANISMS IN THE LASER 

A flow diagram of the processes occurring during lasing is shown in 
figure 1. 

2.1 Absorption of Radiation 

The black body radiation is absorbed at 4.256u over a bandwidth of 0 1 

"■ "" radiati ° n flUX P6r unit ,en 9th entering the laser tube of radius r 

if EC*) d* *r When E(*) is determined by Planck's radiation formula and dk 

' S <U ThS fraCt1 ° n absorbed depends on n a> the absorption cross sec- 
tion. The latter can be calculated: . kWav, „ here A . 424 . 6 „ t|)( 

Einstein coefficient for the transition (OOO).(OOl), k = 1 0 .6„, 4A . 0 . 1 „ 

end o a = 1.2xl0-*«cm*. However, a measured value of the absorption length 

15 1 - 56xl0 ' 3c "’ at 1 Atmosphere (ref. 5), which yields o a . 2.41*10-i»c 
Both values were tried. 

TWO possibilities exist for the calculation of the fraction absorbed 
If the absorption is small, then the fraction absorbed per unit volume can 
be expressed approximately as 


F i = 2 (1 - exp (-(C0 2 ) a ?r)/r 


( 1 ) 


-here the radiation is assumed to cross a distance 2r and the amount per 
omt length is proportional to 2mr and per unit volume, Uvr/wr*. ,f the 
pressure increases, however, the absorption will take place at the outer 
edges of the tube. The excited species can, however, diffuse to the center 
in a time t where t - r^/D and 0 is the diffusion coefficient. Assuming a 
mean frM ° f et 1 torr, a total gas pressure of 1 torr 
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Figure 1. Flow Diagram for C0 2 Black Body Laser 



(ref. 6) and a temperature of 340*, then t - 3 ms, a time comparable with 
other rates. Neglecting diffusion, a second estimate is obtained by 
assuming the radiation is absorbed at the center of the tube between a depth 
x i and x 2 . Then the fraction F absorbed per unit volume is 

F 2 = (exp(-(C0 2 )a a x 1 ) - exp(-(C0 2 )o a x 2 ))/(x 2 -x 1 ) ( 2 ) 

where here Xj = 0.35, x 2 = 0.45 cm. The effect of the outer layer of gas 
shielding the center of the tube at high pressures is now included, but a 
one dimensional, not cylindrical, geometry is assumed. 

A number of computer runs were made using F 2 , however, the results did 
not agree with the experiment (ref. 1). The computer results showed a sharp 

cutoff of power at a few torr and, accordingly, the coefficient Fj was used 
for the following calculations. 

The nimber of absorption events per unit volume is then 
E(A)dA 

S = F i = 5.76 x 10 19 F 1 /(exp(3381/T b )-l) (3) 

where T b is the black body temperature, which gives the number of C0 2 (000) 
molecules raised to the upper laser level C0 2 (001) per cc per second. 

2.2 Heating of the Medium 

Exper imentally it was realized that unless a nitrogen coolant was flow- 
ed past the laser tube, that lasing ceased (ref. 1). The implication was 
that the medium became heated causing the lower level to be filled. A rough 
estimate of the increase in temperature, and its effect on the density of 
the lower level follows. 
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As the efficiency of the laser would probably not be greater than a few 
percent, it is assumed that all the radiation absorbed goes into heating the 
gases. The rate at which heat energy enters a unit volune, dQ/dt, increases 
with the pressure of CO 2 which is included in f\, 2 > 

dQ/dt = E(X ) dx • Fi/4.2 Cals cnrV 1 (4) 

The heat conduction K of a mixture of gases CO 2 , He, Ar depends on the 
gas ratios and the absolute temperature (ref. 3). It is independent of gas 
pressure. 


K = 3.27 x 10- 4 S VT/273 (5) 

/I + 0.71 C0 2 /He + 0.32 Ar/He\ 

where S = ( — I (6) 

\l + 4.55 C0 2 /He + 2.67 Ar/He/ 

It is assumed that the heat is conducted over an effective area A 

through an effective distance d to the wall at a temperature T . As the 

w 

conduction is proportional to vr. the conduction equation yields that 
the temperature in the center of the gas T g is: 

T g * (Tw + 1.8 x 10 4 • (dQ/dt)/ (S ♦ A/d)) 273 (7) 

The expression is approximate as a one dimensional geometry is assuned. For 
example, with T w * 340 K,^ a black body temperature of 1500 K, a total 
pressure of 10 torr (CO^ 20%, He 20%, Ar 60%), tube radius 0.4 cm, then T^ 
approaches 360 K. 
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2.3 Filling of the Lower Laser Level 
The lower laser level which has a density Nj of C0 2 (100) atoms, will be 
filled according to Boltzmann statistics which depend on T However com- 
peting against this will be the depletion processes due to collisions. To 
include the temperature effects they must be expressed as a rate of produc- 
tion of N 1 vs. time. It is assumed that C0 2 molecules in the tail of the 
energy distribution, which have energies greater than = 0.1649eV, the 
energy of the lower laser level, can collide with the rest of the C0 2 
population with a finite probability of producing C0 2 (100). The density of 
molecules with sufficient energy is (CO*) exp ( -E ] /kT) and the rate of 
production from this process is: 


* (C0 2 ) (C0 2 * exp (-Ej/kTg) <*v (8) 

Here k is Boltzmann's constant, a is the cross section for the reaction, v 
the particle velocity. R 1 increases as T g increases. 

2.4 Diffusion Losses 

At low pressures losses of particles by diffusion to the tube walls may 

be important. We assume that the diffusion loss time of C0 2 molecules, 
obeys (ref. 7): 


t D = ( r /2.405) 2 /Dr (9) 

where D is the diffusion coefficient for C0 2 in a mixture of C0 2 , He and Ar, 
(Appendix I) and the 2.405 corresponds to losses in an infinitely long cyl- 



inder. The quantity y is the probabi 1 ity of deexcitation at the wall and is 
somewhere between 0.2 and 1. The diffusion losses apply equally to all 
excited level s of CO2 . 

2.5 The Kinetic Equations 

The CO2 laser (Fig. 1) is a four level system. The absorbed photon 
raises the CO2 (000) level to the upper laser level C0 2 (001) the asynmetric 

stretching mode at 0.2912 eV above ground. Lasing at 10. 6u occurs by a 

transition between the (001) and (100) levels, the latter a synmetric 

stretching mode at 0.1649 eV above ground. 

The upper laser level of density N y is depopulated by collisions with 
COj, He, Ar (rate constants K x , K 2 , K 3 ) to the ground level, and by 

spontaneous emission (Einstein coefficient A). Transitions can also occur 
to the lower laser level (100) (K 7 , Kg, K 12 , K 13 , and the spontaneous 
coefficient A^-j) and the (010) level (K a , Kg, K 14 , K 15 ). 

The (100) level (density N-j ) has a (020) bending mode which is, 

0.0005 eV higher in energy. The conversion (020)^(100) is very rapid, and 
the two are lumped together here. 

The lower laser level is filled by collisions from the upper level (K 7 , 

Kg, K12 , Ki 3 ) and by the spontaneous emission A ul . It is also filled during 
lasing by stimulated transitions of rate constant r. The lower level is 

depleted by collisions with CO2 (K10, Ku) and with He (K^), (no data on Ar 
was available), and the transitions end either on the (010) level or the 
(000). However, at high gas temperatures, the lower laser level will be 
filled at a rate R-j . 

The (010) level of density N m is in turn depleted by collision with 
C02, He, Ar (K 4 , K 5 K 6 ). The reactions are given in more detail in Table 1, 
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Table 1. Rate Coefficients for C 0 2 at 360*K 
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with values from four different sources. In the colunn under reference 9, 
the values are calculated from the expression 

- 1/3 - 2/3 

In K = A + BT + CT (10) 

where the constants A, B, C, are given for the different reactions. The 
values may not be accurate when extrapolated to 360K. References 8, 10 and 
11 give experimental values. 

The unknown quantities are then assumed to be N u , N-j , N m and N the 
10. 6u photon density. In the following equations, the quantities in paren- 
thesis (CC^) etc. are the densities of the species. The rate equations for 
the four species are: 

dN u 

= S - N [(Kj + K7 + Kg + Kg ) (CO2 ) + (K 2 + K12 + Ki4)(He) 

dt u L 

+ (K 3 + K 13 + Kis)ArJ - N^(A + A^)- r - u/ 


dN 1 

r + N u A ul + N u [ (l< 7 + <9)(C0 2 ) + K 12 (He) + K 13 (Ar)] 

dt 

- N«| + ^iiHCOj) + K 16 (He) ] - e ^ x Q + 

dN m 

— - N ] (K 10 + 2K n )(C0 2 ) + N u [(K 8 +K 9 ) (C0 2 ) + K 14 (He)+K 15 (Ar)] 

- N m [k 4 (C0 2 ) + K s (He) + Kg ( Ar ) ] - N ffl/ (13) 

t D 


G-ll 


dN 


( 14 ) 


The stimulated emission rate r is given by: 

r = N a e c [M u - N,] (15) 

and is proportional to the stimulated emission cross section whose value 

was calculated: a g = A ul / 4* 2 c AA 0 , where X Q = 10. 6u, A yl = 0.186 and 

4A g was taken as the bandwidth from Doppler Broadening. Here c is the 

velocity of light. Assuming a gas temperature of 360’, then AA =0 217A 

e 

and = 9.14xl0- 17 cm 2 . 

In equation 15 the quantity 6 is a geometrical factor which takes 

into account that the spontaneous emission is isotropic; G = Zr z /l 2 where 

r and L are the radius and length of the laser. The quantity t c is the 

containment time for photons, T( . = -2L/c In ( ri r 2 ) where r x = 0.99, r 2 = 

0.95 are the reflectivities of the laser mirrors. 

The equations can be simplified by noting that N , N. , N «(C0->) and 

that (CO 2 ) can be regarded as constant. Then the following constants can be 
defined: 


a l = [ (Ki+K 7 +Ke+K 9 )(C02) + (K 2 + K 12 + Km) (He) + (K 3 + K i3 + Ki 5 )(Ar) 

+ A + + 1/T o ] 
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1 

a 2 » [ (Kio+Kil ) (C0 2 ) +K 16 (He) + — )1 

T D 

A' = [ (K 7 +K 9 ) (C0 2 ) + K u (He) + Kia(Ar) + A u] ] 

a 3 = [K u (C0 2 ) + K 5 (He) + K 6 (Ar)] 

a,, = [(K 10 + 2 Ki 1 ) (C0 2 ) ] 

a 5 = [(K 0 +K 9 )(CO 2 ) + Km (He) + K 15 (Ar)] 

Then equations 11 through 14 can be reduced to 


dN 

u 

dt 

S 

- ^ 1 % 

r 

dN 1 

dt 

r 

+ 

- a 2 N-| + R 1 

dN 

rn 

dt 

*5% 

+ &4 N-j 

- a 3 N m 

dN 

dt 

r 

+ G A ul 

"u‘ "/'c 


(16) 


(17) 


(18) 


(19) 
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3. TIME VARYING SOLUTIONS 

Computer graphic solutions of equations 15 through 19 are shown in 
Figure 2. It was assumed for convenience that the pumping signal from the 
black body rose from zero to full value in a time of 100 us, and that the 
signal was proportional to sin 2 2m (t/10- 4 ) (where t is the time in sec- 
onds) up to 100 us; thereafter it was constant. The pumping signal is shown 
in the upper curve. 

In the second plot, N^, N^ and (N^-N-jJxlO are plotted. The densities 
of the upper and lower laser levels increase with time, but at first, when 
the signal is small, N-| > N^. (This can.be seen on amplified plots). How- 
ever, N u overtakes N ] at 28us, and the signal (N (J -N ] ) starts to grow. It 
increases with time until threshold at t = 53 us. Then N u falls and N-j 
rises, and the inversion -N i ) falls. At the same time the photon density 
N increases sharply and then falls to zero as the inversion falls. The 
upper level is then increased again until a second light pulse appears and 
so on. This behavior is characteristic of kinetic equations similar to 
these, where is small; if o g is much larger these violet oscillations 
are smoothed out. 

The inversion quickly assumes a constant value, and at 135 us the value 
of the photon density N also becomes constant. This implies that continu- 
ous lasing is possible. 

4. STEADY STATE SOLUTIONS 

As they stand, equations 16 through 19 are nonlinear. In the steady 

state, however, the stimulated emission rate far exceeds the spontaneous 

rate, and in equation 19, r >> GA ,N or: 

ul u 
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o 100 160 


t(us) 

Figure 2. Time behavior of black body pumped laser. The pumping signal 
was assumed proportional to sin 2 (Tit/2xl0 -4 ) up to lOOus, then 
became constant. The middle diagram shows Nu» N] and (Nu-Ni)xlO. 
The lower diagram shows the 10. 6y photon. 



r 


( 20 ) 



c 


and from 15 and 19 


N u - N 1 


* const 


CT 

e c 


( 21 ) 


The constancy of the inversion N (j -N ] is evident in Figure 2, and the 
value is given by equation 21. Replacing r in equations 16 through 19 by 

N/t the solutions are 
c 


S + R 1 + a 2 /T c a e c 
aj + a 2 - A 1 

t 


N, = N - 
1 u 


Ve c 


( 22 ) 


(23) 


N 


m 


a 5 N u + a i+N-| 


a 3 


a l 

1 - 


a l a 2 


- T c R, a , - 

ai+32-A' 


(T e C( a i + a 2 - A ' ) 


ai+a 2 -A' 


(24) 


(25) 


The effect of gas temperature is included in R 1 and as R ] increases the 
value of N falls. 


The value of N from equation 25 agrees with the value, late in time, 
from the time varying computer solutions. 
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4.1 Steady State Power Output 

The evaluation of N enables the power output of the laser to be 
evaluated: 


1 

p = _ Nc hv T nr 2 (26) 

2 r 

The factor 1/2 takes into account that half the photons are travelling away 
from the output mirror of transmission T * (l-r 2 ), c is the velocity of 

light, hv the energy of the lasing photons of 1.868xlO- 20 0. Experiment- 
ally it was found that the output beam had a radius of 0.3 cm corresponding 

to a Gaussian profile in a 4 cm radius tube, (ref. 13) and for T r = 0.05, 

P = 3.97x10- 12 N (W) (27) 

Plots of power output vs. total pressure p (torr) for different black 
body temperatures are given in Figure 3. The gases were in the ratio C0 2 
16*, He 4*, A r 80*. The output power rises with black body temperature, but 
as p increases, the power falls. The effect is due to the R.| term in the 
equations; on making R i = 0 the output power rose with increasing pressure. 

The variation of output power with black body temperature at a fixed 
total pressure of 7 torr, but with the same gas ratios, is shown in Figure 
4. A linear rise is observed. 

The gas temperature of T^ can be altered by assuming different wall 
temperatures T . Assuming a gas mixture as before of 7 torr (16* C0 2 , 4* 
He, 80* Ar) and a black body temperature of 1450K, then the power output 
vs. T is plotted in Figure 5. It shows that an increase in gas 
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Figure 3. Power output of black body pumped laser vs total pressure 
for various black body temperatures. The gases were 
C0 2 -16%, He-4%, Ar-80°. The wall temperature was 340K 









Figure 6. Efficiency (as defined in equation 29) of the black body 
pumped laser vs total pressure p for various black body 
temperatures. The gases are CO 2 16%, He 4% and Ar 80°. 




temperature greatly reduces the output power indicating that cooling of the 
mediun is essential. 


4.2 Efficiency 

The input power P in is the power absorbed by the laser: 

P in * dA 2irr L (l-exp(-(C0 2 )<^ a 2r) (28) 

and the efficiency n is defined here as output power/the power absorbed 
from the black body radiation: 


n = P/P. n (29) 

(It should be noted that the black body pumped laser has other "efficien- 
cies" besides n, however, these are not considered here.) 

Using equation 29 to define the efficiency. Figure 6 shows plots of n 
vs. p for different black body temperatures. There seem to be high effi- 
ciencies at very low pressures where the output power is small. In the 
regime around 7 torr where experimentally the highest output powers were 
observed, the efficiency is only a few percent. 

CONCLUSIONS 

The results here depend critically on the values of the rate coeffi- 
cients in Table 1. They are insensitive to K 4 , K 5 , K 6 , but are very de- 
pendent on K 9 , K n and *13- The most sensitive parameter seems to be 1 

where a change of a few percent altered the output power by several times. 
Also Tg was a sensitive parameter. 

Unfortunately, the rate coefficient for R-j was not known. The values 
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chosen indicate that about one collision In 10 s of two ground state atoms 
resulted in the creation of one CO 2 (100) which does not seem unreasonable. 

Figures 3, 4, and 5 show reasonable qualitative agreement with 
experimental measurements (ref. 1). However, it did not seem possible to 
get agreement with experimental measurements of power output vs. total 
pressure (ref. 12), at a black body temperature of 1450K where the ratios 
C^iHerAr were varied. In view of this, these results must be treated with 
caution, and must be regarded as preliminary. 
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APPENDIX 


Diffusion of CO 2 When Mixed with He and Ar 


The mean free paths of molecules in their own gas at 1 torr for CO 2 , He 
and Ar are 2.7xl0- 3 cm, 13.1xl0 _3 cm and 4.6xl0~ 3 cm, respectively (ref. 7). 
These values give the cross sectional diameters, and it is then assumed that 
for a collision between unlike molecules that the cross section diameter is 
the arithmetic mean of the individual diameters. Adding the collision 
frequencies of C02-C02> C0 2 -He, C0 2 -Ar, the total collison frequency v for 
CO2 is: 


v 


/ He 

2.37x10- 11 (C0 2 ) 1 + 1.79 + 

' C0 2 



(Al) 


The diffusion coefficient D is then 


D = 


1 v 
3 


1. 928x1 0 6 T/v 


where 


v r = 2405 


i; 


the random velocity of CO 2 . 


G-25 




